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FOREWORD 


Even though Appendix A contains a list of symbols used in this report, additional explana- 
tion is provided. In section 2, d/dt represents the time derivative, but in section 3, in the 
discussion of finite element modelling, 'x' is used. For the scalar product is used and 
for the vector product, *x\ Vector quantities are underlined except in the finite element 
formulation. Some symbols are not listed in Appendix A, but these are used only locally in 
specific subsections and are not referenced again in any other subsection. 
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SUMMARY 


This report is intended to serve two purposes. The first is to give a survey of the theoreti- 
cal background of the dynamic interaction between a non-viscid, compressible fluid and an 
elastic structure. Section 1 presents a short survey of the application of finite element 
method (FEM) to the area of fluid-structure-interaction (FSI). Section 2 describes the 
mathematical foundation of the structure and fluid with special emphasis on the fluid. The 
main steps in establishing the finite element (FE) equations for the fluid structure coupling 
are discussed in section 3. 

The second purpose of this report is to demonstrate the application of MSC/NASTRAN to 
the solution of FSI problems. Some specific topics, such as fluid structure analogy, 
acoustic absorption, and acoustic contribution analysis are described in section 4. Section 
5 deals with the organization of the Acoustic Procedure implemented in MSC/NASTRAN. 
Steps which have to be performed for a complete acoustic analysis are illustrated by a 
flowchart. Section 6 includes the most important information that a user needs for applying 
the Acoustic Procedure to practical FSI problems. Beginning with some rules concerning 
the FE modelling of the coupled system, the NASTRAN USER DECKs for the different 
steps are described. The goal of section 7 is to demonstrate the use of the Acoustic Proce- 
dure with some examples. This demonstration includes an analytic verification of selected 
FE results. The analytical description considers only some aspects of FSI and is not in- 
tended to be mathematically complete. 

Finally, section 8 presents an application of the Acoustic Procedure to vehicle interior 
acoustic analysis with selected results. 


1. INTRODUCTION 

In many areas of everyday engineering the interaction between structure and fluid plays an 
important role. Phenomena of this kind have to be considered if a structure surrounds or is 
embedded within a fluid. The presence of the fluid can significantly alter the behavior of a 
structure, and at the same time, the deformation of the structure alters the loads transmitted 
from the fluid. This means that there exists an interaction in both directions: the structure 
acts on the fluid and vice versa. These phenomena are referred to as fluid-structure interac- 
tion (FSI). 

These phenomena can be divided in subareas, such as: slow viscid flow, laminar flow, tur- 
bulent flow, shallow water flow, meteorology, acoustic phenomena and FSI. 

In a 1982 paper of O.C. Zienkiewicz (ref. 1) one can find a more detailed discussion of 
these subareas. In that paper, research progress and the advantages and drawbacks of the 
FEM for the solution of FSI problems are described. An overview on the use of FEM in 
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fluid dynamics and related areas is presented in the Finite Elements in Fluids proceedings 
(ref. 2). 

The present description deals only with the last two areas listed above, where the fluid is 
considered to be compressible and inviscid. 

There are many problems in everyday engineering where fluid structure interaction is pre- 
sent. Sound radiation into passenger compartments of automobiles and airplanes is a typi- 
cal example. The main purpose of acoustic calculations is to analyze the sound transmis- 
sion into the passenger compartment and the coupling mechanisms between sound con- 
ducted through solids and sound propagation in fluids (e.g., air). Through the use of these 
techniques one can potentially influence the design of an automobile or an airplane by iden- 
tifying structural components that produce high sound pressure levels in the passenger 
compartment. 

The calculation of sound pressure distribution in an interior domain with complex shaped 
boundaries (a passenger compartment of an automobile) is not possible through the use of 
empirical formulas. Instead, numerical solutions of initial and boundary value problems 
for various interior domains with different types of boundaries (acoustic soft and rigid 
boundaries, elastic and absorptive boundaries) are used. Although there are some restric- 
tions, FEM is a very useful and powerful tool for solving these kinds of problems. 

This report presents a unique capability to perform FSI analysis using the finite element 
code MSC/NASTRAN. Through the use of this added capability, MSC/NASTRAN be- 
comes a versatile tool for acoustic and noise control analysis. 

The first work in the acoustic field was done by Gladwell and Zimmermann who presented 
two papers in the mid 1960s (see ref. 3 and 4). These papers include finite element (FE) 
formulations for coupled fluid-structure vibrations. Some years after these basic develop- 
ments, Craggs presented, in a series of papers dated 1970-1973 (see refs. 5,6, and 7), a 
variational formulation for the coupled problem in terms of fluid pressure. 

Very early it was realized that the commercial finite element codes, originally developed for 
analyzing structural mechanics problems, could be used to analyze scalar field problems as 
well (see ref. 8). With FEM it is possible to investigate standing waves and the related fre- 
quencies in any complex-shaped interior domain with rigid boundaries. In order to solve 
eigenvalue and response problems for the coupled fluid-structure systems, it is necessary to 
consider the elastic as well as the absorptive boundary conditions. This task is accom- 
plished, in this procedure, through the use of the NASTRAN DMAP language (Direct Ma- 
trix Abstraction Program) and a single Fortran program. No changes to the source code of 
NASTRAN were necessary. 
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A description of this approach for solving the FSI problems can be found, e.g., in refer- 
ences 9, 10, and 11. Specifically, for the applications to the interior acoustics of automo- 
biles, refer to references 12, 13, 14 and 15. 

Besides the interaction of coupled fluid-structure systems via the elastic boundary condi- 
tion, sound absorption effects at the boundary also have to be considered. In the present 
procedure, a simple normal impedance model was developed. This model is, in general, 
sufficient for describing the absorption effects at the boundary in the low frequency range. 

Using a special boundary condition of D.B. Bliss (ref. 16), an oblique absorption damping 
can be described. An application can be found in reference 17. Also, A. Craggs published 
some papers which dealt with the question of boundary absorption and the description of 
such effects using the FE technique (see refs. 17, 18, and 19). 

The procedure presented here can be embedded in the standard procedure of structural 
analysis. The fluid medium is treated simply as another part of the structure. This allows 
for a simultaneous solution of the two systems, without the need for separate analysis, as is 
presently done in some cases. 

The procedure requires the use of superelement solution sequences, but it is not necessary 
to define superelements. In many cases, though, it is highly desirable to use superelements 
because of computational efficiency. 

The complete procedure consists of five steps beginning with the uncoupled modal analysis 
of the structure and the fluid (with rigid boundaries). Following the eigenanalysis of the 
coupled system, response calculations in the time or frequency domain can be carried out. 
Furthermore, it is possible to perform an acoustic contribution analysis identifying the in- 
fluence of the structure on the sound pressure level at an interior fluid point. 

It should be noted that this capability is particularly applicable to the interior problems; that 
is, for problems in which the fluid is contained by the structure. The exterior problem, 
fluid outside the structure (e.g., submarines) is very difficult to solve with the FEM. In 
this case, the boundary element method is more applicable. 
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2. MATHEMATICAL DESCRIPTION 


In this section, a mathematical description and a mathematical derivation of the FSI phe- 
nomenon are given. Because the dynamics equations for elastic media are well known (one 
can find a derivation in any book on structural or continuum mechanics), the main effort 
will be directed at deriving and describing the mathematics of the fluid and the FSI. 


2.1 Structure 

The structural components can be formulated as 


UUs) = E(r,t) ( 1 ) 

where L is a partial differential operator with respect to time and space. The vector, u s , 
represents the structural displacement vector, which is a function of location r and time t. 
Finally, £, is a time-dependent load on the structure, which can be divided into an external 
time-dependent load F s (r,t) and the influence of the fluid pressure on the structure 
F p (r,t). That load will manifest itself through the surface load vector 

F p (r,t) = p(r,t)n (2) 

where p is the unsteady fluid pressure and n is the outward unit normal vector at the fluid- 
structure interface. 

Note: The subscript 's' denotes structure and f the fluid. 


2.2 Fluid 

The governing equations for the fluid can be obtained from most books on fluid dynamics 
(see ref. 19), but it is useful to offer a complete derivation of these equations in order to get 
a more complete understanding of the problem. Some of the following mathematical ex- 
planations were initially given in reference 20. 

For an arbitrary fluid volume, V(t), the mass balance is given by 

— ff pdV W qdV (3) 

dt\ 3v C > ) M‘) 

where p is the density and q is the added fluid mass per unit volume and time. Reversing 
the order of differentiation and integration yields 

<4> 
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where v is the fluid velocity field. 


Because the volume, V(t), is arbitrary, 

&+pVv = q (5) 

Although, when the system matrices are defined later, a Raleigh-type damping will be taken 
into consideration, for now the viscosity will be neglected. Internal forces acting upon a 
volume, V(t), in the fluid act at the boundary, S(t), of that volume. Furthermore, they are 
perpendicular to the boundary, which is a consequence of the inviscid fluid assumption. 
External forces are the body forces, b, (per mass unit). Conservation of linear momentum 
yields 

-r\ pvdV = -{ pndS + f pbdV + f qvdV (6) 

dt iv ( 0 Js(o M<) Mir - * v 

where v 1 is the velocity of the added fluid mass. 

Interchanging the order of space integration and differentiation with respect to time on the 
left-hand side and applying the Gauss theorem to the first term on the right-hand side of 
equation (6), one can obtain the following equation: 

Jv (l ) + ^ V • = + P b + W q )d y 

As before, the volume will again be assumed to be arbitrary, which leads to 

j t {pY) + P^Y) = -^P + pb + qv q (7) 

Using equation (4), 

d iz \ dp * dv . * dv 

^(pv) = -v + p- = «-p(V.!;)v + p- 

Combining this result with equation (7) results in 

p^ = -Vp + pb-q(v-v q ) (8) 

Equations (4) and (7) are sufficient if the fluid is an incompressible fluid and p is a known 
and time-invariant quantity. These equations can then be solved with respect to p and y. If 
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the fluid is compressible, then one more relation containing p and p is needed. Assuming 
the fluid motion is independent of temperature, the additional equation (equation of state) is 

p = p(p) (9) 


Summarizing the results so far, 

dp 1 

— + pV • v = q 

dt H - 

P^ + Vp = k-q(v-V')> 

P = P{P) 


( 10 ) 


2.3 Linear Fluid 

In the following section, some simplifications of the physical model of the fluid described 
by equations (10) are introduced, because they contain effects that can be neglected for the 
present purpose. 

The variation of the density is divided into two parts. One contains a time-dependent part, 
p, and the other contains the static part, p f . Thus 

p=p f + p 


where Vp, = 0. The first equation of (10) can be rewritten 


dp 

dt 


+ Pf 


v Pfj 


V y = q 


Assuming the variation in time to be small compared to the static value, that is 

M«i 

p f 


leads to 


-fr + Pp v = q 
dt 
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If the second equation of (10) is dealt with in the same way, the new system equations be- 
come 


dp - - 1 

-tL + p f Vv = q 

P f ^ + Vp = b-q(v-v v ) ■ 

P = p{p) 


(ID 


The total time derivatives of p and y are 


^e=^+v.vp 

dt dt ~ H 

dv dv , 

— =v+V'Vv 
dt dt u ’~ 


The convective terms in these expressions and the influence of the added fluid mass, q, on 
the second equation in (1 1) are insignificant under a certain condition. This condition oc- 
curs when the fluid velocity is considerably smaller than the speed of sound. Mathematical 
derivation of that simplification is not shown. 

Under such a condition, from equations (11) 

^- + P/V-v = < 7 (12) 

P f ^ + Vp = P f b (13) 


Together with the equation of state, equations (12) and (13) constitute a physical model for 
a linearized flow. 

A further simplification is possible. As was done in the case of density, the pressure is di- 
vided into two parts, p and pf (similar notation as for p), and the equation of state is lin- 
earized 


P = p f + P~ p f +p 


r dp] 

\ d PJ 


1 p=p i 


where Pf is assumed constant. Introducing a new constant 
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results in 


p = c f p 


Equation (12) is then differentiated with respect to time 


dq 


+ V. p,-i 


St St 


Substituting equation (13) into this expression yields 




In this relation, p is replaced by p from equation (14) and p is replaced by p+Pf implying 

^= C/ ! v i p+v ! P „-p / v*+^ 

From the definition of p, above, the term *P/ is zero and so the final result reduces to a 
wave equation with source terms, 

| £ = + f (15) 

By using equation (14), the wave equation can be expressed in terms of the pressure 

^ = c 2 f V 2 p-c 2 f p f V b + cj^- (16) 


Returning to equations (12) and (13), the gradient operator, V, is applied to the first and 
the partial time derivative to the second of these equations. Solving with respect to y leads 


3 =4V(V-,) + f -*V, 
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It may be noted that 

V(V • y) = V 2 v + V x (V x y) 

For a perfect fluid, as long as the body forces are derivable from a potential function, the 
velocity field, once irrotational, remains irrotational. Thus 

Vxy = 0 

This may also be regarded as a constraint; only conservative solutions are of interest. If 
this constraint is imposed on equation (17), once again, the wave equation is obtained in 
terms of the velocity field. 


d V 2r72 Cf r-j 


( 18 ) 


If the displacement field, u f , of the fluid is introduced, 



Equation (18) can be integrated with respect to time (disregarding the initial state for u f 
and b and the convective term as before) which results in 

^ = + (19) 

A P, 

where 

t 

Q = \qdr 

o 

Hence the displacement field also satisfies the wave equation. 


The assumption that the velocity field is irrotational implies that the displacement field has 
this property also, i.e., 


u f = Vy ( 20 ) 

where y/is a potential function to be introduced. Further, the potential associated with the 
body forces is expressed as 


b-V(p 


( 21 ) 
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Equation (19) yields 


v {^_ c jvv-* + Je}=° 

where the function inside the brackets is independent of position and depends on time only. 
Because an extra function of time can be added to y/ without changing the displacement 
field, u f , with no loss of generality 

= c 2 ^ 2 y/ + <p- ^-Q (22) 

dt P f 

is obtained. Alternative formulations derived in this subsection are 



(23) 


(24) 

<? 2 v 2 c ?2 db cl 

— f = clV 2 v + -^--- L Vq 
dt 1 dt p f 

(25) 

d l~2 f =c 2 f V 2 Uf + b-ivQ 

- P, 

(26) 

= c / vl ¥ + 0 - —Q 
dr p f 

(27) 


In the derivation of these equations, b was introduced as the source of influence inside the 
fluid. Later, the existence of a potential, <p, to the body forces was assumed. The relation, 
b = V<j) , can be introduced into equation (24), which would then lead to 

_ (28) 

where p f <t> can be interpreted as a prescribed pressure quantity applied inside the fluid do- 
main, 
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2.4 Boundary Conditions 

For a complete understanding of sound propagation processes it is obviously not sufficient 
to investigate the field equations for sound propagation. Additionally, the question "What 
happens at the boundary?” is of great interest. 

Assuming the fluid mass inflow, q, in the neighborhood of the boundary is identically 
zero, the fluid close to the boundary can now be treated as a continuum. A consequence of 
the continuum hypothesis is that the displacement field at any time must be a continuous 
mapping of the initial state. Therefore, boundary particles must remain as boundary parti- 
cles and interior particles remain as interior particles at all times. 

The fluid boundary is divided into five different parts according to their properties. The 
first four types of boundary conditions are especially important for the case of interior 
acoustics. 

B 1 : elastic boundary (fluid-structure-interaction/boundary type: Neumann) 

B 2 : open boundary (with prescribed external pressure/boundary type: Dirichlet) 

B 3 : rigid boundary (natural boundary condition/boundary type: Neumann) 

B4 : energy absorbing boundary (boundary type: Neumann) 

B 5 : artificial boundary to imitate radiation to infinity (boundary type: mixed) 

In the following paragraphs, it is desirable to formulate the boundary relations for the 
quantity 



where n is the outward normal vector to the fluid boundary. 

B1 (elastic boundary) - The motion of the structure and the normal component of the 
fluid motion coincide, that is 


UfQ = u sf (29) 

where u f is the fluid displacement vector and u sf is the structural displacement component 
perpendicular to the fluid boundary. The second time derivative of equation (29) yields 


dv d 2 ^ 
dt ~~ dt‘ 


(30) 
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Using equation (13) with p = p f + p (p f = constant) and b = 0 leads to 


dp__ n 

dn Pf dt 2 


( 31 ) 


B2 (open boundary) - The pressure at a point on the surface depends on the height of 
the wave created at the surface and on the external pressure, p e . 


P = P f g(Ufn)+ p t (32) 

where n is a unit vector normal to the surface in the initial state, which is independent of 
time. The quantity g(u f ■ n ) is essentially the weight of the wave created at the surface and 
g is the acceleration due to gravity. 


Differentiating equation (32) twice with respect to time, leads to 

dv _ 1 d 2 (p-p,) 
dt'~ p f g dt 2 


(33) 


dp = 1 d\p-p e ) 

dn g dt 2 


(34) 


B3 (rigid boundary) - On the surface with no waves the displacement perpendicular to 
the surface is assumed to be identically zero, so 

u f n = 0 (35) 


or after differentiating twice with respect to time 


or 




(36) 

(37) 


B4 (absorbing boundary) - This type of boundary is meaningful only in the frequency 
domain. The boundary condition belonging to this type of boundary is discussed in section 


4.1. 
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B5 (artificial boundary approximating radiation to infinity) - If it is required to investi- 
gate sound radiation from an interior to an exterior domain using discrete methods, the ex- 
terior domain must be truncated. This can be achieved in different ways, for example, by 
infinite element techniques or boundary integral methods. 

The main purpose of this report is to present the use of MSC/NASTRAN in fluid-structure- 
interaction problems, especially for the case of an interior fluid domain surrounded by elas- 
tic, rigid, or absorbing walls. The discussion of radiation to infinity is based on the Som- 
merfeld’s method. This method is exact if the incident wave is perpendicular to the bound- 
ary. If the boundary is located a great distance from the sources of disturbances, it is quite 
accurate. The method is also easy to understand. A derivation of the necessary expres- 
sions is given below. 

Assume that the boundary consists of a series of dampers. The force developed in a 
damper is 


d B 


du f 

~cT 


n 


where d B is the damping of the fluid on the boundary and 

du f 

~dT'~ 

is the velocity component perpendicular to the boundary. This term should equal the fluid 
pressure, thus 


P + Pf=d B -jj-n (38) 

Then 

^ = d B ^jfn = d,c}V(V- u,)-n (39) 

where the wave equation given by equation (26) with b = 0 and VQ = 0 is used. Through 
the use of equations (12) and (14), noting the condition that no fluid is added at the fluid 
structure interface, and assuming the initial state to be at rest, leads to 


p f c 2 f Vu f = -p 


(40) 



Combining equations (39) and (40) results in 


it = - dB -Lvp.n (41) 

dt p f 

where Vp • n = dp/dn is the normal derivative of the pressure at the boundary. Far away 
from the boundary the wave propagates at the speed of sound, c f . Therefore, 

— = -CfVp-n (42) 

at 

is just the d'Alembert's solution to the wave equation, so 


d B — p f c f 


which determines the value of the damping constant (acoustic wave resistance Z 0 ). It 
should be noted that this quantity is exact only for the case of plane waves. Equation (42) 
can then be written as 


rr ^P 


1 dp 

Cj dt 


(43) 


or according to equation (13) with b = 0, from equation (43) follows 

dv _ 1 dp 

dt ~ d B dt 


(44) 


Note that the feature of the boundary condition as expressed in equation (43) is to transmit 
a right-angle incident wave. 


2.5 Coupled Initial Boundary Value Problem 

Before continuing to the discrete, i.e., FE formulation of the coupled problem, the follow- 
ing formulation of the problem is presented as a summary of sections 2.1 to 2.4: 
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STRUCTURE (see section 2.1) 


FLUID (see sections 2.2 to 2.4) 


Dynamic equation of elastic media 
(displacement formulation) 


L{u s ) = F{rj) 


where L is a partial differential operator 
with respect to space and time, u s is the 
structural displacement and 


F(r,t) = F s (r,t) + F p (r,t) 


F s (r,t ) is an external time-dependent load 
and E p (r,t ) the fluid pressure-induced 
force onto the structure. 


Additionally, initial and boundary condi- 
tions must be given. 


Wave equation (see eq. 24) 


^ = c}V 2 p-c}p,Vb+c}^ 


where p is the pressure depending on space 
and time, p = p(r,t). 


Boundary Conditions: 


1 . Elastic boundary (B 1 ) 




d 2 u l 


dn Pf dt 
( u sf is the component of the structural 
displacement u s perpendicular to the 
fluid boundary). 


2. Open boundary (B2) 

P = o 

(special case: "acoustically soft"). 


3. Rigid boundary (B3) 


^ = 0 


dn 


(special case: "acoustically hard"). 


4. Absorbing boundary (B4) 
(see section 4.1) 


5. Radiation condition (B5) 


dp 1 dp 


dn 


c f dt 


Additionally initial conditions must be 
given. 


2.6 Coupled Boundary Value Problem 

In order to complete the mathematical formulation of the coupled FSI problem, the steady- 
state version of the coupled problem discussed in section 2.5 is presented. Assume a 
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STRUCTURE 

FLUID 

Dynamic equation of elastic media 

Helmholtz equation 

V 2 p + k 2 p = F f (r) 
where p = p(r), k = (o/c f 
( k : wave number, ft): radian frequency )■ 

Boundary Conditions 

1 . Elastic boundary (B 1 ) 

dp 2 

2. Open boundary (B2) 

P = o 

(special case: "acoustically soft"). 

3. Rigid boundary (B3) 

^ = 0 

dn 

(special case: "acoustically hard"). 

4. Absorbing boundary (B4) 

(see section 4.1) 

dp . 1 

7n=- mP >T, P 

(Z„: normal impedance of the ab- 
sorber). 

5. Radiation condition (B5) 

— = -ikp with k = Oi/c, 
dn 1 

(displacement formulation) 

L(u,) = F(r) 

where Lisa partial differential operator 
with respect to space, u s is the structural 
displacement and 

E(r,t) = F s (r,t)+F p (r,t) 

is an external time-independent load 
and F s (r) the fluid pressure-induced load at 
the structure, which is also time-indepen- 
dent. 


harmonic time dependence for the structural displacement u s (r,t), external load, F s (r,t), 
pressure-induced load F p (r,t), and pressure p(r,t). The sources inside the fluid will be 
summed in F f (r,t) and assumed to have a harmonic time dependence also. For simplic- 
ity, the notations for the time-independent quantities and the operator L are the same as in 
the time-dependent case. 
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3. FINITE ELEMENT MODELLING 


The content of this section relates to the discretization of the equations derived in the previ- 
ous section using FEM. In the fluid domain, the FE formulation is derived using the pres- 
sure as an independent variable. A reformulation is carried out using the weak form and 
Galerkin's method of choosing the test functions. In the first step, some notes referring to 
the discretization of the structure are presented. 


3.1 Structure 

As stated in section 2.1 the structural behavior is governed by 


L(u s ) = F(r,t) (45) 

with 

E(r,t) = F,(r,t) + F p (r,t) 


where F, is an external load and F p is the fluid pressure load on the structure. 

A discretized formulation for the structure yields, in matrix notation, 

[M,]{</} + [Z>,]{t/} +[*,]{£/} « {V,} + {L,} (46) 

where [M s ], [dJ, and [AT,] are the structural mass, damping, and stiffness matrix, respec- 
tively. {£/} is the structural displacement vector, j is the load vector due to external 
structural loads, and {l, J is the load vector due to the coupling effects. 

In the derivation of equation (46) the test functions and the trial functions are of the same 
type, although they may differ from the ones used for the fluid domain. That is why this 
function is designated by a subscript 's', i.e., N s and the function set by {Wj. Thus the 
coupling vector, L f , between the structural and the fluid domain is 

L f = jN,npdS (47) 

B i 


3.2 Linear Fluid 

Referring to equation (24), the following form of the wave equation is used as a starting 
point: 
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(48) 


^ = c}V 2 p + c}^-c 2 f p f Vb 


dt 


Multiplying equation (48) by a test function (or "variation") w=w(x,y,z) and integrating: 
over the total fluid domain volume, V, results in 

jw^dV-c 2 f jw^ 2 pdV =c 2 f lw^dV-cjp ; fwVbdV (49) 


Applying Green's first formula to the second integral, in equation (49) yields 
JwV 2 p dV = Jvv(Vp)-n dS - J(V»v)-(Vp) dV 


where S is the boundary of V and n is the outward normal of S. Inserting this expression 
in equation (49) leads to 


j w ~jPT + c] J (Vw)(Vp) dV = c}\ vv(Vp) • ndS + cj J w^j- dV - c 2 p f j wV-bdV (50) 

1/ &*' 1/ Cl/ 1/ 


dq 


Noting that p = p 0 + p (p 0 = constant) from equation (13) one obtains 

Y7 <?v 

Vp = -Pf-£+P f !> 

Reformulate the first and the third integral on the right side of equation 50. It becomes 
c 2 J.vv(Vp) • ndS - c 2 p f j wV bdV = -c 2 p f Jw~ ndS + cjp f j wb ndS- cjp f j wV • b dV 

S V s *** S V 

= -c 2 p f ^ W ~pt ~ + C fP/j^ W '- ^ 

C ^ I f 


The next step is to split the surface integral in equation (50) into the four boundary parts, 
(Bl, B2, B3, B5) discussed in section 2.4. To that end, note the reformulated integrals 
above: Splitting of the surface integrals leads to 




c ' p 'f * ds= c )pf J dS 


-18r 


(51) 



c}p,f w - ndS= ’jw ( 

B 2 Ul £ b 2 


(52) 

c 2 f p f \w-^--ndS = 0 

i * 


(53) 

4p,lw%.ldS-c,l$dS 

B S B i 


(54) 

Now, equation (50) is transformed to the following form: 



f W l£ ™ + C 'J ( Vw )' ( Vp ) dV = c / p / \ w J dS+ l\ W 

v & V B x 01 8 B 2 

-c f j dS + c’J w ~z~ dV - cjp f j (Vw) • b dV 

B t dt v dt v 

a 2 (p- p .)J 

dt 2 

► (55) 


If there are no body forces and no external applied pressure, p e , equation (55) reduces to 


jw&dV + c I ,j(Vw)(Vp)d v = -cjp r \ 

1 / 1 / P 


w- 




d\ 

dt 2 ™' R r dt 2 




<?<7 


A 


A 


(56) 


The discretization of equation (55) is carried out by expanding the pressure, p, in terms of 
FE basis functions or shape functions, each one associated with a unique nodal point. The 
shape functions in the fluid domain are designated by a subscript 'f, i.e., N f , and the 
function set by {A^}. 

The discretization process of the pressure, p, leads to the expression 

= <57) 

j 

The summation is over the number of shape functions and Pj(t) is the value of the pressure 
at the associated nodal point at time, t. 

Finally, the standard Galerkin formulation with 
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is applied to equation (56) and noting equation (57), gives the discretized form of equation 
(55) 


KH^l + NM + \ K ,V) = ■ -M + K } + M + M (5® 

where {P} is the column vector of the unknown nodal values of the pressure, and 

KW N^',dV + ijN‘ r N‘<IS 

V o B 2 

( O r )‘c r jN‘ f N>dS 

B, 

(*/), = cjj(W‘pNj)dV 

V 

(Ar),- = C f dS 

M, = c il N ‘fj; dV ’ <*»), = c >t \ ( w / ) • ft <fl ' 

v 01 v 

^■- c il N ‘^ ds 

are the elements of the matrices [M/]» [Df \ [^/]> an d the vectors {Lj, {L v }, {Z, fc }, a °d 
{L,} where i is the row and j is the column index. 

The vector {L,} is responsible for the coupling of the fluid with the structure and repre- 
sents the action of the structure on the fluid surface. The counterpart of {Zy} is {Zy } (see 
eq. 47), which describes the action of the fluid on the structure. 

3.3 Fluid-Structure Coupling and Symmetrization 

Following the developments in sections 3.1 and 3.2, a matrix equation of the coupled sys- 
tem using equations (46) and (58) can be written. As mentioned above, the coupling be- 
tween structure and fluid is fixed by the coupling terms {Zy} and {Zy}on the right-hand 
side of equations (46) and (58) respectively. 

| Zy} is a function of the fluid pressure and {Zy} is a function of the structural displace- 
ment. All other vectors on the right-hand side Of equations (46) and (58) are true load 
vectors. 
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The definition of {L,} and {ly} is 


and 


Because 


(L)=c}p,\N‘ t u,dS 

St 

(L f )=jN! s -npdS 

Bx 

u sf = u,-n = Ni-n Uj 


where the nodal index,;, runs through the set of structural trial functions, results in 


where 

W^cjP/j^/Ki-ndS 

b, 

Similarly, 

P=2>/P 

j 

and thus 
where 

(K c \ = \u!,-nN>dS 


Only those structural degrees of freedom that are perpendicular to the fluid boundary need 
to be considered. The same is valid for the fluid degrees of freedom, i.e., only the pres- 
sure nodes along the structure, that are linked to the structure, need to be considered. 


The assembled system of equations is 


X 

Op" 


D, 


'u 

1 

1 

I 

& 
1 

~ir 


i 

y 

i 

M c 

M /JlA 

+ 

0 

1 

cT 

P 

T 

L° 

*/J 

p 


L' + L b + L t \ 


(59) 


It is obvious that M c - cjp f K T c . 

If the second row of equation (59) is multiplied by the factor 

Vx 
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another form of equation (59) can be obtained, where % = c 2 f p f is the bulk modulus. For 
the sake of brevity, the same notation for M f , D f , K f , L q , L b , L e as in equation (59) is used 
which leads to 


' M, O' 

'0' 


'D, 0 Yu 


K, C" 

'U' 


0 

-C T M f 

p 

+ 

0 D,]p_ 

+ 

0 K f 

P 


_ L q + h + 


(60) 


where C (Coupling Matrix) is the matrix -K c from equation (59). 

Consider the following interesting aspects of equation (60). Assume that the acoustic fre- 
quencies are all greater than the maximum structural frequency of interest. In this case, the 
fluid can be considered to be incompressible. This corresponds to letting the bulk modu- 
lus, X » be infinite. This leads to M f and D f => 0. 

It is easy to see that for this case the coupled system can now be completely described in 
terms of structural nodal displacements only. That is, equation (60) reduces to its structural 
part, where an additional mass term reflects the total effect of the incompressible fluid on 
the response of the structure. 

It is interesting to observe that the assumptions of incompressible, small deformation be- 
havior (the latter assumption was made for the acoustic approximation) of the fluid elimi- 
nates the need for a coupled solution. The presence of the fluid is entirely reflected by the 
additional mass term for the structure. 


3.4 Aspects of Solving the Coupled Fluid-Structure Equation 

3.4.1 Symmetrization 

For the solution of the general FSI problem, matrix equation (60) has to be solved But an 
important feature of the above equation is its lack of symmetry, i.e., the system has un- 
symmetric mass and stiffness matrices. This can present some difficulties in extracting the 
eigenvalues and in the calculation of the dynamic response of the system. 

Eigenvalue equation of the coupled system can be obtained if the right-hand side of equa- 
tion (60) is replaced by a zero vector, D t — D f = 0, and harmonic time-dependence of U(t) 
and P(t) is assumed, i.e.. 


U(t) = Uexp(icot) 

and 


P(t) = Pexp(icot) 
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respectively. The vectors U and P are now dependent only on the nodal position. The 
eigenvalue problem of the coupled system leads to the following nonsymmetric equation: 


a*. ci 

/ M ' 0 T 

W' 


roi 

1 

1 o 

" |_-c T mJ 

k 


i 

o 

i — 


(61) 


Equation (61) is costly to solve and it is difficult to obtain eigenvalues for large systems. 
With regard to practical problems, the existing capabilities of MSC/NASTRAN are not suit- 
able for solving systems of the type in equation (61). Furthermore, it can be shown that 
the eigenvalues of the coupled system in equation (61) are real. The nature of the problem 
suggests a symmetric formulation for equation (61). A symmetric form can be obtained 
directly by using the displacement formulation for the fluid (ref. 9). This approach has a 
major drawback because of an increase in the number of unknowns, especially in three-di- 
mensional applications. In the literature one can find many descriptions of symmetrization 
procedures for equations (61) and (60), see references 21-28. The method used in the 
Acoustic Procedure presented here is described in references 25 and 28. The characteristic 
feature of this method is that the symmetrization is accomplished by replacing the vector P 
in equation (60) by a vector Q given by 


P = Q 


This vector is, except for a multiplicative constant, the velocity potential used in fluid dy- 
namics. 


A symmetric form of equation (60) is obtained, if the second row of equation (60) is mul- 
tiplied by -1, integrated with respect to time, t, and finally {P} is replaced by {()}• The re- 
sult is the following system: 



0 

-M, 


D. 





(62) 


where 


G{t) = \{L q + L b + L)dx 
0 


The disadvantage of this symmetrization is, that in case of the eigenvalue calculation, it 
leads to an additional symmetric but fictitious damping matrix, which makes it necessary to 
use complex eigenvalue methods. This situation can be avoided by applying the following 
matrix transform: 
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0 


(63) 


r = 


km: x 

c t m;' 



where I f is a unit matrix with the same size as K f . Multiplying equation (61) from the left 
by equation (63) results in 


T kjuj'k, -kju;'c 

-/» 2 

X o -' 

T V 


' O ' 

[-c t m;X K f +c T M; l c 

— UJ 

0 Mf 

p 

/L J 


0 


(64) 


In practice, the systems represented by equations (62) and (64) are very large. There exists 
a method that can reduce the expense in solving these systems. It involves the use of 
modal reduction. 

3.4.2 Modal reduction 

For very large problems it is computationally attractive to apply modal reduction for the 
purpose of reducing the size of the analysis. In this case, both the fluid and the structure 
are possible candidates for this operation. 

Applying modal reduction to the structure (for the modes of the structure in vacuo) and to 
the fluid (for the modes of the fluid in a rigid container) results in the equations listed in the 
following table: 



Structure (s) 

Fluid if) 

Variables 

U = <P£. 


Modal mass 

m, = <pjM<t> s 

II 

-O- 

V 

Modal damping 

** 

& 

II 

df — <P f D(p f 

Modal stiffness 


•o- 

ii 

Modal excitation 

fs=^A 

ii 

■o- 

Q 


The matrices </> s and (p f represent the modal matrices of structure and fluid respectively. In 
general, both are rectangular matrices, where the number of rows is equal to the number of 
physical degrees of freedom and the number of columns is equal to the number of modes 
used in the analysis. 

Following the transformation process given in the table above, equations (62) and (64) are 
converted to a set of "modally coupled" equations, where the coupling is represented by the 
"modal coupling matrix" 
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(65) 


C,=<p T f C<p s 


The resulting equation of motion for the modal coupled system is 


K °1 

'{/ 

1 

1* n 

'L' 

_L. 

1 

0 

1 

j; 


1 

< 

1 

1 

0 

1 

4* 

1 

Jf. 

T 

1 

1 

Jf. 

T 

1 

1 

o 
1 

Jf. 


L~*/J 


The eigenvalue equation for the coupled modal system is 


* s , 

1 

vn 
1 

,7 

r 

I 

o 

u; 


"O' 

1 

+ 

1 

i 

- CO 

A 

1 

s 

o 

1 

Ur. 


0 


(67) 


Assuming a harmonic time-dependence for structural displacement and pressure, the fol- 
lowing frequency response equation is obtained from equation (66): 


T*, 

o " 

- CO 2 

X 

0 ' 

+ i(Q 

X 

C ■' 

[X‘ 


■ / 's' 

L° 

~ k f. 

0 

~ m r. 

J 


Ur. 


-*/_ 


( 68 ) 


3.4.3 Residual flexibility for the fluid 


If the modal reduction is used, the analyst is alerted to the possible effects of mode trunca- 
tion inherent in the method. Enhancements can be made to the above equations of motion 
by adding terms for the fluid model residuals. This addition compensates for the omission 
of high frequency fluid modes. Normally, high frequency modes respond "statically" and 
residual flexibility accounts for their static behavior (ref. 30). 


The derivation of the residual flexibility terms is done in terms of an unsymmetric formula- 
tion (see eq. (60)). For simplicity the damping matrices in equation (60) are ignored be- 
cause they have no impact on the final result. Likewise, for simplicity, assume that there 
are no body forces, no added fluid mass and ho external time-dependent pressure, i.e., 
b = 0, <7 = 0, p e = 0 . Therefore, from equation (60) one obtains 


'Ms O' 

-C T M, 


U 


K s +C 
0 K r 



( 69 ) 


After transforming the lower part of equation (69) to modal coordinates, the "fluid part" 
yields the following equation: 
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m /(£/ + ["/]£/) ^/^ 


(70) 


where j© 2 ] = nrfkf is a diagonal matrix and ft) 2 are the fluid eigenvalues. The matrices 
m f and k f are modal mass and stiffness matrices of the fluid. The next step is to separate 
the fluid modes in the following manner: 


0 / - [ 0 /| 0 /| 0 /] 

where 

<l>f : "zero" modes with (0° f =0 
<f>‘ f : "elastic" modes with (o e f > 0 

{ p' f : "truncated" modes with 


(71) 


In the same manner the quantities, , m f , and [cu 2 ] can be separated, i.e., 


€/ = 


ii 

sr 


• Kl= 

I 

* V, 

3 

0 ^ 

3- 

1 




i 

s 

^ - 
1 


1 

- ^ 

3 

1 


(72) 


Then, from equation (70) the following is obtained: 


^ f = (m°y<p 0 f T C T U (73) 

= (m € f )~y + (o e ;Y<p; T c T u ( 74 ) 

= (75) 


where the following have been used: case (a) 0 = s 2 U and = s 2 ^° f , case (b) t,‘ s = s 2 £,‘ f , 
and case (c) the condition t0/ 2 ]£/»|/ . 


Using equations (71), (72), (73) through (75), and P - <p f ^ f , the following expression is 
obtained: 


P = <t>° f (m°y 1 <l> 0 f T C T U + <p e f (my\s 2 + (ofY <pfC T 0 + Z‘ f C T U (76) 
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where 


(77) 


z;=f >KK]f ♦: 


represents the residual flexibility. 

Now, equation (76) can be solved for Z‘ f by separating the static from the dynamic part of 
the pressure, i.e., 


P = P‘ + P d 


Set s - 0 (i.e., P = 0) and get 

P a = <p° f {m° f y'(p 0 f T C T U + <p° f (m' f y l [s 2 + 0}‘ f 2 \ X <pfC T V + Z)C T U (78) 

Multiply equation (78) by K f and note that the exact static solution is known from equation 
(69): 


K f P s -C T U 


(79) 


Furthermore, for zero modes 


K f <P°f = 0 


Then from equation (78), the following expression is obtained: 

K f Z' f = -K f <p f (m f y\(oy]<p T f (80) 

By constraining zero modes, equation (80) can be solved directly for Z ‘ f , i.e., 

z‘ f =-<t> f (m f y\ (oyy f (sd 

Further define 

zy = <P°f(m° f y l <f>y m 

Using equation (82), the following is obtained from equation (76) 

/ > = ^ / + z; , C r {/ + z;c r U (83) 

Combine equations (70) and (83) with the upper half of equation (69) and transform to 
structural modes to obtain 
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( 84 ) 


m f + Z'f 

o' 

'f/ 

+ 

’k s + z7 

c' 

As' 


7/ 

t 

1 


_i t _ 


0 

k f_ 

Af. 


0 


where 

Z'f = <t>JCZ‘ f C T <p s and Z°/ = CZ°/C T <P S 


This section is concluded by noting that if there is an excitation in the fluid the residual 
flexibility terms to be added are given by the following two equations 


Z' f = <p T s C(K - 0/rc; 1 (0}Vf)<Ps (85) 

Z/ =0jc(0 / m / - 1 0;)C> J (86) 


After symmetrization, this addition leads to an equation which has the following form: 


m s + Z' f 

0 ’ 

'iV 

+ 

i 

■*- ; 
i 

'L' 

+ 

i 

+ 

N1 
^ s? 

0 ' 

As' 


'/. s' 

0 


if. 


-*/J 

if. 


L 0 

~ k f. 

Jf. 


_~8f_ 


- 28 - 



4. SPECIAL FEATURES 


4.1 Acoustic Absorption 

The effect of an absorbing wall in an acoustic cavity can be modelled in two ways, depend- 
ing on the nature of the absorption. If the absorption is due to a porous material then spe- 
cial porous elements must be used. The porous elements have to model the behavior of a 
flexible porous material under certain assumptions regarding the structure of the material. 
At the present state of the program, such elements are not available. The other possibility is 
to model the absorbing wall, both in the case of a porous material or in a general type of ab- 
sorption, by using the impedance formulation of the absorber. In order to obtain a bound- 
ary condition in this case, the boundary condition for boundary B1 (see section 2.4, eq. 
31) can be rewritten by the use of a relation between pressure, p, and the normal fluid ve- 
locity, v„ = du 3f /dt, at the wall 


to give 



dp . 1 

— = ~i(Op f —p 


dn 


( 88 ) 

(89) 


To use the impedance for modelling absorptive boundaries, the value of Z„ has to be de- 
termined as a function of the radian frequency, (0. In some cases an analytical expression 
can be formulated, whereas in others, the values can be obtained from experimental mea- 
surements. In the latter case, the so-called "standing wave tube" is used, which can be 
used to measure the (normal) impedance Z n (see ASTM El 050-85 A) 

Z n = Z, + iZ 2 (90) 


where Z, is usually called the specific acoustic resistance and Z 2 the specific 
acoustic reactance. 

If the absorber is connected to an elastic boundary, the boundary condition can be ex- 
pressed as 




dn 


icop f \ 


f P ' 
v„ + — 

V 


(91) 


A surface with a "large" impedance, |Z n | => «>, is considered to be "acoustically hard." A 
surface with a "small" impedance, |Z„| => 0, is considered to be "acoustically soft." Relat- 
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ing to equation (91), the first case leads to the elastic condition in section 2.6, where 
v, = dUjf dt , while the second case leads to p = 0, i.e., "pressure release." 

If the enclosed acoustic domain is discretized with finite elements, the equations of motion 
are obtained as described in section 3. However, the damping matrix, D (see equations 
(60) and (62)), has to be complemented by the boundary absorption, D t . It is easy to see, 
that D z is given by 


A = 



(92) 


where A t is the matrix representing those parts of the boundary which are acting as ab- 
sorbers. 

For a lumped approximation, D z is null except for terms corresponding to the grid point lo- 
cation of the absorber material. 

The Area Matrix Program generates the matrix A t . With an additional input of the material 
property of the absorber, it also generates the damping matrix, D z . Usually, the impedance 
property of the absorber is described by the specific acoustic admittance defined by 

/3=|> (93) 

A 


where Z 0 = p f c f is the wave resistance for plane waves in the fluid. It should be noted that 
Z 0 is a real quantity only in the case of plane waves, because there is no phase shift be- 
tween pressure,/?, and the velocity, y, of the fluid particles, if the fluid is inviscid and ho- 
mogeneous. 

Using equation (90), the following expression is derived from equation (93) 


where 


P = y + io 



and 



(94) 

(95) 


where Z„ is defined by equation (90). Therefore, instead of equation (92), the following 
expression is obtained for the damping matrix, D z : 

D, = Ur+io )A, (96) 

A 
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That is, as input of the absorber’s material property, the quantities yand a are used, which 
are defined by equation (95). They determine the specific acoustic admittance defined by 
equation (94). 

The condition for the validity of the normal impedance is that only plane, normal incident 
sound waves will be correctly damped. If this condition is violated then reflections will oc- 
cur. 

It is possible that reflections may, in some cases, corrupt the solution and the results are not 
what the analyst expects. To avoid this situation, it is possible for some problems to esti- 
mate when sound waves, originating from a source, are plane. Nevertheless, it is ques- 
tionable whether the impedance model really does represent the absorbing surface correctly. 
There are examples where this model is completely confusing. 

But these reservations, which can be theoretically proven, are diminished by the experience 
of others especially with regard to applications in automobile interior acoustics (ref. 14). 

An improvement for the simple impedance model implemented in the present Acoustic Pro- 
cedure is possible, if instead of equation (88), the following condition proposed by D.B. 
Bliss (ref. 16) is used 


p + B(co)V 2 B p = v£ n 


(97) 


where V fi is the gradient operator on the absorbing boundary (B4) and (B5). 


B(o)) = 


1 

k(cof 


is the so-called bulk-reacting coefficient and k(co ) is the complex wave number. The 
condition in equation (97) can be applied to describe the boundary absorption effects. This 
procedure was implemented in MSC/NASTRAN by Burfeindt (ref. 15). 


4.2 Acoustic Contribution Analysis 

This section presents a brief description of the calculation procedure for determining the 
contribution of different parts of the structure to the sound pressure level at any point in the 
fluid domain. A more detailed presentation is given in reference 31. 

The modal solution method provides valuable information regarding the structural and 
modal participation in the acoustic response. The acoustic modal participation at the 
interior points, {/^} , can be calculated by 
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( 98 ) 


{/*} - ZfiQfi) 

i=i 

where n f is the total number of acoustic modes used in the analysis. Each vector {f*} rep- 
resents the contribution to the acoustic response from the acoustic mode, or the modal 
acoustic participation factor. 

The structural modal participation factor can be calculated for a harmonic solution 
of radian frequency, ct), as follows. First, the equation of motion is solved for the fluid 
modal vector: 




(99) 


where [£] is the modal coupling matrix defined by equation (65) and 

[z2] = (-C0 2 [m f ]+[k f ]) ' 


( 100 ) 


is the acoustic modal frequency response function, which is included also in equation (67). 
Substituting equation (99) into P = <p f ^ f results in 

= = < 101 > 

7=1 

where n s is the total number of structural modes. Note that each vector P = {/*} represents 
the contribution of the structural mode j to the acoustic response or the structural modal 
participation factor. 

The boundary panel participation for a boundary panel b is basically the collective 
structural modal participation of all the structural nodes on that panel. The [C 6 ] matrix, 
which is the coupling matrix due to panel b, has to be computed and the global coupling 
matrix [C] will be the assembly of all panels in the structure. Once the [C fc ] are known, 
equation (65) can be written as 


KH*/f£[ c *M < 102 > 

b = 1 


which when substituted in equation (102) becomes 
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( 103 ) 


p = = iw 

6= 1 6=1 

where n b is the total number of boundary panels. Note that each {P b } vector represents the 
contribution of boundary panel b to the acoustic response. 


4.3 Structure-Fluid Analogy in MSC/NASTRAN 


Since there are at present no "fluid elements" in MSC/NASTRAN, the analogy between the 
structure and fluid equations has to be used. The following theoretical development is re- 
produced from reference 23. Rewrite equation (16) neglecting the source terms and ex- 
pressing Cf in terms of the bulk modulus % and density p f of the fluid. The following form 
of the wave equation (in Cartesian coordinates) is obtained: 




dy^P/fy) 


d_ 

' \ _c7 

1 d*p 

dz 

,Pi^, 

x* 2 


(104) 


An analogous equation in structural mechanics is the equation for the equilibrium of 
stresses in a particular fixed direction: 


dc ** , dt *y , dt *- 0 
dx dx dx * dt 2 


(105) 


where u SI is the structural displacement in the x direction, and T„ are stress compo- 

nents and p, is the structural mass density. 

In order to establish the acoustic-structural analogy, let 

= p (106) 


Ps=- 


1 dp 

< 7 ^ = — 

pfdx 

1 dp 


= ~U, 




P/dy 

_ 1 dp 

T ‘~7,7 


= -u 


fy 


= ~U f 


(107) 

(108) 
(109) 
(HO) 
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where the second forms of equations (108), (109), and (1 10) are, obtained from equation . 
( 10 ) with s b '= 0 and q = 0 and where, u A is the ^-component of fluid acceleration, etc. 


In order to complete, the analogy, note that, if structural displacement components u^, and* 
u n are set 1 equal to zero,, the general stress-strain relationship,, provided; by. the,: 
MSGVNASTRAN MAT9 material card is 


° xx 


~pn G lA G 16 ‘ 




► = 

G 44 g 45 


7*' 


! 

jym g 66 _ 


.7xx. 


where Gy is an element of the 6 x 6 elastic material matrix, and the strain-displacement rela- 
tionships are 


Therefore, if 


Yxy 

Jxi 


jjfg. 

dx 

dUsx 

dy 

dz . 



( 112 ) 


(113) 


equations (111), (112), and (106) will produce equations (108), (109), and (110) exactly. 
The . other, components of the [Gy] may be set , to any values, including zero, since e„, 
and Yj* are all zero. Equations (107) and (113) indicate the manner in which the MAT9 
card should be filled out. Note that several MAT9 cards may be used to represent the fluids 
where P/ and x vary with position in the fluid. The basic coordinate system should be 
used to define components of displacement and also the material coordinate system. (The 
symbol ,0 (zero) must be inserted for CORDM on the PSOLID card.) 

Note that, from equations (108), (109), and (110), normal NASTRAN stress data recovery 
will produce the accelerations components within the fluid, with a change in sign. 
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4.4 Boundary Conditions and Their Input in MSC/NASTRAN Using the 

Structure-Fluid Analogy 

In section 2.4 boundary conditions are discussed in a more mathematical sense. The pur- 
pose of this section is to present a method of implementing the boundary conditions for the 
fluid with regard to FSI calculations using MSC/NASTRAN. 


1 . At free surfaces (open boundaries), set u x = p = 0 by means of an SPC card. 


2. At rigid walls, take no action. The acoustic boundary condition, 


dp 

dn 


d P + n d P,„ fy-n 

n -Tx +n ’^ +n -Tr° 


(114) 


where n x ,n y , and n t are the direction cosines of the normal, implies an analogous 
structural boundary condition 




+ «,t x> + « 1 t jb = 0 


(115) 


where f x is the x-component of traction at the boundary. Satisfaction of equation 
(115) requires only that no structural forces be applied to grid points on the bound- 
ary. 

3. At a surface where the pressure is a known function of time, set 

= P( 0 


using the large mass method described in the MSC/NASTRAN Application Man- 
ual. 

4. At a surface where the normal component of the displacement, u fn , velocity, ii fn , or 
acceleration, u fn , is known, apply a grid point load 

= A Mt) (116) 


where A, is the area associated with the grid point i. Equation 1 16 follows from 
equations (108), (109), (110), and (115). Note that u fn , u fn , and ii fn must be in- 
ward to the fluid. 

5. At a surface where the pressure and the normal velocity are linearly related, i.e., 
where 
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(117) 


p = -Wu fH 


connect a CDAMP2 element to u sx with a damping coefficient 



(118) 


where, again, ^ is the area associated with the grid point i. 

6. Boundary conditions for other cases can easily be worked out. For example, if 

P = -Hu f „ (H9) 

connect a scalar mass, M = A Jfi to the boundary grid point i. 
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5. DESCRIPTION OF THE ACOUSTIC PROCEDURE 


5.1 Introductory Remarks 

This section presents the summary of the capabilities of the Acoustic Procedure to perform 
FSI analysis using MSC/NASTRAN. It is specifically aimed at the acoustic and noise 
control analysis for the interior acoustic domains, e.g., passenger cabins in automobiles 
and aircraft. 

The goal of this development is to treat the fluid medium simply as another part of the 
structure. This allows for a simultaneous solution of the two systems, without the need for 
a separate analysis, as is presently done in some cases. The only difference between the 
standard MSC/NASTRAN dynamic analysis and the Acoustic Procedure is the addition of 
some Direct Matrix Abstraction Programs (DMAP) and a single FORTRAN program. 

The procedure requires the use of superelement solution sequences, but it is not necessary 
to define superelements. In many cases, though, it is highly desirable to use superele- 
ments. Since MSC/NASTRAN currently does not have special fluid elements, Structure- 
Fluid Analogy (see section 4.3) is used. This capability is particularly applicable to the in- 
terior problems; that is, for problems in which the fluid is contained by the structure. The 
exterior problem, fluid outside the structure (e.g., submarine) is very difficult to solve with 
FEM. In this case boundary element methods are more suitable. 


5.2 Organization and Flowchart of the Acoustic Procedure 
The analysis capabilities that have been developed include the following: 

1 . Uncoupled normal modes analysis of the structure and the fluid, SOL 63 

2. Normal modes analysis of the "modal coupled" fluid-structure system, SOL 70 (see 
section 3.4) 

3 . Modal frequency response analysis of the coupled system, SOL 7 1 

4. Modal transient response analysis of the coupled system, SOL 72 

5. Absorption damping for the fluid medium (attention: only normal impedance, see 
section 4.1) 

6. Calculation of the structure grid point participation factors in the generation of the 
pressure at any specified point in the fluid (see section 4.2) 

The Acoustic Procedure consists of five NASTRAN runs. Figure 1 illustrates the sequence 
of these runs during the process of an acoustic calculation. 
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r RUN 4 ~ 

ACOUSTIC CONTRIBUTION 
ANALYSIS 

l DMAP J 


Figure 1. Acoustic procedure flowchart. 


The most costly run is the first run (see fig. 1). This run contains the FE models of struc- 
ture and fluid, the definition of the coupling and absorbing areas, and the absorber proper- 
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ties. The coupling matrix and the absorber matrices are formed by the Fortran program 
AMP (Area Matrix Program), which is also involved in Run 1, see figure 1. Run 2 is the 
cheapest run and the user has to decide if it should be executed. To better understand the 
coupling effects, Run 2 should be carried out. The response Runs 3.1 and 3.2 can be per- 
formed immediately, after Run 1. Run 4 (Acoustic Contribution Analysis) is carried out as 
a restart of Run 3.1. 

To get an idea of the CPU time, three examples of different size, i.e., with different num- 
ber of degrees of freedom (DOF), are presented. 


Example 1: This problem is a simple piston/tube model described in section 7.1. It con- 
sists of only 25 structural nodes and 16 QUAD4 elements. It has 1275 fluid 
nodes and 800 HEXA8 elements. 


Example 2: The following problem is described in Section 7.2 and it is a two-dimensional 
model composed of 17 nodes and 16 BAR elements for the structure and 313 
nodes and 128 QUAD8 elements for the fluid. 

Example 3: This example involves a solution of a complex automobile structure composed 
of 6 superelements. There are 46000 grids and 51000 assorted elements in 
the structure model and 800 grids and 700 solid elements in the fluid. 


Table 1. CPU time (sec) for acoustic analysis (CRAY YMP) 

Run 



Example 3 

1 

33.3 

10.9 

15066.0 

2 

2.9 

2.3 

633.6 

3.1 

8.4 

6.1 

1954.4 

3.2 

5.8 

3.4 

- 

4 

10.5 

1.0 

47.6 


In figure 1, the runs which have to be performed during an acoustic analysis were pre- 
sented. Figure 2 shows the essential steps in acoustic analysis. These steps include the 
modelling of the structure and fluid, defining the areas of coupling and absorption, fre- 
quency and time response calculations, and so on. 
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Figure 2. Possible analysis steps in the acoustic procedure. 
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6 . DESCRIPTION OF THE USER INTERFACE 


6.1 Modelling Aspects 


6.1.1 Fluid 

In the building of a finite element model certain aspects must be considered. The acoustic 
problem is a dynamic problem, hence the variation of the pressure field, both in time and 
space, has to be incorporated into the model. For the sake of simplicity, the discussion is 
restricted to the spatial variation of the pressure. From the basic Helmholtz equation, it can 
be shown that for a one-dimensional problem, 

p(x) = a, cos kx + a 2 sin kx 

where k = c o/c f . To obtain the behavior indicated in this equation, the finite element mesh 
has to be able to reproduce a sinusoidal variation. Since this variation is dependent upon 
the frequency, because of the wave number, k, it is clear that the size of the elements must 
be compatible with the wavelength of the highest frequency of interest. 

Given an upper frequency limit for a specific problem, the size of the elements can be de- 
duced. If ND nodes are required to discretize a wavelength, then the distance between any 
two nodes should be 


D = 


2it 

ND-k 


( 120 ) 


and the maximum size of an element side is of course 2 • D. Depending on the type of 
analysis and the desired accuracy, the value of ND can be selected as 

6 < ND < 12 (121) 

If the mesh is too coarse, the stiffness of the fluid will be exaggerated. It should be noted 
that the rule of thumb given above is applicable to any dynamic analysis problem. 

6. 1 .2 Coupling of structure and fluid 

For the case in which the fluid interacts with a flexible structure, the considerations in sec- 
tion 6.1.1 must be applied to the finite element model of the structure as well as to the 
modelling of the interface coupling between the fluid and structure. Special attention has to 
be given to the interface where, in general, two systems with large differences in wave 
length are coupled to each other, see figure 3. 
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Structure 



Fluid 


Figure 3. Comparison of wave length. 

The most obvious approach in this situation is to use the same fine mesh in both the fluid 
and the structure. But for most situations, this would result in many unnecessary nodes in 
the fluid. Thus, it is desirable to use a different mesh size in the fluid and in the structure 
even at the interface. This approach gives rise to a discretization enror because the effect of 
the structure upon the fluid is smoothed (assuming the situation in fig. 3). However, since 
this effect is not very important for the strength of interaction, it is acceptable. An interface 
similar to that shown in figure 4 would thus be possible. 

Structural nodes 




Fluid nodes 


Figure 4. Discretization at the interface between fluid and structure. 

It should be noted, though, that the above approach is by no means generally valid and the 
equations (120) and (121) must be considered for the fluid and the structure (ref. 32). 

6.1.3 Use of superelements 

One of the unique capabilities of this procedure is the capability to use superelements. In 
many large problems, especially for dynamic analysis, the use of superelements can signif- 
icantly reduce the cost of the analysis. There is no limitation on the use of superelements 
for fluid-structure coupling as long as the fluid model is in the residual structure. Any su- 
perelement can be in contact with the fluid, regardless of the number of superelement levels 
used. Currently, secondary superelements are not supported . Any superelement that is in 
contact with t he fluid must be a primary superelemcnt* 
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6.2 Fluid Media 


6.2. 1 Grid point definition 

The fluid media geometry is described by the use of the standard GRID and element con- 
nectivity cards. The fluid grid points should be distinguished from structural points by 
numbering the fluid points in a different range. Thus, partition between the two systems is 
easier to specify. The fluid grid points have "stiffness" and "mass" properties at only one 
of the six DOF at a grid. The user can specify the active DOF by constraining all other 
DOFs except that which is to be used. Typically, the T1 DOF is used (T1 is the 
NASTRAN notation for the DOF in ^-direction). Where the displacement coordinate sys- 
tem for the fluid grid point references a cylindrical or spherical coordinate system, it is re- 
quired that the T3 DOF (DOF in z-direction) must be the active DOF. 

6.2.2 Element definition 

The fluid medium can be modeled with any MSC/NASTRAN elements (in the three-di- 
mensional case, e.g., HEXA, PENTA, and TETRA elements). It is recommended that the 
user specify the consistent mass for the fluid. This leads to better results in practice. 


Caution: On the PSOLID card (property definition for solids) Integration 
Option 2 should be used for VERSION V65. 


6.2.3 Acoustic absorption 

The absorption capability has been provided and is defined by the DMIG input. DMIG 
matrices are generated by the AMP. Provision has been made to define frequency depen- 
dent absorption. For transient response analysis boundary absorption cannot be used since 
the derivation shown in section 4. 1 is valid only in the frequency domain. 

6.2.4 Loading 

The loading of the fluid is generally analogous to that used in structural mechanics analysis. 
Several types of loading can be specified for the fluid model: 

1 . Constant enforced pressure at the grid points 

2. Frequency or time dependent enforced pressure at the grid points 

3 . Acoustic source characterized by a volumetric flow rate 

The pressure can be enforced in the same manner as structural base excitations, i.e., "large 
mass" approach or Lagrange multiplier method. The volumetric flow rate is specified via 
standard structural load specification. Owing to the svmmetrization of equations of motion 
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in section 3.4. the forces applied directly to the fluid first have to be integrated in time and 
then multiplied bv -1.0. The resulting load is applied in the usual way. 


6.3 Uncoupled Normal Modes Analysis of Structure and Fluid (Run 1) 

This run is based on solution sequence, SOL 63, of MSC/NASTRAN. In general, the su- 
perelement capability is available to the user of the acoustic package. The input listing of 
Run 1 is shown in Appendix Bl. 


Caution: If superelements are used, the following limitations (mentioned in 
section 6.1.3) must be observed: 


1 . The whole fluid model must be in the residual structure. 

2. Any superelement that has any interior grid points that are in contact with the fluid 
must be in primary superelement. It cannot be an identical or a mirror image, i.e., 
secondary superelement. 

3. Fluid grid points must be "reasonably" coincident with an adjacent structure grid 
point. Coincidence is defined, see section 6.3.4. 1, by the AMP parameter, TOL 1. 
(Remember: An unconnected fluid surface grid point is assumed to have a 
dp/dn = 0 boundary condition, i.e., rigid wall). 

6.3.1 Executive Control Deck 

The DMAP alter must be inserted following the SOL 63 call. This instruction is in general 
automatically fulfilled by the actual machine-dependent procedure, i.e., no user action is 
necessary. If the analysis is performed in multiple steps, the DMAP alter must be present 
during the residual structure processing. It is recommended that the alter program be in- 
cluded for processing of all upstream superelements. There are no other special instruc- 
tions. 

6.3.2 Case Control Deck 

During the residual structure processing the user is required to do the following: 

1. Specify a loadset = xx to select an LSEQ card in the Bulk Data Deck, lseq 
will be described in the discussion of the Bulk Data Deck. 

2. Define set yy = i thru J such that all fluid grid points are in the range of i to J. 
This set can include any number of nonexistent grid points. 
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3. Specify disp (...) = all for all superelements that have any interior grid points 
in contact with the fluid. 


4. Define two subcases for Fluid-Structure Interaction. Note that there is ordinarily 
only one subcase for the residual structure processing. The first subcase contains 
the partn = yy, where yy references the set yy defined above and a method 
card selecting an eigenvalue extraction method for the structure only. The second 
subcase contains only the method card for selecting the eigenvalue extraction 
method for the fluid. 


6.3.3 Bulk Data Deck 

The structural model is defined in the usual way. The fluid model is defined in the same 
way with one exception. The material property cards define the bulk modulus and the 
density of the fluid. For two-dimensional elements use the MAT2 and for three-dimen- 
sional elements use the MAT9 cards. Due to the analogy described in section 4.3, the 
density of the fluid is 


and the G„ entries are 


VZ 

ypf 


Use of the parameter param, coupmass requires some consideration. This parameter 
selects the consistent mass formulation. If the user does not wish to use the consistent 
mass for all the elements in the model, the parameter can be specified in the Case Control 
Deck. Specifically, if the superelements are used then the coupmass can be specified for 
any superelement. In any case the param, coupmass should be used for the residual 
since the fluid is in the residual structure. 


Caution: The parameter specification must be placed in the first subcase 
and not in the second subcase of the residual. 


The two unique remaining items for the fluid model are the specification of the pressure 
load data to help define the 1) coupling areas (=> coupling matrix) and 2) absorbing areas 
(=> absorption matrix) of the fluid surface. The definition of these areas takes place in the 
two-dimensional case by force cards and in the three-dimensional case by PLOAD2 or 
PLOAD4 cards. 
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Caution: The pressure magnitude must be 1.0 and its direction must be in- 
ward to the fluid. If one wishes to use PLOAD2 cards the areas have to be 
described by uniformly oriented surface elements, where the orientation is 
prescribed by the direction of the surface normal using, e.g., the right-hand- 
rule. If an error is made in the definition of the direction, it can be corrected 
later in the AMP Input (see section 6.3.4. 1). 


The first pressure (force) load defined must be applied to the surface of the fluid in contact 
with the structure. 


The absorption areas are defined by specifying pressure loads on the surfaces at which ab- 
sorption material is in contact with the fluid. Note that the absorption area may be internal 
to the fluid, such as a seat in the automobile acoustic cavity. Each pressure load with a dif- 
ferent Set Id (Id: Identification) represents a separate absorption area. The order number for 
these areas is specified on the lseq card. 


The format of the lseq card is 


I LSEO 


SID 


AID 


LID 


The sid is selected by the loadset = sid in the Case Control Deck. The aid entry 
defines the order of the load vectors representing the total surface area and the absorption 
areas. The lid represents the Set Id of the applied pressure (force) loads. Each lid is 
associated with a specific aid. Therefore, the absorption areas are numbered by their as- 
sociation with the aid entry. 


Caution: The pressure load applied to the surface of the fluid in contact 
with the structure (=> coupling matrix) must have the lowest numbered 
aid. 


Example: 

Let S B be the total (discretized) surface of the fluid, S BC the total surface of the fluid cou- 
pling with the structure and I BC i through J Ba the element numbers of S B . Assume S BC =S B , 
i.e., the whole surface is coupling. Further, let S BA be the total absorbing area of the fluid 
surface with S BA c S B consisting of m different subareas S BA i, S BA2 ,...,S BAm . The ele- 
ment numbers belonging to these subareas are I BA i through J BA1 , I BA2 through J BA2 ,..., 
^BAm through J BAm . 
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PL0AD2 and LSEQ Input 


Coupling surface 

P LOAD 2 

REMM 

1.0 

Em 

THRU 








LSEQ 

SID 

BJSIRIH 

issam 









Absorbing areas 

PLOAD2 

WSEEM 

1.0 

Ibai 

THRU 

Jbai 

P LOAD 2 

LID B a2 

1.0 

IbA2 

THRU 

JbA2 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

P LOAD 2 

■MSB 

1.0 

wmm 

THRU 

■iiif 







LSEQ 

SID 

AIDbai 

BBDSffifBI 



LSEQ 

SID 

WEBSSM 

LID ba2 



♦ 

• 

• 

• 



• 

• 

• 

• 



♦ 

• 

♦ 

• 



LSEQ 

SID 

E"2HI 

E 

> 

3 




In this example each absorbing subarea S BA1 , S BA2 ,... > S BAm is also a coupling area, i.e., it 
is S BAi c S BC (i = 

In this example, the assumption was: S BC =S B . If there are subareas of S B that are not cou- 
pled, then S BC c S B and according to this situation the element numbers, which determine 
the coupling area, have to be changed on the PLOAD2 card for this case. 

Other combinations of coupling and absorbing areas have to be handled in a similar way. 

6.3.4 Area matrix program input/output 

Before any of the response analysis can be performed the AMP must be executed. The 
purpose of AMP is to generate the coupling matrix (or matrices) between the structure and 
the fluid. It also generates the absorption matrix (or matrices) if they are to be included in 
the analysis. 
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6.3.4. 1 Input files 


There are two input files for AMP that by default are defined to be Fortran Units 14 and 15. 
The data which are assigned to these files are generated by MSC/NASTRAN SOL 63 
Normal Modes Analysis of the structure and the fluid. These data consist of the following 
MSC/NASTRAN data blocks: 

OUTPUT2 PFS, CSTM, PGF, SEID, PGA//0/V, Y, CSTMOUT=ll $ 

OUTPUT2 GPLS, SILS, B6PDTS //0/V, Y, GEOMOUT=12 $ 

During the SOL 63 run (Run 1, see fig. 1) five data blocks are written to Fortran Unit 11 
and three data blocks are written to Fortran Unit 12 by default. As can be seen in the above 
DMAP statements these default values can be redefined by specifying parameters CST- 
MOUTandGEOMOUT. 

During the AMP run, data blocks PFS, CSTM, PGF, SEID, and PGA are assigned to For- 
tran Unit 14. Data blocks GPLS, SILS, and BGPDTS are assigned to Fortran Unit 15. 
The contents of these data blocks are as follows: 


Name 

Contents 

PFS 

This is the partition vector that separates the fluid grid points from the 
structure grid points in the residual. 

CSTM 

Coordinate System Transformation Matrices. 

PGF 

This file contains the areas associated with each fluid grid point that is 
on the interface with the structure model. 

SEID 

Contains the Superelement ID's if there are superelements in the analysis. 

PGA 

This file contains the absorption data. Each column of this matrix defines 
a different absorption area. 

GPLS 

Grid Point List Table. 

SILS 

Scalar Index List Table. 

BGPDTS 

Basic Grid Point Definition Table that defines the locations of all the grid 
points in the model. 


The direct input to the AMP program is input at the terminal or Unit 5. The input is rela- 
tively simple and will be described by the order of appearance. The input is assumed to be 
in free format, meaning that the data are separated by a comma or a blank. 
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RECORD NO. 1 - TOL1, TOL2 (real, real) _ 

TOL1 If the fluid and structure grid points are within the distance less than TOL1 

then they are considered to be coincident. 

TOL2 In general, the area matrix terms have comp onents in x,y, and z directions 
(A,, A y ,A t }. The total area then is simply A t = ^A^ + A^ + A^. If for any of the 
area components, AJA t < TOL2 then A i will be ignored, i.e., A i = 0 (t = 
x,y,z), Usually, TOL2 is equal to 0. 


RECORD NO. 2 - IFF, FF (integer, real) 

IFF 

Identifies which of the structural DOFs will be used as the pressure DOF. 

If the Tt DOF is used as pressure DOF, then IFF = i ( i = 1,2,3), e.g., T1 is 
the DOF in x direction. 

FF 

Defines the multiplier for the area matrix terms. Used primarily to reverse the 
sign of the surface normal, which must be inward to the fluid, if one made an 
error. 

Usually, FF is equal to 1.0. 


RECORD NO. 3 - ANS (character) 

ANS Y or N, i.e., Yes or No. If Yes, then the user will be requested in the next 
record to specify the search region for the determination of fluid to structure 
grid points match. In some cases, by specifying the search region, the com- 
putational effort may be reduced. In most cases this has been found not to 
be very significant. 


RECORD NO. 4 - XL,XU,YL,YU,ZL,ZU (reals), required only if ANS=Y 
The above six numbers define a rectangular region in the basic coordinate system. They 
provide the lower (L) and upper (U) bounds for each coordinate direction X, Y, and Z. 


The following records are required only if there are absorption areas. The AMP will print 
out the number, NCOL, of unique absorption areas that it has found from the PGA data 
block. 


RECORE 

NO. 5 - (NR(I), 1 = 1, NCOL) (integers) 

NR(I) 

IfNR(I) < 0 the absorption data will not be generated for I'th region. 

If NR(I) > 0 the absorption data will be generated for I'th region. 

If all NR(I) < 0 no further input is required. 
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RECORE 

NO. 6 - NF (integer) 

NF 

Specifies the number of frequencies at which the absorption data will be pro- 
vided. This is required if absorption is defined as a function of frequency. If 
a constant value of absorption is desired, set NF = 1. 


RECORI 

NO. 7 - (F(I), I = 1, NF) (reals) 

F(I) 

Frequencies at which the absorption data will be provided. Note that the ab- 
sorption data for all retained regions must be supplied for all the frequencies 
even though some regions may not have frequency dependent absorption 
properties. 


RECORD NO. 8 - RHOC (real) 

RHOC Value representing the product of the density, p f , and the speed of sound, c f , 
1 for the fluid medium (see section 4.1, constant Z„ in eq. 93). 


The following NF records are repeated for every absorption area that is active. At most 
this number will be NCOL. Note, that in Record 5 any or all absorption areas can be ig- 
nored. 


RECOftD NO. W 




RECORD NO. 9+(NF-l) - CR(NF), CI(NF) (reals) 


CR(I) 


Real part of the specific acoustic admittance. 


CI(I) 


Complex part of the specific acoustic admittance. 
See eqs. 93, 94, and 95 in section 4.1. 


6.3.4. 2 Output files 

Three files are generated by the AMP program. By default, they are written to Fortran 
Units 16, 19, and 20. Data contained in these files are as follows: 

* UNIT 16: MSC/NASTRAN DMI cards that define the coupling matrices, describ- 
ing the coupling area of the fluid surface. There will be as many different names 
for the DMI matrices as there are superelements that have any internal points con- 
necting to the fluid model. If there are absorption areas, then this file will also con- 
tain the DMIG cards that describe the damping properties for each absorption area. 
All absorption areas are grouped into one DMIG matrix. Since absorption matrices 
can be frequency dependent, there will be as many unique DMIG matrices as there 
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are frequencies at which the the absorption properties are defined. This file is re- 
quired in SOL 70, 7 1 , and 72. 

• UNIT 19: MSC/NASTRAN Case Control Set card definition, 

SET 911 where n{s represent fluid grid points that are on the sur- 

face of the fluid mesh and in contact with the structure. This file will be used in the 
calculation of the Grid Point Participation Factors. 

• UNIT 20: This file contains some DMAP statements. These statements redefine 
the names of the DMI matrices (area matrices) that are stored in the data base. It 
changes the names from that of "DMIxxtx 0 0" to "AREA 0 xxxx," where xxxx is 
some integer value and the other two values are the data block qualifiers. This file 
is required in the SOL 70, 71, and 72. 

6. 3 .4. 3 Absorption property interpolation 

Generally, absorption properties are frequency dependent (see section 4.1). The depen- 
dency is significant and cannot be ignored. To accommodate this phenomenon, a provision 
has been made to allow the user to specify any curve that defines the properties as a func- 
tion of frequency. There is a capability in the FRRD2 module (Frequency Response for 
Aeroelastic Analysis) that uses a spline technique to interpolate any complex matrix quan- 
tity. The user is required only to specify matrices of complex quantities as a function of 
frequency. This capability originally was intended for the interpolation of aerodynamic ma- 
trices. Owing to this fact, a very special form of spline interpolation is used. The user 
must understand this spline form, so that the absorption function can be properly defined. 

This interpolation technique is best illustrated in figure 5. 



0 1 2 3 4 5 6 


Figure 5(a). Interpolation technique. 
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In this figure a simple line, defined by two points, is shown. The absorption is a complex 
quantity and this line represents the same variation versus frequency for both the real and 
the imaginary components. Figure 5(b) illustrates the behavior of MSC/NASTRAN spline 
interpolation. 



Figure 5(b). Behavior of MSC/NASTRAN spline interpolation. 

For the real component the zero slope at zero frequency is assumed, i.e., the curve is sym- 
metric about the y-axis. For the imaginary component the assumption is made that the 
function value is zero at zero frequency, i.e., the curve is antisymmetric about the y-axis. 
Obviously, the resultant interpolation of the real and complex quantities does not produce 
the correct result. To remedy this problem, additional input is required. This problem can 
be avoided by specifying a slope of the curve at the first frequency value. The slope can be 
defined implicitly by specifying two points spaced closely together as shown in figure 5(c). 
The resultant interpolation is shown in figure 5(d). 
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0 1 2 3 4 5 6 

Figure 5(c). Slope defined implicitly by specifying two points spaced closely together. 



Figure 5(d). Resulting interpolation. 


Caution: The absorption should not be defined at zero frequency because 
MSC/NASTRAN divides the function value by frequency when interpolat- 
ing the imaginary part of the absorption quantity. 


6.4 Normal Modes Analysis of the Modal Coupled Fluid-Structure System 

(Run 2) 

This analysis is performed only if the coupled Fluid-Structure-Interaction (FS1) modes are 
required. In this case both the structure and the fluid are represented by their respective 
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modal coordinates, calculated during Run 1 (SOL63). This is, compared to Run 1, a rather 
quick and simple calculation, since the problem size is generally very small. Even though 
SOL 70 is used, the DMAP program allows the user to calculate the real eigenvalues of the 
coupled system if there are no actual damping matrices. In that case a different symmetric 
formulation is used where the coupling terms are no longer in the damping matrix (see sec- 
tion 3.4,, e.g., eq. (64)). Before executing the MSC/NASTRAN analysis the AMP must be 
run in order to prepare coupling (DMI input) and also, if required, absorption (DMIG in- 
put) data. The input to AMP for this case is 


RECORD NO: 1: TOL1 = 1., TOL2 = 0. 

RECORD NO: 2: IFF = 1, FF = 1. 


RECORD NO: 3: ANS = N 

RECORD NO: 5: NR(I) = -1 (I=l,...,NCOL) means: 

-1,-1,...,- 1 (NCOL "times" -1 or NCOLM) 


The Record 5 implies that no absorption matrices will be generated. Input listing of Run 2 
is shown in Appendix B2. 

6.4.1 Executive Control Deck 

1 . It is required that the downstream data base from SOL 63 be assigned to the read- 
only data base, DBSET 15. 

•NASTRAN DBSET 1= (DB01 , DB15) , DBSET 2= (DB01 , DB15) 

• NASTRAN DBSET 15=(DB15) 

2. Include a small DMAP program (DMAP.DAT) generated by AMP, fetched from 
UNIT 20. This usually is done by a merge utility, which varies with the computer 
system used. 

3. Include one of two DMAP alter programs; one is used for the real and the other is 
used for the complex eigenvalue analysis. The real eigenvalue analysis capability 
has been provided for the case with no damping terms. 

Both DMAP programs are inserted following the SOL 70 call. This instruction is, in gen- 
eral, automatically accomplished by the actual machine-dependent procedure and needs no 
additional user action. 
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6.4.2 Case Control Deck 

1. Specify set xx = 0, semg = xx. 

2. Select the usual output requests, i.e., structural or fluid grid points for, e.g., plot- 
ting deformed shapes. 

3. Select eigenvalue extraction method, either method = yy or cmethod = zz, 
depending on the selection of real or complex eigenvalue analysis. 

Where xx, yy, and zz are some user-specified integer numbers. 

6.4.3 Bulk Data Deck 

Specify the following parameters: 

1. PARAM, DLOAD , -1 

2. PARAM, RE SDUAL , -1 

Caution: see parameter description NASTRAN USER’S MANUAL Vol n. 

3. param, lmodes, xx; where xx > sum of modes for the fluid and the structure 
from Run 1 (SOL 63). param, lfreq and param, HFREQ must not be used. 

4. param, small, EPS; where EPS is used to filter the modal area matrix terms. If 
any area matrix term, A tj , divided by the largest L4,yl is smaller than EPS, then that 
Ajj is discarded. If (A, y .It. A ijmax ) A tj = 0. Default for EPS is l.E-6. 

Additionally, the user must specify the eigenvalue extraction method, either EIGR or eigc 
depending on the type of solution. Using the computer system merge utility, include the 
contents of file on UNIT 16. This file includes the area matrix defined on DMI/DMIG 
cards. 


6.5 Modal Frequency Response Analysis (Run 3.1) 

This is probably the most commonly used analysis technique forFSI calculations. Like the 
coupled modes solution, this analysis is also restarted from SOL 63 data base (Run 1). Be- 
fore executing this analysis, the AMP must be run to generate the required files (DMI, 
DMIG, DMAP.DAT) as described in section 6.4. The input to AMP for this case is, e.g.. 


RECORD NO: 1: TOL1 = 1., TOL2 = 0. 

RECORD NO: 2: IFF = 1, FF = 1. 

RECORD NO: 3: ANS = N 
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NO ABSORPTION 

NR'(I) = -1 (r=l,...,NCOL) means: 

-1,-1,...,- 1 (NCOL "times" -1 or NCOL*- 1) 


ABSORPTION: 

NR(I) = 1 (1=1 NCOL) means: 

1,1,..., 1 (NCOL "times" l)orNCOL*l 


Input listing of Run 3.1 is shown in Appendix B3.1. 


6.5.1 Executive Control Deck 

1 . It is required that the downstream data base from SOL 63 be assigned to the read- 
only data base, DBSET 15. 

• NASTRAN DBSET 1= (DB01 , DB15) , DBSET 2= (DB01 , DB15) 

• NASTRAN DBSET 15=(DB15) 

2. Insert a small DMAP program (DMAP.DAT) generated by the AMP. 

3 . Insert a DMAP Alter Program for Frequency Response. 

Both DMAP programs are inserted following the SOL 71 statement. This statement is, in 
general, satisfied automatically by the machine-dependent procedure and needs no addi- 
tional user action. 

6.5.2 Case Control Deck 

No- special instruction must be observed. 

Exception: If pressure response curves at fluid grid points GPFi (t = 1,2,...) are 
required, the user must use the following Case Control cards: 

SET pp = GPF/,GPF2,... 

VELO = pp 


Plot request: 

XYPLOT VELO / GPFI (TlRM) . 

Where pp is some user specified integer number. 
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6.5.3 Bulk Data Deck 


Specify the following parameters: 

1. PARAM, DLOAD, -1 

2. param, norf, -l; no residual flexibility for the fluid. 

3. PARAM, LMODES, XX; see section 6.4.3) 

4. param, SMALL, EPS; see section 6.4.3) 

There are two possible methods for specifying the modal damping for the structure and 
fluid. The first method is to use the TABDMP1 Bulk Data card. In this case, the fluid 
modal damping is defined in the third quadrant, i.e., frequency and damping values must 
be input with a negative sign. The second method is to input the modal damping via the 
DMI cards. Damping can be specified for each structure and fluid mode, as output from 
Run 1, via DMI cards. One DMI card is assigned to each mode that has modal damping. 
The number of the DMI card must be the same as the number of that mode. Using the DMI 
input, the DMAP program generates modal damping matrices TGDMPS and TGDMPF for 
structure and fluid respectively, which are added to the total modal damping matrix of the 
modal coupled system. 

The following example illustrates the modal damping input via DMI cards: Let n s be struc- 
tural modes and nj- fluid modes calculated in Run 1 . Some of the modes are damped and 
some are not damped. 

Fluid: 


DMI, TGDMPF, 0, 2, l,l„l,n/ 



Mode No. 


Modal Dmp. 

DMI 

TGDMPF 

1 

i 

0.03 

DMI 

TGDMPF 

2 

i 

0.03 

DMI 

TGDMPF 

3 

i 

0.03 

• 

♦ 

• 

• 

• 

• 

• 

• 

• 

• 

DMI 

TGDMPF 

n 

i 

0.03 


- 57 - 




Structure: . 


DMI, TGDMPS, 0, 2,1, 

* 



Mode No. 


Modal Dmp. 

DMI 

TGDMPS 

1 

i 

0.06 

DMI 

TGDMPS 

2 

i 

0.03 

• 

* 

• 

• 

• 

• 

• 

• 

• 

• 

DMI' 

TGDMPS 

9 

i 

0.0 

DMI 

TGDMPS 

10 

i 

0.02 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

DMI 

TGDMPS 

21 

i 

0.0 

DMI 

TGDMPS 

22 

i 

0.04 

• 

• 

• 

• 

• 

• * 

• 

• 

• 

• 

DMI 

TGDMPS 

n 

i 

0.0 


Caution: The modal damping is input as % of critical modal damping, 
where 3% is specified as 0.03. 


6.6 Modal Transient Response Analysis (Run 3.2) 

In many respects, this analysis is very similar to the frequency response analysis so it will 
not be described in detail. The difference in the two analysis methods are the same as for 
any other: non-FSI analysis, i.e., FREQ versus TSTEP, RLOADi versus TLOADi, etc. 
Furthermore, a special DMAP Alter Package is used and no absorption damping is al- 
lowed. 


6.7 Structural Grid Point Participation (Run 4) 

Before executing this run, the frequency response analysis, Run 3.1, must have been suc- 
cessfully completed. The output from this run is a very useful diagnostic tool in the acous- 
tic analysis. With this output the user can pinpoint the source that causes a pressure rise or 
fall at any point in the interior of the fluid domain. Since one is usually concerned with the 
structure surrounding a fluid, this information can be beneficial in the redesign of the part 
of the structure that causes a noise problem. 
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6.7.1 Executive Control Deck 


1 . It is required that the downstream data base from SOL 63/SOL 71 be assigned to 
the read-only data base, DBSET 15. 

NASTRAN DBSET 1= (DB01, DB14, DB15) 

NASTRAN DBSET 2= (DB.01, DB14 , DB15) 

NASTRAN DBSET 15= (DB14, DB15) 

2. Insert a DMAP alter program 

The Grid Point Participation DMAP Alter Program is automatically inserted by the machine 
dependent procedure and no user action is necessary. 

6.7.2 Case Control Deck 

1. Insert the file SET.DATA (SET 911) from UNIT 19, which contains the fluid 
points that are coupled to the structure. 

2. disp = 911 

3 . Input of frequencies for which the participation factor should be calculated 

SET XX = Fl,F2,F3, . . ,Fn 

OFREQ = xx 

The input frequencies Fl,F2,F3, . . ,Fn should be approximately equal to the 
frequencies used in Run 3.1 (SOL 71). 


6.7.3 Bulk Data Deck 

1 . User Instruction for calculating the participation factor at fluid grid point GID/ (cid 
is the DOF which determines the pressure; usually it is CID = l (see description of 
Run 1). 

DMIG, PDOF , 0, 6, 1, 0 

DMIG, PDOF, GIDi, CID, , GIDi, CID, 1 . 0 

2 . If one only wants the real displacement output set 
P ARAM, REAL, 0 
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6.7.4 Output 


The output from this run requires some explanation because it is not immediately obvious. 
This output describes the contribution of the structural grid points to the generation of the 
pressure anywhere in the interior of the fluid domain. Obviously, the only structural grid 
points that can contribute to the pressure in the fluid domain are the points on the interface 
to the fluid. The AMP establishes the relationship between the fluid and structural points 
on the interface. In general, one structural grid point is "connected" to a closest single fluid 
grid point. In other words, there is a one-to-one correspondence between the fluid and 
structural grid points on the interface. The goal of this output is to provide a numerical as 
well as a graphical representation of structural grid participation factors. 

In most typical applications of this acoustic procedure (such as aircraft or automobile noise 
studies), the structural mesh will usually be much finer than the corresponding fluid mesh, 
therefore, only some of the structural points on the interface will be "connected" to the fluid 
surface. Also, the structural mesh on the interface may contain a number of different types 
of elements. It is very difficult to graphically (contour plots, etc.) present the participation 
information on such a mesh. The fluid mesh is usually composed of three-dimensional el- 
ements and it is easy to plot contours on the outside surface of the fluid domain. To facili- 
tate the graphical display, the output of the structural grid participation factors is transfered 
to the corresponding fluid grid points as produced by the AMP. The results of this opera- 
tion can be seen in section 8. 

The output of this analysis looks like the usual SORT1 "Complex Displacement Vector" 
output for frequency response. It is generated for a user-specified interior fluid point, 
GID«, at each frequency, F, (/ = 1,2,. ..,NF), input in the Case Control Deck (see above). 
The contribution to the pressure at interior fluid point, gid /, influenced by the participation 
points, is given. The participation points are the structural points on the interface, but in 
the output they are labeled by the corresponding fluid point on the interface as explained 
above. 


COMPLEX DISPLACEMENT VECTOR 
(REAL/IMAGINARY) 

Participation 
Point ID 

T1 

T2 

T3 

R1 

PGP; 

RE(T1) 

RE(T2) 

RE(T3) 

R1 

0'=1,2,.,N C ) 

IM(T1) 

IM(T2) 

IM(T3) 

0 

• 

• 

• 

• 

• 

• 

• 

• 

♦ 

• 

• 

• 

• 

• 

• 
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Nc is the number of fluid coupling points on the interface and they are listed in SET 911. 

1 . Listed under the heading T1 is the real (RE(T1)) and the imaginary (IM(T1)) com- 
ponents of all participating points PGP, (t=l,...,N c ). If all the terms are added, the 
actual pressure at frequency F, for the selected fluid point, gid; , cid, will be ob- 
tained. 

2. Under the heading T2 each term in T1 is rotated in the complex plane such that the 
real term is in the direction of the resultant pressure vector (for point GID/) and the 
imaginary term will be at 90" to that vector (see fig. 3). The new coordinate system 
is labeled by ( ’ ). 



Note the sign of RE(T2): ± means "contribute" to pressure p 

- means "subtract" from pressure p. 


3 . Under the T3 heading the output is simply the vector 



which gives the fraction of the contribution for each point. Note that only the real 
term of T2 contributes to the resultant pressure. iPjl is the magnitude of the pres- 
sure vector. 


4. Under the R1 heading, the output is the same as the T3 except that each term is 
normalized by the largest real term. 
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7. DEMONSTRATION OF ACOUSTIC CALCULATIONS AND 

THEORETICAL VERIFICATION 


7.1 A One-Dimensional Example 
7.1.1 Description of the model 

Generally, it is very difficult to obtain theoretical solutions for practical problems of FSI. 
Such solutions are necessary to evaluate the results and the accuracy of the finite element 
analysis. A simple Fluid-Structure System for which an analytic solution exists will be 
used to verify the finite element results and to demonstrate most of the capabilities of the 
acoustic procedure presented here. Owing to his expertise in the mathematics of acoustics, 
all the equations in this section were derived and solved by the second author. 

The system consists of a straight air-filled tube with quadratic cross section. One end of 
the tube is closed by a rigid wall (|Z„| = <»), or by an absorbing wall (0 < |Z„| < °°), and the 
other end has a moveable piston supported by a spring-damper-system. Z„ is the normal 
impedance of the boundary. Additionally, the piston surface can be chosen to be absorb- 
ing. A more general situation is to introduce flexible walls on both ends of the tube (with 
or without absorbers) which can be done by a simple spring-damper system. 

The natural frequency of a one-degree-of-freedom spring-damper-piston system is assumed 
to be close to that of the air with a fixed piston in order to create an interaction between the 
air and the structure. In view of response calculations, the piston can be loaded by an ex- 
ternal force. The model is shown in the figure 6. 
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The following table gives the physical properties which are necessary to describe the cou- 
pled system. 


Table 2. Physical properties j 

Spring-damper-piston system (Structure) ! 

area of piston cross section 

A 

0.000625 m 

spring stiffness 

K, 

7474.75 N/m 

mass of piston 

K 

0.01 kg 

viscous damping 


0.50 kg/s 

external load 

r. 

2.1885 N 

| Air-filled tube (Fluid) | 

length of tube 

i 

1.25 m 

area of tube cross section 

A 

0.000625 m 

speed of sound (air) 

c f 

344 m/s 

density of air 

Pf 

1.205 kg/m 

bulk modulus 

X 

142594.17 N/m 

specific acoustic resistance 

z, =ReZ„ 

228 kg/m 2 s 

specific acoustic reactance 

Z 2 = Im Z„ 

-1456 kg/m 2 s 


Remark: Because / » 4a , it is sufficient to describe the wave propagation in the tube by a 
one-dimensional wave equation along the longitudinal axis, x-axis, of the tube. 

7.1.2 Theoretical solution 

The equation of motion of the damped spring-piston system is given by 

Mfi(t)+D,u(t) + K,u(t) = F(t) (122) 

M t> D s , and K t are mass, viscous damping, and stiffness of the system. The displacement 
of the piston in Jt-direction is u. The exciting force, F(t), consists of an external part, F s (t)> 
acting on the piston and a part, F p (t), induced by the fluid pressure on the piston surface. 
Discussion is restricted to the steady-state case, i.e., piston displacement, pressure, and 
exciting force are all assumed to be harmonic in time. Therefore, the following boundary 
value problem, equation (123), has to be solved: 
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Structure 

Fluid 

(K - (0 2 M;+i(oD,)u = F S + F p 
with' F p = -Ap( 0). 

{ 

< 

i 

^p(x) + k 2 p(x) = 0 (k = co/c f ) 

Boundary Conditions: 

1. * = 0: 
dp 2 

. .dp f . i V 

2. x = l: -Z- = a -mpf—p , 

V L n ) 

with a = 0 : rigid, a = 1 : absorbing. 

Z n is the impedance of the absorber (see 
section 4.1) 


In the following, several formulas will be presented for 

1. coupled eigenvalues 

2. pressure distribution in the tube along the x-axis belonging to these eigenvalues 

3. expressions for the piston displacement, u = u(v), and pressure, p(x,v), in the 
tube at location x 

All these formulas will depend on the excitation frequency. The derivation of the formulas 
is elementary, therefore, details can be omitted. A more detailed mathematical investigation 
of some FSI questions can be found in reference 33. 

7. 1.2.1 Coupled Eigenanaly sis 

With D, = 0, F, = 0, and dp/dx = 0 at x = /, from equation (123) follows the eigenvalue 
problem of the the spring-piston-tube system. In order to obtain the eigenvalue equation, 
first obtainrsolutions of the fluid boundary value problem of the form 

p(x,k ) = C,cosfcc + C 2 sin fee 

where C, and C 2 have to be determined from the boundary conditions, elastic at x = 0 and 
rigid at x = /. 

Using k = to jc f , co = 2rcv, and Z 0 = c f p /? the following expression is obtained 


p(x,v)=2nZ 0 uv 


2* ) 

-1 

. . \ 

2nl 

f . 2 nv 


r 2 nv 

tan — v 

< 

tan v 

sin x 

+ 

cos X 

c / J 


c f J 

c f J 


l c f )\ 


(124) 
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and from that 


p(0,v) = 2nZ 0 uv 


' 2 nl V ’ 

tan v 

\ c f J 


(125) 


Combining the structural part of equation (123) with equation (125) yields a transcendental 
equation for the frequencies, v= v- (j = 1,2,3,...), of the coupled system, 


2 nl ZA 
tan v = — — 


c f 2 nM. v 2 -v. 2 


(126) 


where y, is the frequency of the spring-piston system (see eq. 127). Table 3 shows the 
frequencies of the uncoupled and coupled structure-fluid system. The calculation is based 
on the data shown in table 2. 


Table 3. Coupled and uncoupled frequencies (in Hz) 

Spring-Piston 

Tube 

Coupled 

System 

(air-filled) 

System 

[ 137.6 

137.6 

128.3345 



147.1835 


275.2 



412.8 



550.4 



Note that the frequencies of the spring-piston system and the air-filled tube with rigid 
boundary at x = 0 and x = / were calculated from 


v = 



and v i " ) = n 2 i (" = 1> 2 » 3 »--) 


(127) 


It can be seen from table 3 that the interaction between the spring-piston system and the air 
in the tube, with rigid boundary at x-l, leads to a considerable shift in the frequency, v,, 
of the spring-piston system and of the fluid. Clearly, as j increases, the frequency 
shift, Vj - 0 = 3,4,5,...), tends to zero. Moreover, the coupling increases the fluid 

frequencies, i.e., v ( j~ x) < v. if j > 2. 

In figure 7, equation (126) is solved graphically by the intersection of two curves 

yyv)=uJ^v = (128) 

Cf 2 JlM s V — V s 
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in the case v, = v^. This clearly shows the effect of coupling described above. 



Figure 7. Graphical determination of coupled frequencies. 

In the uncoupled case, i.e., rigid boundary at x - 0 and at x = /, the stationary pressure dis- 
tribution is given by 


p“(x) = const, cos 




n 


= 1,2,3. 


(129) 


In comparison with equation (129) for the coupled case, the following pressure distribution 
is obtained (see eq. 124): 


Pj(x) = const. Vj 


- , ^ 

2nl 

-1 

2 nl \ 

( 2 nv ' 


f 2wv )\ 

tan v. 

A 

tan v, 

sin L x 


cos L x • 

c / J 


c / J 

c f J 


l c f )) 


\J = 1,2,. (130) 


As above, v l f n) ( n = 1,2,3,...) are the frequencies of the standing waves in the rigid 
bounded tube and v y (j = 1,2,3,...) represent the coupled eigenfrequencies. Both groups 
of frequencies are listed in table 3. 

Because of condition 1. in equation (123), the coupling also changes the pressure gradient 
at x = 0. For the first two coupled modes 
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4a <o, & >o (i3i) 

dx i= 0 dx x =o 

In comparison, for the uncoupled case at the rigid boundary, x = 0, the pressure gradient is 
null. 

It is obvious that the position of the pressure node, i.e., values for which the pressure is 
zero, is also influenced by the coupling. These values can be calculated from equation 
(130) and compared to the corresponding pressure node positions of the standing waves in 
the tube with rigid boundaries at x = 0 and x = /. There is a significant shift of the pressure 
node position, x) (j = 1,2), especially for coupled modes 1 and 2, compared to 1/2, which 
is the pressure node position of the first standing wave in the rigid tube. These effects are 
shown in figure 8. 



Figure 8. Pressure distribution, p ; (x), for coupled modes j =1,2, compared to 
PT(x), of the first standing wave (n = 1) in a rigid tube. 


Coupled case 

:(•) 

£ <0, 
dX ,=o 

(..) f >0 

dx j-o 

Uncoupled case ("uc") 

: (*♦*) 

is. =0 

dx x=0 
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7. 1 .2.2 Frequency response analysis 


To get the response (piston displacement u and pressure p) of the coupled system, the 
boundary value problem, equation (123), has to be solved. Starting with a formulation for 
the pressure function, equation (124), an elementary calculation yields u( v) for the piston 
displacement and p(x, v) for the pressure in the tube. The following two expressions re- 
sult: 






and 


p(x, v) = 2 nZ v 


/«. 




( 


ft 


_ fjM. 

„ 2 nv .2 nv 

F cos x + sin ; 


(132) 


A 


V 


7 ) 


. 2nv 
+ q.T. cos x 


-H 


{r Inv c 
^ r 2 cos x - <*. 


V 


2nv . 2 nv 
r cos x + sin x 

Cf C f ) 


respectively, where 


^ v ‘- v >T%r r ' v 




D+AZJ \ 


•o'- 2 


2nM s 

f — fi- 

Jm * 4?r 2 M, 


r,= 


r,= 


1 r, it 2 2 \i • 4 nl 4 nl 

-[1 - a [a + y ) J sin v+ acrcos v 


f \ 2 

1 . 2 nl 2 nl 

sin + aa cos v 


v */ 


2 2 

+ ay 


' 2 nl V 

cos V 


7 ) 

ay 


2 ~/ 


. 2nl 2 nl 

sin — + aa cos — 


V <7 


7 ) 


\2 

( ~ t \ 2 

2 2 

2tt/ 

+ a y 

cos V 

J 

^ c / > 


(133) 


(134) 


(135) 


a = 


j 0 : boundary at x = l is rigid 

[l : boundary at x = / is absorbing 


a— 36. 


r - Z o h 

Z 2 +Z 2 ’ Z 2 +Z 2 " 

Z 0 = p f c f , Z, = ReZ„, Z 2 = ImZ n 


(136) 
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It is useful to present the response results, i.e., piston displacement, w(v), and pressure, 
p(x, v), in the tube, in terms of absolute values, |u( v)|, | p(x, v)|, and of their phase shifts, 
cp s (v) and (p f {x, v), due to the exciting force 

F,{t) ~ F s ' exp(t • 2nvt) 


From equation (132) one obtains 


and 


|m(v)| = 


Sm, 


(«?+«) 

tan <p.(v) = -— 


1/2 




(137) 


(138) 


Analogously equation (133) yields 


\p{x,v)\ = 2nZ 0 v 


Si 


M. 






\ 


-|2 


_ 2nv .2 JIV 

i cos x + sin x 

v c f c f ) 


, e T' 2nv 

+ $,r cos — x 


2^V 


^r 2 cos x - % 


( 2 nv . 2nv N 

T cos x + sin x 

v c / c / yj 


>V2 


(139) 


and 


e ^ 2?TV , 

cos x - 


tan <p f (x, v) = 


„ 2;rv . 2nv ^ 
i,cos jt + sin x 

c f c f ) 




_ 2nv . 2;rv 
Fj cos x + sin x 

c s c f J 




2 nv 


(140) 


The solution for a rigid boundary at x = / can be obtained by setting a = 0 (see eq. (136)) 
in equations (134) and (135). Then the coefficients T, and T 2 are reduced to 


r, = l/tan— v, r 2 = o 


(141) 


It is straightforward to derive the response formulas for this case from equations (137) 
through (140) using equation (141); therefore, the details will not be shown. 

The following figures include the graphs of the expressions for the piston displacement 
|u( v)|, the pressure \p(x, v)|, dependent on v with fixed x. Furthermore, the phase shifts 
<p t ( v) and (p f (x , v) are also shown. 
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Displacement (nun) 



Figure 9. Piston displacement. 


Phase angle (Deg) 



Figure 10. Phase shift (p t (v ) of the piston motion. 
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Pressure (N/mm 2 ) 



Figure 11. Pressure response curves at x = 0 (i.e., piston surface). 
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Figure 12. Pressure response curves at* = 0.15 m. 
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Figure 13. Phase shift (p f (x, v) of the pressure response. 


Figures 9 through 13 exhibit the characteristic phenomena caused by the interaction of the 
piston and fluid in the tube. The discussion of these phenomena is restricted to the follow- 
ing remark: In the frequency range from v^/2 to 3 v^/2 , neglecting absorption at x = /, 
from figure 10, the phase shift, (p,(v), for the coupled case jumps from 180° to 0° at 
v = v^. On the other hand, figure 9 shows that u( vjP| = 0. That means the piston does 
not move if the excitation frequency is v = vj?* . This effect occurs because the pressure-in- 
duced force, F = -Ap{ 0, v) (see eq. (123)), is moving in the direction opposite to the ex- 
citing force, F s (t) = F s zxp(i27ivt) and IfJ = \F s (t)\ if v = . 


The pressure on the piston surface depends on the piston displacement expressed by the 
following equation: 


p(0, v) = 2tcZ 0 T x v u( v) (142) 

which follows from equation (133) with x = 0, T 2 = 0 and equation (132). The factor, Tj, 
given by equation (141), in equation (142) causes p(x, 0) 4T u(v) for v*p/2 < v < and 
p(x,0) i-i u(v) for v| !) < v<3v| 1) /2. For the first two coupled frequencies the terms 
m( v) and p{ 0, v) are purely imaginary. An analogous statement is valid for x > 0. 

Table 4 shows the phase relation between exciting force, F s (t) at t = 0, the piston dis- 
placement, u( v), and the pressure, p( 0, v), on the piston surface. 
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Table 4. Phase relation between F s (t), n(v),and p( 0,v) 


■ 

/ 2 < v < v. 

v= V, 

V, < V < V { p 

v=vp 

Re u 

>0 

=0 

<0 

=0 

Im u 

<0 

<0 

<0 

=0 

Re p 

<0 

=0 

>0 

=0 

Im p 

>0 

>0 

>0 

=0 


\J 



L 

p 


\4 


_P_ 


\ 



~~7 

U 


p=o 

u=0 


v? < v < v 2 

v= v 2 

v 2 < v<3vJ 1) /2 


Re u 

>0 

=0 

<0 


Im u 

<0 

<0 

<0 

Re p 

>0 

=0 

<0 

Imp 

<0 

<0 

<0 





p 


p 



1* 

r 

'P 




In closing, the pressure distribution, \p(x, v)|, in the tube as a function of jc at selected fre- 
quencies is presented (see figs. 14 and 15). There is a significant difference in the location 
of the minimal pressure which depends on parameter v . Positions, x, can be determined by 
using equation (139). 
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Pressure (N/mm 2 ) 



Figure 14. Pressure in the tube for ( v = 100, 128.33, 137.6, 147.2, and 200 Hz) 
without absorption at x = /. 


Pressure (N/mm 2 ) 



Figure 15. Pressure in the tube for ( v = 100, 128.33, 137.6, 147.2, and 200 Hz) 
with absorption at x-l. 

7. 1 .2.3 Transient response analysis 


This type of response analysis is more complicated than the one investigated in section 
7. 1.2.2. The mathematical description of the solution procedure requires significant elabo- 
ration, therefore it will not be presented. 
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7.1.3 Finite element approximation and comparison with theoretical results 

The coupled fluid-structure system introduced in section 7.1 is approximated by a suitable 
finite element model, which can be a one-, two-, or a three-dimensional. The theoretical de- 
scription and solution of the coupled problem is based on the assumption that the coupled 
system is one-dimensional. 

This assumption is sufficient because the condition, l »y[A, where / is the length of the 
tube and A is the area of the cross section, is fulfilled. Therefore, it is appropriate to de- 
scribe the coupled spring-piston-tube system by a one-dimensional finite element model. 
In reality most problems are three-dimensional. Consequently, a three-dimensional model 
is used for the FE approximation. It is left to the reader to describe the coupled system by a 
one-dimensional finite element model. 

In the following sections, the finite element model of the spring-piston-tube system is de- 
scribed first. Subsequently, the finite element results for the usual types of analysis are 
presented. These analyses consist of 1) uncoupled and coupled eigenanalysis of the spring- 
piston-tube system, and 2) frequency response analysis. In many cases, transient response 
analysis can be performed by transforming the problem to the frequency domain via the 
Fourier transform. Then the same solution scheme can be used. 

The finite element results were calculated by applying modal and direct (or physical) cou- 
pling of structure (spring-piston system) and fluid (air in the tube) and compared to theoret- 
ical results, which were presented in section 7.1.2. Equations (67) and (68) were used for 
the modal coupling and equations (62) and (64) for direct coupling (see section 3.4). 

A detailed description of the NASTRAN runs can be found in Appendix B. NASTRAN 
user decks (Executive, Case and Bulk Data Deck) are shown for all runs. The background 
can be found in section 6 (sections 6.3 through 6.7). Various cases of boundary absorp- 
tion input data are presented. 

7 . 1 . 3 . 1 Finite element model 

The fluid in the tube is modelled by 800 8-node hexa elements. In this case, they are the 
NASTRAN CHEXA8 elements. The appropriate FE mesh consists of 1275 nodes. 
Boundary conditions need only be introduced at the ends of the tube. At the rear end, the 
boundary is either a rigid wall or an absorbing wall. Only the latter case needs an addi- 
tional input. The absorbing wall is defined in the same way the coupling surface of the 
fluid (air) is defined, by 16 NASTRAN CQUAD4 plate elements. Both surfaces belong to 
the fluid model. The piston, combined with a spring, represents a flexible wall (coupling 
surface) at the front boundary of the tube. The finite element model of the spring-piston 
system consists of a central node with a CONM2 element, representing the mass of the 
piston, and connected to ground via the CELAS2 element. The remaining 24 nodes define 
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the piston surface and are connected to the center point by MPC conditions. Damping is 
provided by a CDAMP2 element. 


Viscous damper 
(CDAMP2, 4 ) 


Absorbing wall 


Coupling wall 




Fluid-filled tube (p^. ,c^ ) 


Piston (CONM2, Mg ) 
with rigid surface 
using MPC's 


Spring (CELAS2, K s ) 


Figure 16. Finite element model of the spring-piston-tube system. 

7 . 1 . 3 . 2 Uncoupled and coupled eigenanalysis 

NASTRAN input decks are shown in Appendixes Bl, B2. In tables 5 and 6 the analyti- 
cally computed frequencies of the structure (spring-piston system without damping) and 
fluid (air-filled tube) are compared to FE results. 


Table 5. Eigenfrequenc 

sysl 

ies of the spring-piston 
tern 


1f(Hz> 

M s 

(Hz)/ FE (NASTRAN) 

1 137.6 

137.6 
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Table 6. Frequencies of the first four standing 
waves (longitudinal) 

v ( /° -n— (Hz) 

f 21 

v f* (Hz)/NASTRAN 

137.6 

137.6226 

275.2 

275.3812 

412.8 

413.4115 

550.4 

551.8498 


Table 7 demonstrates the variation of frequencies with respect to different FE-model dis- 
cretization. Frequencies of the first four longitudinal waves versus four different mesh 
sizes are shown. The analytic values are listed once again. 


Table 7 Frequencies of tubes with different mesh sizes 

| 

275 nodes 
160 hexa 

650 nodes 
400 hexa 

1275 nodes 
800 hexa 

2525 nodes 
1600 hexa 

137.6 

138.1667 

137.6907 

137.6226 

137.6057 

275.2 

279.7465 

275.9251 

275.3812 

275.2452 

412.8 

428.1957 

415.2490 

413.4115 

412.9529 

550.4 

586.9007 

556.2116 

551.8498 

550.7625 


Table 8 shows the frequencies of the coupled system calculated from equation (126) and 
from the finite element model using modal and direct coupling between the structure and 
fluid. 


Table 8. Frequencies of the coupled system 

v in Hz 

v in 1 

tfz/FE 

"theoretical" 

direct coupled 

modal coupled 

128.3345 

128.3430 

128.4463 

147.1835 

147.1987 

147.3315 


276.2529 

276.2570 


413.9049 

413.9070 


552.2022 

552.2039 
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As mentioned in section 7.1.2, the coupling shifts the pressure node position and changes 
the pressure gradient at x = 0 compared to the uncoupled case (rigid boundary at x = 0), see 
figure 8. These effects can also be shown by using the finite elements to describe the dy- 
namics of the coupled system. The quality of this approximation depends on the mesh size 
of the FE-model as well as the type of coupling. 

Table 9 shows the dependence of the coupled frequencies on the mesh size of the FE- 
model. In this case only the mesh size of the fluid model was varied. Furthermore, the 
computations were carried out using only modal coupling with four fluid modes. 


Table 9. Air-filled tube coupled frequencies versus mesh size 


"theoretical" 
Results (Hz) 

Finite Element Results (Hz) 


275 nodes 
160 hexa 

650 nodes 
400 hexa 

1275 nodes 
800 hexa 

2525 nodes 
1600 hexa 

128.3345 

128.6230 

128.4690 

128.4463 

128.4406 

147.1835 

147.6648 

147.3728 

147.3315 

147.3212 


280.6546 

276.8048 

276.2570 

276.1200 


428.7420 

415.7506 

413.9070 

413.4470 


587.3251 

556.5737 

552.2039 

551.1145 


Assuming that the mesh size of the FE-model is "good enough," the best approximation is 
obtained by "direct coupling," which means that the structure and fluid are coupled in the 
physical coordinates before any modal transformation, see section 3.4. In realistic prob- 
lems, like automobile interior acoustics, where the FE-models of a car body and the air in 
the passenger compartment are very large, modal reduction must be applied, (see also sec- 
tion 3.4), to reduce computation time. As a consequence of this action, the quality of the 
coupling approximation is reduced. The analyst has to select the number of structure and 
fluid modes needed for the analysis. That is, how many modes are required to obtain a 
sufficiently accurate approximation for the physics of the coupled problem? The user 
should be careful at all times in applying modal reduction. 

The result of direct and modal coupling for the pressure distributions corresponding to the 
first two coupled modes is shown in figures 17 and 18. For the case of direct coupling, the 
deviation from the theoretical result (see figs. 17 and 18) is imperceptible. In comparison, 
the result based on modal coupling with a reduced fluid modal matrix shows a considerable 
deviation. 
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Figure 17. First coupled mode pressure distribution computed by using direct and modal 
coupling. 



Figure 18. Second coupled mode pressure distribution computed by using direct and modal 
coupling. 

7 . 1 . 3 . 3 Frequency response results 

This section presents the results for the boundary value problem, equation (123). The re- 
sults were calculated using the finite element method and compared to analytic results (see 
figs. 19, 20, and 21). See Appendix B3.1 for the NASTRAN input of the frequency re- 
sponse run. 
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Figure 19. Piston displacement response. 
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Figure 20. Pressure response at x = 0. 
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Figure 21. Pressure response at x = 0.15 m. 


7. 1.3.4 Acoustic contribution analysis 

During this analysis the influence of the fluid surface, which is coupled to the surrounding 
structure, is determined. This is performed in Run 4, which is a restart of Run 3.1. 

A description of the NASTRAN input deck can be found in section 6.7. The NASTRAN 
input for Run 4 of the coupled system, presented in figure 16, is listed in Appendix B4. 
Contribution of fluid surface points, which are coupled to structure, to the pressure at any 
interior fluid point "GiDi" is only of formal interest so that the results of Run 4 can be 
shown. Description of the output and its interpretation is given in section 6.7.4. 

The following figure shows the FE-model of the piston and the coupling surface of the 
fluid with numbered grid points. The numbers of the fluid coupling points (the so-called 
"participation points" PGP i (i = 1,2,..., 25), see section 6.7.4 can be found in table 10. 
The "participation factors" of these grid points for frequencies v = 128 Hz and v = 147 
Hz are listed in tables 10 and 1 1, respectively. 
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Figure 22. Coupling surface of the fluid and piston. 
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Table 10. Participation factors for point 553 at 128Hz 


COMPLEX DISPLACEMENT VECTOR 
(REAL/ IMAGINARY) 


POINT ID. 

TYPE 

T1 

T2 

T3 

R1 

1 

G 

-3.434522E-06 
-1 . 172005E-04 

1 . 172508E-04 
2 . 710505E-20 

1.562499E-02 
3 . 612052E-18 

2 . 499992E-01 
0.0 

2 

G 

-6.869054E-06 
-2 . 344014E-04 

2 . 345020E-04 
0.0 

3 . 125002E-02 
0.0 

4 . 999991E-01 
0.0 

3 

G 

-6 . 869051E-06 
-2 . 344013E-04 

2 . 345019E-04 
5 . 421011E-20 

3 . 125001E-02 
7.224104E-18 

4 . 999989E-01 
0.0 

4 

G 

-6 . 869055E-06 
-2 . 344015E-04 

2 . 345021E-04 
2 . 710505E-20 

3 . 125003E-02 
3 . 612052E-18 

4 . 999992E-01 
0.0 

5 

G 

-3 . 434528E-06 
-1 . 172007E-04 

1 . 172510E-04 
1 . 355253E-20 

1 . 562501E-02 
1.806026E-18 

2 . 499996E-01 
0.0 

6 

G 

-6 . 869044E-06 
-2 . 344011E-04 

2 . 345017E-04 
2.710505E-20 

3.124998E-02 
3. 612052E-18 

4.999984E-01 

0.0 

7 

G 

-1.373808E-05 
-4 . 688019E-04 

4 . 690031E-04 
5.421011E-20 

6.249992E-02 
7 . 224104E-18 

9 . 999963E-01 
0.0 

8 

G 

-1 . 373806E-05 
-4 . 688011E-04 

4 . 690024E-04 
0.0 

6.249982E-02 

0.0 

9 . 999947E-01 
0.0 

9 

G 

-1.373808E-05 

-4.688019E-04 

4 . 690031E-04 
5 . 421011E-20 

6 . 249991E-02 
7 . 224104E-18 

9 . 999962E-01 
0.0 

10 

G 

-6 . 869044E-06 
-2 . 344011E-04 

2 . 345017E-04 
5 . 421011E-20 

3 . 124997E-02 
7 . 224104E-18 

4 . 999984E-01 
0.0 

11 

G 

-6 . 869044E-06 
-2 . 344011E-04 

2.345017E-04 
2 . 710505E-20 

3 . 124998E-02 
3 . 612052E-18 

4 . 999984E-01 
0.0 

12 

G 

-1.373805E-05 
-4 . 688009E-04 

4 . 690022E-04 
0.0 

6 . 249979E-02 
0.0 

9 . 999942E-01 
0.0 

13 

G 

-1.373803E-05 
-4 . 688003E-04 

4 . 690015E-04 
5 . 421011E-20 

6.249970E-02- 

7.224104E-18 

9 . 999928E-01 
0.0 

14 

G 

-1 . 373805E-05 
-4.688007E-04 

4 . 690020E-04 
5 . 421011E-20 

6.249976E-02 
7 . 224104E-18 

9 . 999938E-01 
0.0 

15 

G 

-6 . 869036E-06 
-2.344008E-04 

2 . 345014E-04 
5 . 421011E-20 

3 . 124994E-02 
7 . 224104E-18 

4 . 999978E-01 
0.0 

16 

G 

-6 . 869063E-06 
-2.344017E-04 

2 . 345023E-04 
5 . 42101IE-20 

3 . 125006E-02 
7 . 224104E-18 

4 . 999998E-01 
0.0 

17 

G 

-1.373812E-05 
-4 . 688032E-04 

4 . 690045E-04 
5 . 421011E-20 

6.250009E-02 
7 . 224104E-18 

9 . 999991E-01 
0.0 

18 

G 

-1.373813E-05 
-4 . 688036E-04 

4 . 690049E-04 ' 
0.0 

6.250015E-02 

0.0 

1 . 000000E+00 
0.0 

19 

G 

-1.373812E-05 
-4 . 688031E-04 

4.690043E-04 
5. 421011E-20 

6.250007E-02 
7 . 224104E-18 

9. 999988E-01 
0.0 

20 

G 

-6.869055E-06 
-2 . 344014E-04 

2 . 345021E-04 
0.0 

3.125002E-02 

0.0 

4 . 999992E-01 
0.0 

21 

G 

-3 . 434535E-06 
-1.172010E-04 

1 . 172513E-04 
1 . 355253E-20 

1.562505E-02 

1.806026E-18 

2.500001E-01 

0.0 

22 

G 

-6 . 869094E-06 
-2 . 344028E-04 

2 . 345034E-04 
2 . 710505E-20 

3.125020E-02 
3 . 612052E-18 

5 . 000020E-01 
0.0 

23 

G 

-6 . 869106E-06 
-2.344032E-04 

2 . 345038E-04 
2 . 710505E-20 

3 . 125026E-02 
3 . 612052E-18 

5 . 000029E-01 
0.0 

24 

G 

-6 . 869093E-06 
-2 . 344028E-04 

2 . 345034E-04 
0.0 

3 . 125020E-02 
0.0 

5 . 000020E-01 
0.0 

25 

G 

“3 . 434535E-06 
“1 . 172010E-04 

1 . 172513E-04 
1.355253E-20 

1 . 562505E-02 
1 . 806026E-18 

2 . 500001E-01 
0.0 


- 83 - 


- 0.000028 
- 0.000033 
— 0 . 000038 
- 0.000043 
- 0.000048 
- 0.000053 
— 0 . 000058 
- 0.000063 
- 0 . 000067 
- 0.000072 
- 0.000077 
- 0 . 000082 
- 0.000087 
- 0.000092 
- 0.000097 
- 0.000102 
- 0.000107 
- 0.000112 


17 

16 

15 

14 

13 

12 



3 

2 

1 



Figure 23. Grid point participation for fluid point 553 ( x = 0.15) at 128 Hz (Re(T3) from 
table 10). 
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Table 11. Participation factors for point 553 at 147 Hz 


COMPLEX DISPLACEMENT VECTOR 
(REAL/ IMAGINARY) 


POINT ID. 

TYPE 

Tl 

T2 

T3 

R1 

1 

G 

-1 . 901086E-06 
1 . 260128E-04 

1 . 260271E-04 
1 . 423015E-19 

1 . 562499E-02 
1 . 764271E-17 

2 . 499992E-01 
0.0 

2 

G 

-3 . 802178E-06 
2 . 520259E-04 

2 . 520546E-04 
2.710505E-19 

3 . 125002E-02 
3 . 360516E-17 

4 . 999991E-01 
0.0 

3 

G 

-3 . 802177E-06 
2 . 520258E-04 

2 . 520545E-04 
2 . 710505E-19 

3.125001E-02 
3 . 360516E-17 

4 . 999989E-01 
0.0 

4 

G 

-3.802179E-06 
2 . 520259E-04 

2 . 520546E-04 
2 . 574980E-19 

3 . 125003E-02 
3 . 192490E-17 

4 . 999992E-01 
0.0 

5 

G 

-1.901090E-06 

1.260130E-04 

1 . 260273E-04 
1.287490E-19 

1.562501E-02 
1 . 596245E-17 

2 . 499996E-01 
0.0 

6 

G 

-3 . 802173E-06 
2 . 520255E-04 

2 . 520542E-04 
2 . 846031E-19 

3 . 124997E-02 
3 . 528542E-17 

4 . 999984E-01 
0.0 

7 

G 

-7 . 604342E-06 
5 . 040508E-04 

5 . 041081E-04 
5 . 421011E-19 

6.249991E-02 
6 . 721032E-17 

9 . 999963E-01 
0.0 

8 

G 

-7 . 604330E-06 
5 . 040500E-04 

5 . 041073E-04 
5 . 692061E-19 

6 . 249981E-02 
7 . 057084E-17 

9 . 999947E-01 
0.0 

9 

G 

-7 . 604342E-06 
5.040508E-04 

5 . 041081E-04 
5.421011E-19 

6 . 249991E-02 
6. 721032E-17 

9 . 999962E-01 
0.0 

10 

G 

-3.802173E-06 
2 . 520255E-04 

2 . 520542E-04 
2 . 710505E-19 

3 . 124997E-02 
3 . 360516E-17 

4 . 999984E-01 
0.0 

11 

G 

-3.802173E-06 
2 . 520255E-04 

2 . 520542E-04 
2.439455E-19 

3 . 124997E-02 
3.024465E-17 

4 . 999984E-01 
0.0 

12 

G 

-7 . 604326E-06 
5. 040497E-04 

5 . 041071E-04 
5 . 421011E-19 

6.249979E-02 
6 . 721032E-17 

9.999942E-01 

0.0 

13 

G 

-7 . 604316E-06 
5 . 040490E-04 

5 . 041064E-04 
5 . 692061E-19 

6.249970E-02 
7 . 057084E-17 

9 . 999928E-01 
0.0 

14 

G 

-7 . 604323E-06 
5.040496E-04 

5 . 041069E-04 
5 . 963112E-19 

6 . 249976E-02 
7 . 393136E-17 

9 . 999938E-01 
0.0 

15 

G 

-3 . 802169E-06 
2 . 520252E-04 

2 . 520539E-04 
2 . 981556E-19 

3 . 124994E-02 
3 . 696568E-17 

4 . 999978E-01 
0.0 

16 

G 

-3.802183E-06 
2 . 520262E-04 

2 . 520549E-04 
2 . 574980E-19 

3 . 125006E-02 
3 . 192490E-17 

4 . 999997E-01 
0.0 

17 

G 

-7 . 604363E-06 
5 . 040522E-04 

5 . 041096E-04 
5 . 421011E-19 

6 . 250009E-02 
6 . 721032E-17 

9 . 999991E-01 
0.0 

18 

G 

-7 . 604370E-06 
5 . 040527E-04 

5 . 041100E-04 
5.421011E-19 

6.250015E-02 
6 . 721032E-17 

1.000000E+00 

0.0 

19 

G 

-7 . 604361E-06 
5.040520E-04 

5 . 041094E-04 
5 . 692061E-19 

6.250007E-02 
7 . 057084E-17 

9 . 999988E-01 
0.0 

20 

G 

-3 . 802179E-06 
2 . 520259E-04 

2 . 520546E-04 
2 . 710505E-19 

3.125002E-02 
3 . 360516E-17 

4 . 999992E-01 
0.0 

21 

G 

-1.901094E-06 
1 . 260132E-04 

1 . 260276E-04 
1 . 355253E-19 

1 . 562505E-02 
1 . 680258E-17 

2 . 500001E-01 
0.0 

22 

G 

-3 . 802200E-06 
2 . 520273E-04 

2 . 520560E-04 
2.846031E-19 

3 . 125020E-02 
3.528542E-17 

5.000020E-01 

0.0 

23 

G 

-3 . 802207E-06 
2.520278E-04 

2 . 520565E-04 
2 . 710505E-19 

3 . 125026E-02 
3 . 360516E-17 

5 . 000029E-01 
0.0 

24 

G 

-3 . 802200E-06 
2 . 520273E-04 

2 . 520560E-04 
2 . 846031E-19 

3 . 125020E-02 
3 . 528542E-17 

5 . 000020E-01 
0.0 

25 

G 

-1.901094E-06 
1 . 260132E-04 

1.260276E-04 

1.423015E-19 

1 . 562504E-02 
1 . 764271E-17 

2 . 500001E-01 
0.0 
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-0.000028 
-0.000033 
-0.000038 
-0 . 000043 
-0.000048 
-0.000053 
-0.000058 
-0 . 000063 
-0.000068 
-0.000073 
-0.000078 
-0.000083 
-0.000088 
-0 . 000093 
-0.000098 
-0.000103 
-0.000108 
-0.000113 


Figure 24. Grid point participation for fluid point 553 ( x = 0.15) at 147 Hz (Re(T3) from 
table 11). 

It is very clear that grid points, PGP/ (/ = 7,8,9,12,13,14,17,18,19) (see figure 22), have 
the most influence on the pressure level at fluid interior point 553, (x = 0.15). This is also 
shown by table 10 with figure 23 and table 1 1 with figure 24. 

7. 1.3.4 Transient response 

This analysis is not performed here. It can be easily done considering the remarks in sec- 
tion 6.6. 



7.2 Two-dimensional fluid-structure system 

In section 7.1 a one-dimensional, even though it was solved as three-dimensional, example 
problem was described in detail. In this section a simple two-dimensional system for 
which a solution exists will be used to illustrate all the capability in this development. De- 
tails of the theoretical solution can be found in reference 31. Although this problem is two- 
dimensional, all the capability is directly applicable to any three-dimensional system. 
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The unique aspect of this problem is the use of water as the fluid medium. In this particular 
case one can expect to see a greater effect of the fluid on the structure than in the previous 
example where the fluid medium was air. Through the use of these two sample problems 
one can observe that FSI is dependent on the type of the two interacting mediums. In prac- 
tice, the structure will always have influence on the fluid, whereas the fluid may or may not 
have a significant effect on the structure. This effect depends, to some extent, on the fluid 
density, or the relative weight of the fluid and the structure. 

7.2.1 Problem Description 

A simply supported beam has one side in contact with a finite acoustic fluid, water, as 
shown in figure 20. A unit width of the fluid is considered. 


P = 0 







§ . 

1 

3.048 m 


(10 ft) 


^ Bem 


'SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSj 



6.096 m 
(20 ft) 


1 


Figure 25. Beam with finite acoustic fluid 
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The physical properties of the beam and fluid used are listed in the following table. 


Beam I 

Weight Density (p s ) 

7.682E4 N/m 3 

(0.283 lb/in 3 ) 

Young's Modulus (E) 

2.068E11 Pa 

(30E6 psi) 

Moment of Inertia (I) 

0.1675E-4 m 4 

(40.24 12in 4 ) 

Length (L) 

3.048 m 

(10 ft) 

Fluid 1 

Weight Density (pf) 

9802.25N/m 3 

(62.4 lb/ft 3 ) 

Acoustic Speed (c f ) 

1524 m/sec 

(5000 ft/sec) 

Height of the Fluid (H) 

6.096 m 

(20 ft) 


7.2.2 Theoretical solution 

The coupled modal frequencies of the fluid-beam system are to be determined. The equa- 
tion of motion of the beam is 


m^+ El^ = q(x,t) 

dt 2 dx* 

where m is the mass of the beam per unit length, E is Young's modulus of the beam, / is 
the section moment of inertia, and q(x,t) is the distributed load on the beam. 

The governing acoustic field equation in Cartesian coordinates can be written as 


d 2 p. + ^p + d 2 p _ \ Ip 


dx 2 dy dz 2 


c) dt 2 


(143) 


The boundary conditions for the beam and fluid are 

y = o ) 


Beam: 

^ = 0 

at x = 0 and x-L 

(144) 


dx 2 \ 



Fluid: 

p = 0 

atx = 0,Landy = H 

(145) 


The fluid boundary condition implies that the fluid has a free surface on three sides. Of 
course, this is an impossible real condition, but it presents a solvable mathematical prob- 
lem. 
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At the Beam-Fluid interface: 


(146) 


?P = _ n £> 

dy Pf dt 2 


Solving equations (143) thru (146), the following frequency equations can be obtained: 


, 2 Q) 2 (nn'f n J , p f vmh(a n H)} nt ( nn\ 

for “- = t' I tJ >0> “T — } uJ 1 ' 


(147) 


for <0; let then 


o)Am + 


n = 1,2,3.' 


The mode shapes of the acoustic fluid pressure are expressed by 


for 


for 


at > 0, and let a 2 = -a 2 


(148) 


2 

a 2 > 0, q> (I (jc,y) = ^^ L ^sin^^j(sina„y-tana n //cosa„y) (149) 


®*{ x ’y) = Pf(I)n ^sin^^-j(sinh a n y - tanh a n H cosh a n y ) (150) 


The frequency equations (147) and (148) are numerically solved for a) n and the acoustic 
fluid pressure modes are calculated with the given set of input data. 

7.2.3 Two-dimensional Nastran finite element model 

The two-dimensional finite element model for FSI analysis consists of the beam and the 
fluid. Owing to symmetry about the center of the beam and fluid, only half of the beam 
and the fluid are modeled, using 16 BAR elements and 128 QUAD8 elements, respectively. 
The fluid-beam finite element model is shown in figure 26. MSC/NASTRAN MSGMESH 
program was used to generate the finite element models and the listing of the input data 
deck is shown in Appendix Cl. 
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Excitation Location 



Beam 

Figure 26. Beam/fluid - FE mesh. 

To illustrate the use of superelements, the beam was modeled with one superelement and a 
residual structure. The grid point on the plane of symmetry is in the residual structure and 
all other points are in superelement 1. 

7.2.4 Normal Modes Analysis - SOL 63, 70 

To understand the basic dynamic characteristics of the system and to compare with theoreti- 
cal results, the normal mode analyses for the following cases were performed: 

• normal modes of the structure only 

• normal modes of the fluid only 

• normal modes of the combined system 

7.2.4. 1 Beam-only results 

Five modes were computed below 5000 Hz. The first four natural frequencies and mode 
shapes are shown in figure 27. These are the typical mode shapes one would expect from a 
pinned-pinned beam. 
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NATURAL FREQUENC 
BERM 
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1ST MODE <50.0 HZ > 

2ND MODE <450.0 HZ) 



3RD KOOE <1250.0 HZ) 

4TH MODE <2450.0 HZ) 


Figure 27. Mode shapes of the beam only. 

7. 2.4. 2 Fluid with rigid boundary results 

Thirteen acoustic pressure modes were computed below 1000 Hz. The first four modes are 
displayed in figure 28, plotting the pressure variation along the x-axis. The first mode 
along the line of symmetry is identified as a vertical acoustic quarter wave mode, the sec- 
ond mode as a 3/4 wave mode, the third as 1-1/4 wave mode, etc. 


C-Sb 











Figure 28. Mode shapes of the fluid only. 

7. 2.4. 3 Coupled fluid-structure modal results 

This analysis was performed using SOL 70. The input deck is shown in Appendix C2. 

Five dry structure modes and 13 fluid modes were included in the analysis of the coupled 
system. The lowest modes of the acoustic pressure are shown in figure 29. The second 
mode at 270.14 Hz is slightly longer than a quarter wave. The third, Fifth, and sixth modes 
are close to 3/4, 1-1/4, and 1-3/4 waves, respectively. The results also indicate that the 
structural motion shifts the dominant acoustic frequencies. The finite element results are 
compared with the theoretical results based on the equations (147) and (148). Excellent 
agreement is obtained. The comparison of the natural frequencies is shown in table 12. 
Note that the fundamental beam frequency of 50 Hz has been shifted to 40.3 Hz in the 
coupled solution. This illustrates a significant effect of the fluid on the structural behavior. 
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Table 12: Comparison of coupled natural frequencies between 
theoretical solution and finite element results 


Structure only 


5.0000E+01 


4.5000E+02 


1.2500E+03 


2.4501E+03 


4.0504E+03 



Fluid only 

Coupled 

NASTRAN 

2.5692E+02 

4.0318E+01 

3.1151E+02 

2.7014E+02 

3.9880E+02 

3.2894E+02 

5.0188E+02 

4.1455E+02 

6.1270E+02 

4.2373E+02 

7.2758E+02 

5.1554E+02 

7.4855E+02 

6.2504E+02 

7.6877E+02 

7.3956E+02 

8.0763E+02 

7.5634E+02 

8.4474E+02 

7.8106E+02 

8.6249E+02 

8.2183E+02 

9.3032E+02 

8.5742E+02 

9.632 1E+02 

8.7754E+02 


Coupled exact 


3.9810E+01 


2.6997E+02 


3.2818E+02 


4.032 1E+02 


4.1401E+02 


5.1540E+02 


6.2519E+02 


7.3976E+02 


7.5689E+02 


7.7955E+02 
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7.2.5 Modal frequency response analysis - SOL 7 1 

It is difficult to obtain a frequency response theoretical solution for this problem. This cal- 
culation was performed only to further illustrate the Acoustic Procedure. Several different 
analyses were performed, each with the intent of demonstrating a different feature of the 
procedure. Two areas of absorption were defined on the interface of the beam and the 
fluid. Only the first area of absorption was considered in the following analyses. The first 
area started at x = 0 and continued to x = // 4. The other area started at x = //4 and ended at 
x = 1/2. The analyses that were performed and the respective inputs are listed below: 

• Without absorption, with residual flexibility effects. Appendix C3. 1 

• Uniform absorption, without residual flexibility effects. Appendix C3.2 

• Linearly varying absorption, without residual flexibility effects. Appendix C3.3 

• Linearly varying absorption, with residual flexibility effects, Appendix C3.4 

It should be noted that all these calculations were done as restarts from SOL 63 data base. 
Before each run, the Area Matrix Program was executed to create the appropriate input ma- 
trices. 

The load was applied to the structure and its location is shown in figure 26. Because of as- 
sumed symmetrical boundary condition, there is a symmetric load implied on the right-hand 
side. Also note that the load is applied to the interior point of superelement 1. Owing to 
this fact, one must make use of the LOADSET card in the Case Control Deck and LSEQ 
card in the Bulk Data Deck. For further input details refer to Appendix C. 

Some of the selected results are shown in figure 30 through 33. As can be expected, the 
peak pressure occurs at the first coupled mode (primarily a beam mode) and at the second 
coupled mode (primarily a fluid mode). From figures 30 and 31 it can be seen that absorp- 
tion does not affect the response of the first mode whereas the second peak is significantly 
reduced. 
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Figure 30. Fluid-beam system: no absorption damping - NORF=l. 



Figure 31. Fluid-beam system: uniform absorption damping - NORF=-l. 

Since the calculations represented by figures 32 and 33 have the same damping, without 
and with residual flexibility effects, one can evaluate the effect of residual flexibility. In 
this case the effect is quite small and is only noticeable near the second peak. Unfortu- 
nately, it is very difficult to predict whether residual flexibility will have an important role. 
Since it is not expensive to calculate, it should be included in all cases. 
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7.2.6 Fluid-beam grid participation factor calculation 

To complete the last step in the analysis sequence, a sample grid point participation calcula- 
tion was carried out for this problem. Following the instructions in section 6, the input 
data were prepared which essentially consisted only of selecting the fluid grid point for 
which the structural grid point participation was required. In this case, point 21216 was 
selected. Also the frequencies for which the calculation is to be made have to be selected. 
Only one frequency, 40 Hz, was chosen which corresponds to the first peak in the fre- 
quency response curves. The listing of the input deck can be found in Appendix C4. table 
13 presents a selected output of points. As expected, the maximum contribution comes 
from the point closest to the middle of the beam since the first mode has a peak motion at 
the center. 

Table 13 Participation factors for point 21216 at 40 Hz 


COMPLEX DISPLACEMENT VECTOR 
(REAL/IMAGINARY) 


POINT ID. 

TYPE 

T1 

T2 

T3 

Rl 

20001 

G 

-1.109855E-02 

3 . 724580E-02 

2 . 018416E-03 

1.335892E-02 



3 . 555394E-02 

1 . 053499E-04 

5 . 709099E-06 

0.0 

20003 

G 

-9.526432E-02 

3 . 198780E-01 

1 . 733476E-02 

1 . 147304E-01 



3.053643E-01 

8 . 489896E-04 

4 . 600826E-05 

0.0 

20005 

G 

-2 . 410586E-01 

8 . 103863E-01 

4 . 391628E-02 

2 . 906606E-01 



7.737055E-01 

1 . 851216E-03 

1.003207E-04 

0.0 

20007 

G 

-4 . 124141E-01 

1.389196E+00 

7 . 528301E-02 

4 . 982618E-01 



1 . 326569E+00 

2 . 317275E-03 

1 . 255773E-04 

0.0 

20009 

G 

-5.732617E-01 

1 . 936381E+00 

1.049359E-01 

6. 945202E-01 



1 . 849580E+00 

1. 55894 9E-03 

8 . 448224E-05 

0.0 

20011 

G 

-6 . 990637E-01 

2 . 368622E+00 

1 . 283599E-01 

8. 495518E-01 



2 . 263113E+00 

-3 . 561608E-04 

-1.930099E-05 

0.0 

20013 

G 

-7.806732E-01 

2 . 651931E+00 

1.437129E-01 

9.511660E-01 



2 . 534422E+00 

-2.497307E-03 

-1 . 353335E-04 

0.0 

20015 

G 

-8 . 196042E-01 

2 . 788084E+00 

1.510913E-01 

1 . OOOOOOE+OO 



2 . 664897E+00 

-3.828808E-03 

-2.074899E-04 

0.0 
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8. VEHICLE ACOUSTIC ANALYSIS 


In recent years, the designers of airplane and automobile structures have been increasingly 
concerned with the noise in passenger compartments. These concerns occur because of 
government regulations and competitive pressures. In modern automotive designs, for ex- 
ample, there is a greater possibility of "boom" and "harshness" noises because of the use 
of lightweight structures and elimination of fully isolated chassis frames. In the past these 
problems were attacked by experimental testing methods alone. This process is time-con- 
suming, expensive, and is performed only when the prototype models are available. It is 
also difficult to pinpoint the sources of the noise problems. 

The main purpose of the development of the acoustic procedure is to provide the design 
engineer with a noise analysis program that can be used in the preliminary design stage. It 
is unlikely that an absolute quantitative noise prediction can be obtained for a particular de- 
sign owing to the large number of uncertainties in the structure model and the attendant 
complications in finite element modeling of such structures. Nevertheless, using this pro- 
cedure, a relative comparison can be made between two competing designs. The number 
of prototype tests will be significantly reduced by using this approach. Combining this 
procedure with testing should lead to an effective and timely method for solving these diffi- 
cult problems. 

This chapter is devoted to an example analysis of a realistic automobile structure. The 
choice of this particular example was motivated by the second author who spent most of his 
professional career at Daimler Benz. Most of the steps in the acoustic procedure have al- 
ready been discussed in detail in previous sections. The main purpose in this section is to 
illustrate the practicality and simplicity of the procedure. One of the most important goals 
of this procedure is to present acoustic analysis as a simple extension of the typical struc- 
tural dynamic analysis where the fluid domain is treated as if it were just another part of the 
structure with different properties. Therefore, the acoustic analysis can be performed early 
in the design process by a structural dynamic engineer rather than being treated as a stand- 
alone acoustic problem. 


8.1 Structural Model 

The hidden line plot of the structural model is shown in figure 34. The internal structural 
details could not be shown because most of the picture would be black due to the fineness 
of the finite element model. This structural model consists of 26144 grid points. In this 
particular case the whole model is subdivided into 36 superelements and one residual 
structure. In the MSC/NASTRAN terminology, this model consisted of three superelement 
levels. There are 33 tip superelements, 3 collector superelements, and one residual struc- 
ture. For the purpose of normal modes calculation, the Guyan reduction method was used 
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to approximate the modal behavior of the grid points on the interior of each superelement. 
The number and the type of elements used in the structural model are shown in table 14. 



Figure 34. Vehicle structural model. 


Table 14 Structural model 

Element Type 

Number 

CBAR 

338 

CELAS1 

786 

CQUAD4 

19424 

CTRIA3 

8225 


The total number of structural grid points in the residual structure is 663. From the acous- 
tic analysis point of view the only additions to this model were the grid points and elements 
of the air in the passenger compartment. 

An example of a different and smaller structural model is shown in figure 35. This exam- 
ple illustrates a coarser modelling approach for the structure. This problem was solved 
without the use of superelements. 
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Figure 35. Coarse vehicle structural model. 


8.2 Fluid Model of the Passenger Compartment 

The model representing the air on the interior of the structure is shown in figure 36. This 
model consists of 973 fluid grid points, 394 CHEXA, and 494 CPENTA elements. Note 
that each fluid grid point has only one degree of freedom associated with each fluid grid 
point. It is easily seen that the fluid model is an extremely small addition to the total struc- 
tural model. Clearly, the fluid model is much coarser than the structural model. The mesh 
size was determined based on the frequency of the fluid that had to be accurately repre- 
sented. In this particular analysis, the maximum frequency of interest in the fluid model is 
only 200 Hz. In addition, it is known from previous experimental measurements that the 
structural model is capable of predicting reasonably accurate response only up to 100 Hz. 
Knowing the speed of sound in the air, and assuming a need for eight elements per wave- 
length, the resulting mesh is obtained. 

The only additional required input for the fluid model is the application of the unit pressure 
to the outside surface of the fluid. This is used by the Area Matrix Program to determine 
the coupling points between the structure model and the fluid model. If absorption were to 
be specified, additional unit pressure loads would have to be applied to the fluid on the sur- 
faces that are in contact with the absorbing materials. 

The fluid model corresponding to the structural model in figure 35 is shown in figure 37. 
The size of this model is approximately the same as the one in figure 36. This points out 
the fact that even though the structural models can be significantly different in size, the fluid 
mesh is very similar since the volume that is modelled is approximately the same. 


I NiiMilGNAUJt tOttl 
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Figure 36. Vehicle fluid model. 



Figure 37. Vehicle fluid model (corresponding to coarse str. model). 


8.3 Modal Analysis of the Complete Model, Structure and Fluid 

This analysis is a standard SOL 63 in MSC/NASTRAN. Only the key portions of the input 
deck are shown in Appendix Dl. The Bulk Data Deck could not be listed since it contains 
nearly 100,000 records. The only unusual aspect of this calculation is the use of two sub- 
cases for the residual structure in the Case Control Deck, otherwise it is the same as any 
other structural normal modes calculations. Note that the component modes were not 
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calculated for any superelement since there is no method card request for any superele- 
ments. Guyan reduction was used for each superelement and only the modes of the resid- 
ual structure were calculated. This is performed in two subcases. In the first, the modes of 
the residual for the structure model are calculated and in the second, the modes of the fluid 
model are calculated. A total of 200 modes were calculated for the structure with the first 
flexible mode at 7 Hz and the last mode at 150 Hz. Only 10 modes were extracted for the 
fluid with the highest mode at 200 Hz. For illustration purposes, figures 38 and 39 present 
several views of the first longitudinal fluid mode at approximately 70 Hz. 



Figure 38. Orthographic and side views of the first longitudinal fluid mode (70Hz). 



Figure 39. Top and front views of the first longitudinal fluid mode (70Hz). 


The modes of the coupled model were not calculated for this example, simply because they 
were not required for any further calculation and there was no special interest in their 
eigenvalues. 
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8.4 Response Analysis During Harmonic Excitation 

The loading conditions for a typical automobile structure are too numerous to list. Obvi- 
ously, the input from the road is one of the interesting loading conditions. Input from the 
engine at the mounting points to the frame is another condition that can cause significant 
structurally borne noise in the passenger compartment. 

In this section, a calculation of the acoustic response due to the drive shaft unbalance will 
be illustrated. Ultimately, the value of any analysis is directly related to its ability to corre- 
late with actual measurements. A large number of variables can affect the comparison. Be- 
fore any comparisons can be made for acoustical response, one must ensure that the struc- 
tural responses, such as natural frequencies, accelerations response, etc., correlate with 
measured data. 

For complex structures such as automobiles, that is a very difficult task. There is a signifi- 
cant existing effort, with a generic name of System Identification, whose goal is to assist 
the analyst in correlating analytical to experimental results. This is accomplished by adjust- 
ing the analytical model to better match the test data. Inherent in this approach is the as- 
sumption that the test results are correct. This assumption is not valid in many cases, but it 
is the best that can be done. 

A condensed input deck for this particular analysis is listed in Appendix D3.1. Figure 39 
presents the measured pressure versus the calculated pressure at the rear passenger's ear 
due to the drive shaft unbalance. 


Pressure (dB) 



Figure 40. Comparison of the test and analysis pressure response at the 
rear passenger's ear - drive shaft unbalance. 
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Considering the above remarks, the comparison is very good. There is a 5dB difference in 
the results for frequencies above 70 Hz, but this is primarily due to the diminishing ability 
of the structural model to match the measured natural frequencies above 70 Hz. 

Clearly data such as these are valuable to a design engineer in correcting any potential noise 
problems. As stated previously, even if the comparison were not very good, relative com- 
parisons between different designs are useful in selecting the best design. 


8.5 Grid Point Participation Factor Calculation 

One of the most important aspects of the acoustic procedure is the ability to calculate the 
structural grid point participation in the generation of noise at a particular point in the fluid 
domain. With this information the design engineer can readily pinpoint the source of the 
structurally borne noise. Once the location of the disturbance is known, there are a number 
of available options for reducing the noise level. These include structural modifications, 
addition of absorbing materials, active noise suppression and other techniques. Participa- 
tion factors change with excitation frequency and sometimes certain parts of the structure 
generate noise at one frequency and subtract from the total noise at another frequency. This 
leads to conflicting requirements that can make the ultimate solution difficult to find. 

The ear locations of the driver and passenger are the points at which the pressure values 
were calculated during the harmonic response. Pressure responses at the rear passenger's 
and the driver's ear are shown in figures 41 and 42, respectively. 


Pressure (dB) 



Figure 41. Pressure response at the rear passenger's ear. 
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Pressure (dB) 



Figure 42. Pressure response at the driver's ear. 

It can easily be seen from figure 41 that there is a significant pressure at 64 Hz for the rear 
passenger. Similarly, from figure 42, there is a pressure peak at 86 Hz for the driver. The 
peak near 100 Hz is ignored since the validity of the model is questionable in that range. 

Therefore, to find the source of the pressure, grid participation factor calculation is carried 
out at 64 Hz for the rear passenger and for the driver at 86 Hz. For purpose of clarity, the 
contours of the participation factors are shown on the outside surface of the fluid (see sec- 
tion 6.7.4). The results are seen in figures 43 and 44. From figure 43 one can identify the 
source of noise for the rear passenger to be primarily the rear window. Clearly, the source 
of noise for the driver is the vibration of the front window. 

For a typical design process, more plots of this kind would have to be generated for other 
peaks in the pressure response curves. Fortunately, it is very inexpensive in CPU time and 
man hours to generate large amounts of grid participation factor data. The more difficult 
and time-consuming task is to find a feasible solution from a number of inevitably conflict- 
ing results. 
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Figure 43. Grid participation factors for the rear passenger’s ear at 64 Hz. 
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Figure 44. Grid participation factors for the driver’s ear at 86 Hz. 
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APPENDIX A : Notations 


\ 

£ 5 

b 

C 

C b 


d B 

d r 

d, 

D f 

D, 

D, 


F. 

E-f 

E P 

Es 

fs 

Sf 

/.. 



.matrix represenring absorbing parts of the fluid boundary S 

fluid boundaries 

.body force acting on the fluid 

coupling matrix 

.submatrix of the coupling matrix due to boundary panel b 
.speed of sound 

.damping of the fluid on the boundary B 

.modal damping matrix of the fluid 

.modal damping matrix of the structure 

.damping matrix of the fluid 

.damping matrix of the structure 

.boundary absorption matrix (normal impedance Z ) 

.index indicating fluid quantities 

.total load on the structure 

.load onto the fluid due to sources inside the fluid 

.pressure-induced load onto the structure 

.external load onto the structure 

.modal transformed structural load 

.modal transformed fluid load (symmetric form) 

.imaginary unit 
.indices 

.unit matrix (fluid) 

.wave number = (o/c^j 

.modal stiffness matrix of the fluid 

.modal stiffness matrix of the structure 

.fluid stiffness matrix 

.structure stiffness matrix 

.differential operator 

.fluid load vector due to body forces 

.fluid load vector due to external time-dependent pressure 

.structure load vector due to fluid interaction 

.fluid load vector due to structure interaction 
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n- 

v 

M c 

M f 

Ms 

m f 

m,. 

n .. 

n„ • 


n, ., 

*/ 

Ns- 

ND 

A 

p.. 

P o- 

p.. 

Ps- 

P. . 

q .. 

Q . . 


5 .. 

r... 

r.. 

[f 




C/ 


V, 


structure load vector due to external time-dependent load 

fluid load vector due to added fluid mass 

Fluid- structure coupling mass matrix 

fluid mass matrix 

structure mass matrix 

.modal mass matrix of the fluid 

.modal mass matrix of the structure 

.unit normal vector (components: n„ n v n z ) 

.total number of boundary panels 

.total number of fluid modes 

.total number of structure modes 

.trial function in the fluid domain 

.trial function in the structure domain 

.number of nodes required to discretize a wavelength 

.fluid pressure, total value 

.fluid pressure, reference value 

.fluid pressure, ambient value 

.external time-dependent pressure 

.column vector of the nodal values of the pressure 

.added fluid mass per unit volume and time 

.added fluid mass per unit volume, Q = qd 

.vector ( x,y,z ), x,y,z spatial coordinates 

.index indicating structural quantities 

.fluid boundary 

.time variable 

.transformation matrix for symmetry transforms 
.transpose of matrix 

.fluid displacement field, u f = (« A , u fy , u /z j 
.structure displacement field, u s = («„, u^, 

.structure displacement component perpendicular to the fluid boundary 
(uj = u f n) 

.structure displacement vector 
.fluid velocity field 

.fluid velocity component perpendicular to the fluid boundary 
.volume of the fluid domain 
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Z 0 

Z„ 

z, 

Z, 

y t yOt 



22 

P 

y xy' y ix 

<t> 

<>/ 

<Ps 

v 

A 

P 

Pr •••• 

P 

C 

«/ - 

Sx 

0) 

"/ 

J 

J 

s 

1 

V 

V* 

V 2 


scalar test function 

wave resistance (for plane waves), Z 0 = p / c / 

acoustic impedance 

specific acoustic resistance 

.specific acoustic reactance 

.residual flexibility terms of the fluid 

.acoustic modal frequency response function 

.specific acoustic admittance (p = y + ia) 

.strain components in x direction 

.body force potential 

.fluid modal matrix 

.structure modal matrix 

.displacement potential field 

.fluid density, total value 

.fluid density, reference value 

.fluid density, ambient value 

.stress components in x direction 

.modal coupling matrix 

.fluid modal variable 

.structural modal variable 

.radian frequency (coupled system) 

.fluid radian frequency 

.volume integral 

.surface integral 


.partial surface integrals (/=1,. ..,5) 

f d d d 


.gradient operator, V = , — , — , 

ydx dy dz 

•gradient operator on the fluid boundary 


Laplace operator, V 2 = A 
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APPENDIX B : INPUT DECKS FOR SECTION 7.1 
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APPENDIX Bi 


Run 1 Input - One-Dimensional Piston/Tube Model 


★★★★★★★★★★★★★★ 
** R U N 1 ** 


NASTRAN SYSTEM (105) =2 

NASTRAN SYSTEM (7) =10, BUFFSIZE=4609 $ 

ID FLUID STRUCTURE 
SOL 63 
TIME 999 
DIAG 8,13,20 
$ 

$$$$$$$$$$$ 

♦READ CDMAP 

$$$$$$$$$$$ 

$ 

CEND 

TITLE= XMPLO 

SUBTITLE=SPRING-PISTON-TUBE SYSTEM / UNCOUPLED E I GENANAL Y SIS 
ECHO=SORT 
$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ 

$ USER INSTRUCTIONS : 

$ 

$ 

$ 1. MUST INCLUDE A SET I = J THRU K AND PARTN = I IN THE 

$ CASE CONTROL DECK WHERE, 

$ 

$ I SOME INTEGER VALUE 

$ J THRU K = ALL FLUID GRID POINTS AND SPOINTS 

$ 

$ 2. TWO METHOD CARDS MUST BE SPECIFIED. FOR THE STRUCTURE 

$ EIGENVALUE EXTRACTION THE METHOD CARD IS PLACED IN THE 

$ FIRST SUBCASE FOR THE RESIDUAL STRUCTURE. FOR THE FLUID 

$ THE METHOD CARD IS PLACED IN THE SECOND SUBCASE FOR THE 

$ RESIDUAL STRUCTURE. GENERALLY, THE TWO METHOD CARDS 

$ WILL POINT TO DIFFERENT EIGRL CARDS IN THE BULK DATA 

$ DECK SINCE THE STRUCTURE AND THE FLUID CAN HAVE QUITE 

$ DIFFERENT EIGENSPECTRUMS . 

$ 

$ 3. ONE MUST ALSO SPECIFY A LOADSET CARD IN THE CASE CONTROL 

$ DECK THAT POINTS TO THE LSEQ CARD IN THE BULK DATA DECK. 

$ THE FIRST LOAD ON THE LSEQ CARD MUST BE THE PRESSURE 

$ LOAD ON THE OUTER BOUNDARY OF THE FLUID. THE DIRECTION 

$ OF THE PRESSURE MUST BE INWARD. ADDITIONAL LOADS MAY BE 

$ SPECIFIED TO IDENTIFY THE AREAS THAT WILL BE TREATED BY 

$ THE ABSORPTION MATERIAL. EACH AREA MAY HAVE A DIFFERENT 

$ COMPLEX IMPEDANCE . 

$ THE LOADS FOR THE FREQUENCY RESPONSE ANALYSIS MUST BE 

$ SPECIFIED IN SOL 71 WHEN THE ANALYSIS IS RESTARTED. 

$ 

$ 4. ONE MUST REQUEST AT LEAST DISP(PLOT) = ALL FOR ALL THE 

$ SUPERELEMENTS THAT ARE IN CONTACT WITH THE FLUID 
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$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ 

SEALL=ALL 

$ 

DISP (PLOT) =ALL 

$ 

$ FLUID POINTS 
SET 66 = 1 THRU 1275 
LOADSET =101 
$ 

SPC=11 
SUBCASE 1 
SUPER 10 
METHOD=101 
$ 

SUBCASE 2 
PARTN = 66 
METHOD=201 
$ 

SUBCASE 3 
METHOD=202 
$ 

$ 

BEGIN BULK 
PARAM AUTOSPC YES 
PARAM COUPMASS1 
PARAM EPZERO l.E-6 

PARAM FACTOR 0 

PARAM GRDPNT 0 

PARAM DBNBLKS 100000 

$ 

$ 

$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
$$$$$$$$$$$$$$$$$$$$$$$$$$$ PART I $$$$$$$$$$$$$$$$$$$$$$$$$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$ STRUCTURE $$$$$$$$$$$$$$$$$$$$$$$$$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ 

$ *************** EIGENANALYSIS : SPRING-PISTON-SYSTEM ***************** 

$ 

SPOINT 100001 

SEQSET1 10 0 100001 

$ 

$ 

$ 1. EIGENVALUE EXTRACTION / STRUCTURE SE10 
EIGRL 101 1 

$ 

$ 2. EIGENVALUE EXTRACTION / STRUCTURE SE0 
EIGRL 201 2 

$ 

$ 

$ 

£ ************************ SUPPORT CONDITION ************************* 

$ (FREE DOF : TRANSLATION IN XI DIRECTION) 

$ 

SPC1 11 23456 10001 THRU 10025 
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23456 11111 


SPC1 11 
$ 

$ 

$ ********** SUPERELEMENT-ASSIGNMENT OF PISTON GRID POINTS *********** 
SESET 10 10001 THRU 10025 

$ 

$ 

$ 

5 ****************** THE SPRING - PISTON - SYSTEM ******************* 
$ ************** (MATERIAL DATA, GEOMETRY AND TOPOLOGY) *************** 
$ 

$ MATERIAL DATA / PISTON 


PSHELL 

10 

111 

1.0 

111 

111 

MAT1 

in 

2.06E5 


.314 l.E-15 


CONM2 

£ 

101802 

10013 

0 

0.01E-3 


$ SPRING (STIFFNESS IN 

XI -DIRECTION) 


CELAS2 

c* 

101801 

7.47474911111 

i 


9 

CBAR 

101803 

88 

11111 

10013 0. 

1 . 

PBAR 

88 

111 

l.E+2 

l.E+4 l.E+4 

2.E+4 


$ 

$ VISCOUS DAMPING 


CDAMP2 

$ 

GRID 

C 

101804 5.E-4 

11111 

1 


11111 

0. 

0 0.0 

0.0 

9 

$ GRID 
GRID 

POINTS / PISTON 
10001 

SURFACE 

0. 

0. 

0. 

GRID 

10002 

0. 

-6.25 

0. 

GRID 

10003 

0. 

-12.5 

0. 

GRID 

10004 

0. 

-18.75 

0. 

GRID 

10005 

0. 

-25. 

0. 

GRID 

10006 

0. 

0. 

6.25 

GRID 

10007 

0. 

-6.25 

6.25 

GRID 

10008 

0. 

-12.5 

6.25 

GRID 

10009 

0. 

-18.75 

6.25 

GRID 

10010 

0. 

-25. 

6.25 

GRID 

10011 

0. 

0. 

12.5 

GRID 

10012 

0. 

-6.25 

12.5 

GRID 

10013 

0. 

-12.5 

12.5 

GRID 

10014 

0. 

-18.75 

12.5 

GRID 

10015 

0. 

-25. 

12.5 

GRID 

10016 

0. 

0. 

18.75 

GRID 

10017 

0. 

-6.25 

18.75 

GRID 

10018 

0. 

-12.5 

18.75 

GRID 

10019 

0. 

-18.75 

18.75 

GRID 

10020 

0. 

-25. 

18.75 

GRID 

10021 

0. 

0. 

25. 

GRID 

10022 

0. 

-6.25 

25. 

GRID 

10023 

0. 

-12.5 

25. 

GRID 

10024 

0. 

-18.75 

25. 

GRID 

10025 

0. 

-25. 

25. 


$ 

$ RIGID BAR ELEMENTS / PISTON (RIGID PISTON SURFACE!) 

RBE2 14001 10013 1 10001 10002 10003 10004 

10005+XYZ1111 
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+XYZ1111 

10006 

10007 

10008 

10009 

10010 

10011 

10012 

10014+XYZ1112 







+XYZ1112 

10015 

10016 

10017 

10018 

10019 

10020 

10021 

10022+XYZ1113 







+XYZ1113 

<; 

10023 

10024 

10025 





9 

$ PLATE : 

ELEMENTS / 

PISTON 

(HERE ONLY FORMALLY USED! 

!) 


CQUAD4 

10001 

10 

10001 

10002 

10007 

10006 


CQUAD4 

10002 

10 

10002 

10003 

10008 

10007 


CQUAD4 

10003 

10 

10003 

10004 

10009 

10008 


CQUAD4 

10004 

10 

10004 

10005 

10010 

10009 


CQUAD4 

10005 

10 

10006 

10007 

10012 

10011 


CQUAD4 

10006 

10 

10007 

10008 

10013 

10012 


CQUAD4 

10007 

10 

10008 

10009 

10014 

10013 


CQUAD4 

10008 

10 

10009 

10010 

10015 

10014 


CQUAD4 

10009 

10 

10011 

10012 

10017 

10016 


CQUAD4 

10010 

10 

10012 

10013 

10018 

10017 


CQUAD4 

10011 

10 

10013 

10014 

10019 

10018 


CQUAD4 

10012 

10 

10014 

10015 

10020 

10019 


CQUAD4 

10013 

10 

10016 

10017 

10022 

10021 


CQUAD4 

10014 

10 

10017 

10018 

10023 

10022 


CQUAD4 

10015 

10 

10018 

10019 

10024 

10023 


CQUAD4 

10016 

10 

10019 

10020 

10025 

10024 



$ 

$ 

$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ PART II $$$$$$$$$$$$$$$$$$$$$$$$$$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ FLUID $$$$$$$$$$$$$$$$$$$$$$$$$$ 
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
$ 

$ ******************* EIGENANALYSIS : "TUBE" (SEO) ********************** 
EIGRL 202 6 

$ 

$ 

$ 

$ "PRESSURE ANALOGY METHOD" * 

$ 

$ COMPUTATION OF THE SOUND FIELD IN A CAVITY USING THE SIMILARITY BETWEEN 
$ THE (SCALAR) ACOUSTIC WAVE EQUATION AND A COMPONENT OF THE STRESS EQUA- 
$ TION, E.G. THE X-COMPONENT. 

$ THE FOLLOWING CONDITIONS HAVE TO BE FULFILLED: 

$ 

$ UX "=" P 

$ 

$ 1. SPC1-CARD: UY = UZ = 0 

$ 

$ 2. MAT 9 - INPUT: 

$ 

$ COEFFICIENTS OF THE ELASTICITY MATRIX "G" : 

$ 

$ G14 = G16 = G46 = 0 

$ 

$ Gil = G44 = G66 = 1/RHOF 

$ 

$ DENSITY RHO (ON MAT 9 CARD) : 

$ 
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RHO = 1/K 


= RHOF * ' CF * CF) 


MAT 9 -CARD : 


11 

Gil 

.0 

.0 

.0 

.0 

.0 

.0 +L3 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 +L4 

G44 

.0 

.0 

.0 

.0 

G 66 

1/K 



MATERIAL DATA: - DENSITY OF AIR: RHOF = 1.205 EXP (-12) T/ (MM*MM*MM) 

- SPEED OF SOUND: CF = 3.44 EXP (5) MM/S 

NOTE: THESE STRANGE UNITS HAVE TO BE USED IF THE UNIT OF LENGTH IS "MM" 
(MILLIMETER) INSTEAD OF "M" (METER) . 


5 ***************** FORMAL FLUID SUPPORT CONDITIONS ******************** 

$ 

$ DISPLACEMENT UX "=" PRESSURE P 
SPC1 11 23456 1 THRU 1275 

$ 

$ 

$ **************************** AIR-FILLED TUBE ************************* 

$ 

$ MATERIAL DATA 
$ "AIR" 

PSOLID 1 11 0 2 

$ 

$ "UX = P" 

MAT 9 11 829.88+9.0 .0 .0 .0 .0 .0 

+L3 .0 .0 .0 .0 .0 .0 .0 .0 

+L4 829.88+9.0 .0 .0 .0 829.88+97.012909 


$ GRID POINTS / FLUID 
GRID 1 


GRID 1275 1250. -25. 

$ 

$ VOLUME ELEMENTS / TUBE (FLUID) 
CHEXA 111 
+ 277 276 


CHEXA 800 1 

+ 1275 1274 
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$ 

$ 

$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
$$$$$$$$$$$$ PART III $$$$$$$$$$$$$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
$$$$$$$$$$$$ DEFINITION OF COUPLING AND ABSORBING SURFACES $$$$$$$$$$$$$ 
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
$ 

$ FRONT SURFACE (FLUID) 


CQUAD4 


20001 

2 

1 

2 

7 

6 

CQUAD4 


20002 

2 

2 

3 

8 

7 

CQUAD4 


20003 

2 

3 

4 

9 

8 

CQUAD4 


20004 

2 

4 

5 

10 

9 

CQUAD4 


20005 

2 

6 

7 

12 

11 

CQUAD4 


20006 

2 

7 

8 

13 

12 

CQUAD4 


20007 

2 

8 

9 

14 

13 

CQUAD4 


20008 

2 

9 

10 

15 

14 

CQUAD4 


20009 

2 

11 

12 

17 

16 

CQUAD4 


20010 

2 

12 

13 

18 

17 

CQUAD4 


20011 

2 

13 

14 

19 

18 

CQUAD4 


20012 

2 

14 

15 

20 

19 

CQUAD4 


20013 

2 

16 

17 

22 

21 

CQUAD4 


20014 

2 

17 

18 

23 

22 

CQUAD4 


20015 

2 

18 

19 

24 

23 

CQUAD4 

c 


20016 

2 

19 

20 

25 

24 

y 

$ REAR 

SURFACE (FLUID) 






CQUAD4 


30001 

3 

275 

274 

524 

525 

CQUAD4 


30002 

3 

274 

273 

523 

524 

CQUAD4 


30003 

3 

273 

272 

522 

523 

CQUAD4 


30004 

3 

272 

271 

521 

522 

CQUAD4 


30005 

3 

525 

524 

774 

775 

CQUAD4 


30006 

3 

524 

523 

773 

774 

CQUAD4 


30007 

3 

523 

522 

772 

773 

CQUAD4 


30008 

3 

522 

521 

771 

772 

CQUAD4 


30009 

3 

775 

774 

1024 

1025 

CQUAD4 


30010 

3 

774 

773 

1023 

1024 

CQUAD4 


30011 

3 

773 

772 

1022 

1023 

CQUAD4 


30012 

3 

772 

771 

1021 

1022 

CQUAD4 


30013 

3 

1025 

1024 

1274 

1275 

CQUAD4 


30014 

3 

1024 

1023 

1273 

1274 

CQUAD4 


30015 

3 

1023 

1022 

1272 

1273 

CQUAD4 

$ 


30016 

3 

1022 

1021 

1271 

1272 

$ 

$ FORMAL FRONT WALL 






P SHELL 

2 

22 

1 

. E-5 22 


22 


MATl 

* 

22 

l.E-5 


1 . E- 

-5 



V 

$ FORMAL REAR WALL 






PSHELL 

3 

33 

1 

. E-5 33 


33 


MATl 

33 

l.E-5 


l.E- 

-5 




$ 

$ 

$ *************** DEFINING COUPLING AND ABSORBING SURFACES ************* 
$ ' 

$ FRONT SURFACE ("COUPLING" AND "ABSORBING" / "NON-ABSORBING" ) 
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PLOAD2 

200 

1. 20001 

THRU 

20016 

PLOAD2 

c 

201 

1. 20001 

THRU 

20016 

$ REAR 

SURFACE 

("ABSORBING") 



PLOAD2 

$ 

202 

1. 30001 

THRU 

30016 


$ NOTE: NO "PLOAD2 200" ==> RIGID SURFACE 
$ NO "PLOAD2 201,202" IF NO ABSORPTION DAMPING 

$ 


$ 

LOADSET 



$ 

ID 



LSEQ 

101 

151 

200 

$ 

LSEQ 

101 

301 

201 

LSEQ 

101 

302 

202 


$ 

$ 

$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ PART IV $$$$$$$$$$$$$$$$$$$$$$$$$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$$$$$$$$$$$$$$$$$$$$$$$$$ AREA MATRIX PROGRAM INPUT $$$$$$$$$$$$$$$$$$$$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 


$ 

$ EXPLANATION OF 
$ 

$ l.,0. 

$ 

$ 1 , 1 . 

$ 

$ N 
$ 

$ - 1,-1 
$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 1 
$ 

$ 

$ 

$ 

$ 

$ 

$ 200 . 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 


"AREA MATRIX INPUT" (SEE BELOW) : 

=> TOL1, TOL2 (REAL, REAL) 

=> IFF, FF (INTEGER, REAL) 

=> ANS = Y OR N (Y NOT USED) 

=> NR (1 ) , NR (2) 

NR(I) , 1=1,2 "=" NO. OF ABSRPTION REGIONS (INTEGERS) 

NR (I) <0 : NO ABSORPTION DATA WILL BE GENERATED 
FOR I'TH REGION 

NR(I) >0 : ABSORPTION DATA WILL BE GENERATED 
FOR I'TH REGION 


=> NF (INTEGER) 

SPECIFIES THE NUMBER OF FREQUENCIES AT WHICH THE AB- 
SORPTION DATA WILL BE PROVIDED. THIS IS REQUIRED IF 
ABSORPTION IS DEFINED AS A FUNCTION OF FREQUENCY. 

IF A CONSTANT VALUE OF ABSORPTION IS DESIRED, 

SET NF=1. 

=> F (I) , 1=1, NF (REALS) 

F (1) =200 . 

F (I) : FREQUENCIES AT WHICH THE ABSORPTION DATA WILL 
BE PROVIDED. 

NOTE: THE ABSORPTION DATA FOR ALL RETAINED REGIONS 
MUST BE SUPPLIED FOR ALL THE FREQUENCIES EVEN THOUGH 
SOME REGIONS MAY NOT HAVE FREQUENCY DEPENDENT AB- 
SORPTION PROPERTIES. 
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</></></></></></></></></>-</></></></></></></></></></></> 


4.145-7 => RHOC: WAVE RESISTANCE (FOR PLANE WAVES IN THE FLUID) 

RHOC = CF * RHOF 

CF: SPEED OF SOUND IN THE FLUID (HERE: AIR) 

CF=3 . 44*EXP (5) MM/S 
RHOF: DENSITY OF THE FLUID (HERE: AIR) 

RHOF=l . 205*EXP (-12) T/ (MM*MM*MM) 


2.28-7,-1.456-6 => CR(I),CI(I) (REALS) 

CR(I) : REAL PART OF ACOUSTIC ADMITTANCE 
Cl (I): IMAGINARY PART OF ACOUSTIC ADMITTANCE 

CR(I)=RHOC*ZR(I) / (ZR(I) *ZR(I) +ZI (I) *ZI (I) ) 

Cl (I) =-RHOC*ZI (I ) / (ZR (I ) *ZR (I ) +ZI(I)*ZI(I) ) 

ZR : REAL PART OF NORMAL IMPEDANCE 
ZI : IMAGINARY PART OF NORMAL IMPEDANCE 


ENDDATA 

/EOF 

1.. 0. 

1 . 1 . 

N 

- 1,-1 

1 

200 . 

4.145-7 

2.28-7,-1.456-6 


i 
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APPENDIX B2 


Run 2 Input - One-Dimensionat Piston/Tube Model 


** R U N 2 ** 


NASTRAN SYSTEM(105) =2, BUFFSIZE=4609 $ 

NASTRAN DBSET 1= (DB01, DB15) , DBSET 2= (DB01, DB15) $ 

NASTRAN DBSET 15=(DB15) $ 

ID FLUID STRUCTURE 
SOL 70 
TIME 999 
DIAG 8,13,20 
$ 

$$$$$$$$$$$ 

♦READ CDMAP 

$$$$$$$$$$$ 

$ 

CEND 

TITLE = XMPL02 

SUBTITLE = SPRING-PISTON-TUBE SYSTEM / COUPLED E I GENANALY SIS 
ECHO = NONE 

$ 

SET 99 = 0 
SEMG = 99 
SPC = 11 
$ 

$DISP=ALL 
DISP (PLOT ) =ALL 
$ 

$ SELECTION OF STRUCTURAL POINTS FOR PLOTTING 

$ 

$ SET 66 = 10001 THRU 10025 
$ DISP = 66 
$ 

SUPER = ALL 

$ 

SUBCASE 11 

LABEL = RESIDUAL STRUCTURE 
METHOD =11 

$ 

$ 

BEGIN BULK 
PARAM, LMODES, 7 
PARAM, RESDUAL, -1 
PARAM, DLOAD, -1 
$ 

$ REAL EIGENANALYSIS 

$ 

EIGR, 11,MGIV, , , , 7, , 1 . -9, +E1 
+E1,MAX 
$ 

$ COMPLEX EIGENANALYSIS 

$ 
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</>«/> <o *</></> -c/> 


+EC1 


EIGC, 12, HESS, MAX, , , , , , 
+ EC1, ,,,,,, 7 


READ DMI 


ENDDATA 
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APPENDIX B3.1 : Run 3.1 Input • One-Dimensional Piston/Tube Model 


** R U N 3.1 ** 

★ ★★★★★★★★★★★★★★A- 

NASTRAN SYSTEM (105) =2 

NASTRAN SYSTEM (7) =10, BUFFSIZE=4609 $ 

NASTRAN DBSET 1= (DB01, DB15) ,DBSET 2= (DB01, DB15) $ 

NASTRAN DBSET 15=(DB15) $ 

ID FLUID STRUCTURE 
SOL 71 
TIME 60 
DIAG 8,13,20 
$ 

$$$$$$$$$$$ 

*READ CDMAP 

$$$$$$$$$$$ 

$ 

CEND 

$ 

TITLE = XMPL03#1 

SUBTITLE = SPRING-PISTON-TUBE SYSTEM / FREQUENCY RESPONSE 
ECHO=SORT 
$ 

SET 99=0,10 
SELG =99 
SELR =99 
$ 

SPC = 11 

$ 

$ 

$ SELECTION OF STRUCTURAL AND FLUID POINTS FOR PLOTTING 

$ 

SET 60 = 10013 
SET 61 = 13,553 
$ 

DISP(SORT2, PHASE, PLOT) = 60 
VELO(SORT2, PHASE, PLOT) = 61 
$ 

$ 

LOADSET =10 

$ 

SUBCASE 1 
SUPER =10 
$ 

SUBCASE 2 
DLOAD = 1 
FREQ = 2 

SDAMPING = 3 

$ 

$ 

OUTPUT (XYPLOT) 

PLOTTER NAST 
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$ 

SEPLOT 0 
PENSIZE = 4 
XPAPER =10. 

YPAPER =10. 

XGRID = YES 

YGRID = YES 

UPPER TICS = -1 
RIGHT TICS = -1 
YVALUE PRINT SKIP = 1 
XDIVISIONS = 12 
YDIVISIONS = 10 
XTITLE = FREQUENCY (HZ) 

YTITLE = PRESSURE 

$ 

TCURVE = PRESSURE AT FLUID GRID POINT 13 
XYPLOT, XYPUNCH VELO / 13(T1RM) 

$ 

TCURVE = PRESSURE AT FLUID GRID POINT 553 
XYPLOT, XYPUNCH VELO / 553(T1RM) 

$ 

$ 

SEPLOT 10 
PENSIZE = 4 
XPAPER =10. 

YPAPER =10. 

XGRID = YES 

YGRID = YES 

UPPER TICS = -1 
RIGHT TICS = -1 
YVALUE PRINT SKIP ■= 1 
XDIVISIONS = 12 
YDIVISIONS =10 
XTITLE = FREQUENCY (HZ) 

YTITLE = AMPLITUDE IN MM 

$ 

TCURVE = PISTON DISPLACEMENT / GRID POINT 10013 
XYPLOT, XYPUNCH DISP/10013 (T1RM) 

$ 

$ 

BEGIN BULK 
PARAM, NORF, -1 
PARAM, SMALL, 0 . 

PARAM, LMODES, 7 
PARAM, DLOAD, -1 
$ PARAM, RESDUAL, -1 
$ 

$ EXCITATION AT PISTON CENTER 10013 

$ 

FORCE, 11, 10013, ,2.1885,1.0, , 

RLOAD1 1 12 13 

LSEQ, 10, 12, 11 

$ 

TABLED1 13 
+TAB 

+TAB 100. 1. 300. 1. ENDT 

$ 

$ MODAL DAMPING INPUT 
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$ 

TABDMP1 3 

+CF2 

+CF2 -1000. .0 

+CF3 

+CF3 ENDT 

$ 

$ FREQUENCY DOMAIN 

$ 

FREQ1 2 100 

$ 

$ 

$ 

♦READ DMI 


$ 

$ 

$ 

ENDDATA 


0 . 


0.1 


.0 


1000 


.1 


.0 1000 . .0 
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APPENDIX B3.2 : Run 3.2 Input - One-Dimensional Piston/Tube Model 


★ ★★★★★★★★★★★★★★■A 

** R U N 3.2 ** 


NASTRAN SYSTEM (105) =2 

NASTRAN SYSTEM (7) =10, BUFFSIZE=4609 $ 

NASTRAN DBSET 1= (DB01, DB15) , DBSET 2= (DB01, DB15) $ 

NASTRAN DBSET 15=(DB15) $ 

ID FLUID STRUCTURE 
SOL 72 
TIME 60 
DIAG 8,13,20 
$ 

$$$$$$$$$$$ 

*READ CDMAP 

$$$$$$$$$$$ 

$ 

CEND 

$ 

TITLE = XMPL03#2 

SUBTITLE = SPRING-PISTON-TUBE SYSTEM / TRANSIENT RESPONSE 
ECHO = SORT 
$ 

SET 99=0,10 
SELG = 99 
SELR =99 
$ 

SPC = 11 

$ 

$ 

$ SELECTION OF STRUCTURAL AND FLUID POINT FOR PLOTTING 

$ 

SET 60 = 10013 
SET 61 = 13,553 
$ 

DISP(SORT2, PHASE, PLOT) = 60 
VELO(SORT2, PHASE, PLOT) = 61 
$ 

$ 

SUBCASE 1 
SUPER =10 
$ 

SUBCASE 2 
DLOAD = 1 
TSTEP = 2 
$ 

$ 

OUTPUT (XYPLOT) 

PLOTTER NAST 

$ 

SEPLOT 0 
PENSIZE = 4 
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XPAPER =10. 

YPAPER =10. 

XGRID = YES 

YGRID = YES 

UPPER TICS = -1 
RIGHT TICS = -1 
YVALUE PRINT SKIP = 1 
XDIVISIONS =12 
YDIVISIONS =10 
$ 

XTITLE = FREQUENCY (HZ) 

YTITLE = PRESSURE 

$ 

TCURVE = PRESSURE AT FLUID GRID POINT 13 (X=0.0) 
XYPLOT, XYPUNCH VELO / 13(T1RM) 

$ 

TCURVE = PRESSURE AT FLUID GRID POINT 553 (X=150MM) 
XYPLOT, XYPUNCH VELO / 553(T1RM) 

$ 

$ 

SEPLOT 10 
PENSIZE = 4 
XPAPER =10. 

YPAPER =10. 

XGRID = YES 

YGRID = YES 

UPPER TICS = -1 
RIGHT TICS = -1 
YVALUE PRINT SKIP = 1 
XDIVISIONS =12 
YDIVISIONS =10 
XTITLE = FREQUENCY (HZ) 

YTITLE = AMPLITUDE IN MM 

$ 

TCURVE = PISTON DISPLACEMENT / GRID POINT 10013 
XYPLOT, XYPUNCH DISP/10013 (T1RM) 

$ 

$ 

BEGIN BULK 
PARAM, NORF, -1 
PARAM, SMALL, 0 . 

PARAM, LMODES, 7 
PARAM, DLOAD, -1 
$PARAM, RESDUAL, -1 
$ 

$ 

$ EXCITATION AT GRID POINT 11111 
DLOAD 1 1.0 1.0 11 

$ 

TLOAD1, 11, 111, , ,112 

DAREA 111 11111 1 1.E+0 

$ 

TABLED1, 112, ,,,,,,, +TABL1 

+TABL1, -1 . 0, 0 . 0, 0 . 0, 0 . 0, 0 . 005, 1 . 0, 0 . 0195, 1 . 0, +TABL2 
+TABL2, 0 . 02, 0 . 0, 1 . 0, 0 . 0, ENDT 
$ 

$ TIME DOMAIN 
TSTEP, 2, 2000, 1 . 0-4, 1 
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</> -c/> </> 


$ 

*READ DMI 


ENDDATA 



APPENDIX B4 : Run 4 Input - One-Dimensional Piston/Tube Model 


*♦ R U N 4 ♦* 

★★★★★★★★★★★*★★ 

NASTRAN SYSTEM (105) =2 
NASTRAN BUFFSIZE=4609 $ 

NASTRAN DBSET 1= (DB01, DB14, DB15) , DBSET 2= (DB01, DB14, DB15) $ 

NASTRAN DBSET 15= (DB14, DB15) $ 

ID TEST, PART 
TIME 30 
DIAG 8,13,20 
$ 

$$$$$$$$$$$ 

♦READ CDMAP 

$$$$$$$$$$$ 

$ 

CEND 

$ 

TITLE = XMPL04 

SUBTITLE = PRING-PISTON-TUBE SYSTEM / GRID POINT PARTICIPATION FACTOR 
ECHO = SORT 
$ 

$ 

$ 

$ INSERT SET 911 (ID'S OF FLUID GRID POINTS THAT AREA ON THE FLUID 
$ SURFACE OR ALL) 

♦READ SET 

$ 

$ 

$ 

DISP = 911 

$ 

SET 66 = 128.0,147.0 
OFREQ = 66 
$ 

BEGIN BULK 

$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 


$ $ 

$ USER INSTRUCTIONS: $ 

$ 5 

$ CASE CONTROL DECK - $ 

$ $ 

$ SET j = FLUID GRID ID'S THAT ARE ON THE SURFACE OR $ 

$ = ALL $ 

$ DISP = j $ 

$ $ 

$ SET k = fl, f2, f 3, . . . fn $ 

$ OFREQ = k $ 

$ WHERE f 1, f 2, f 3, . . . fn ARE THE FREQUENCIES AT WHICH $ 

$ THE GRID PARTICIPATION FACTOR IS TO BE CALCULATED $ 

$ $ 
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$ THE REST OF THE CASE CONTROL DECK MUST BE THE SAME AS THAT 
$ IN SOLUTION 71 RUN. I.E. THE SUBCASE STRUCTURE FOR THE 
$ RESIDUAL MUST BE THE SAME. 

$ 

$ BULK DATA DECK - 

$ 

$ DMIG, PDOF, 0, 9, 1, 0, , , 1 
$ DMIG, PDOF, l,l,,GIDi,CID, 1.0 
$ 

$ IF ONE WANTS ONLY THE REAL DISPLACEMENT OUTPUT SET THE 
$ PARAM, REAL, 0 
$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ 

PARAM, REAL, -1 

$ 

DMIG, PDOF, 0, 6, 1, 0 
DMIG, PDOF, 553, 1, , 553,1, 1 . 0 
$ 

ENDDATA 
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APPENDIX C : INPUT DECKS FOR SECTION 7.2 
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APPENDIX Cl 


Run 1 Input - Two-Dimensional Fluid/Beam Model 


NASTRAN REAL=0, HICORE=300000, PREFOPT=2 $ 

ID FLUID STRUCTURE 
SOL 63 
TIME 30 
DIAG 8,13 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ 

$ USER INSTRUCTIONS: 

$ 

$ 1 . MUST INCLUDE A set i = j thru k AND partn = i IN THE 

$ CASE CONTROL DECK WHERE, 

$ 

$ i = some integer value 

$ j thru k = all fluid grid points and spoints 

$ 

$ 2. TWO METHOD CARDS MUST BE SPECIFIED. FOR THE STRUCTURE 

$ EIGENVALUE EXTRACTION THE METHOD CARD IS PLACED IN THE 

$ FIRST SUBCASE FOR THE RESIDUAL STRUCTURE. FOR THE FLUID 

$ THE METHOD CARD IS PLACED IN THE SECOND SUBCASE FOR THE 

$ RESIDUAL STRUCTURE. GENERALLY, THE TWO METHOD CARDS 

$ WILL POINT TO DIFFERENT eigrl CARDS IN THE BULK DATA 

$ DECK SINCE THE STRUCTURE AND THE FLUID CAN HAVE QUITE 

$ DIFFERENT EIGENSPECTRUMS . 

$ 

$ 3. ONE MUST ALSO SPECIFY A loadset CARD IN THE CASE CONTROL 

$ DECK THAT POINTS TO THE lseq CARD IN THE BULK DATA DECK. 

$ THE FIRST LOAD ON THE lseq CARD MUST BE THE PRESSURE 

$ LOAD ON THE OUTER BOUNDARY OF THE FLUID. THE DIRECTION 

$ OF THE PRESSURE MUST BE INWARD. ADDITIONAL LOADS MAY BE 

$ SPECIFIED TO IDENTIFY THE AREAS THAT WILL BE TREATED BY 

$ THE ABSORPTION MATERIAL. EACH AREA MAY HAVE A DIFFERENT 

$ COMPLEX IMPEDANCE . 

$ THE LOADS FOR THE FREQUENCY RESPONSE ANALYSIS MUST BE 

$ SPECIFIED IN SOL 71 WHEN THE ANALYSIS IS RESTARTED. 

$ 

$ 4. ONE MUST REQUEST AT LEAST DISP (PLOT) = ALL FOR ALL THE 

$ SUPERELEMENTS THAT ARE IN CONTACT WITH THE FLUID 

$ 

$ PARAMETERS : 

$ 

$ PARAM, CSTMOUT, 11 

$ PARAM, GEOMOUT , 1 2 

$ 

$ THESE PARAMETERS SPECIFY OUTPUT2 FORTRAN FILE UNITS 

$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ 

$ INSERT DMAP PROGRAM HERE 

$ 

CEND 
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TITLE = FLUID-STRUCTURE INTERACTION 
SUBTITLE = NORMAL MODES ANALYSIS 
ECHO = NONE i 

SEALL = ALL 
SPC = 1 
LOADSET = 11 
$ 

$ DEFINE THE PARTITION VECTOR FOR THE FLUID 

$ 

SET 66 = 20000 THRU 30000 

$ 

DISP (PLOT) = ALL 
SUBCASE 1 
SUPER = 1 
METHOD = 1 
SUBCASE 11 

LABEL = STRUCTURE MODEL 
PARTN =66 
METHOD = 2 
SUBCASE 12 

LABEL = FLUID MODEL 
METHOD = 3 
BEGIN BULK 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ 

$ DEFINE THE STRUCTURE MODEL 

$ 

EGRID, 1 
EGRID,2, ,60. 

$ 

GRIDG, 1, , 1345, 16, -1, -2 
CBARG, 1, 1, 1, , , 1 . 

PBAR, 1,1,7.845,40.24 
MAT1, 1, 3 . +7, , .3,7.324-4 
$ ' 

SPC1, 1,12, 10000 
SPC1, 1, 6, 10016 
$ 

$ DEFINE THE EIGENVALUE EXTRACTION DATA FOR THE 
$ STRUCTURE 
$ 

EIGRL, 1, , , 8 

EIGR, 2,MGIV, , , , 5, , 1 . -8, +E1 
+E1,MAX 
$ 

$ DEFINE THE SUPERELEMENT FOR THE STRUCTURE 

$ 

SESET, 1, 10000, THRU, 10015 
SPOINT, 11001, THRU, 11020 
SEQSET1, 1, 0, 11001, THRU, 11020 
$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ 

$ DEFINE THE FLUID MODEL 

$ 

EGRID, 3, ,60., 240. 
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EGRID, 4, , 0 . , 240 . 

GRIDG, 2, , 23456, 16, -1, -2, -3, , +G2 
+G2, 32,-4 

CGEN,QUAD8, 101,1,2 
PSHELL 1 11 1. 

MAT2, 11, 1 . 060+4,0 . , 0. , 1 . 060+4, 0. , 1 . 060+4,2 . 96-6 


SPCG, 1, 2, 1, AD, CD 
PARAM, COUPMASS, 1 
$ 

$ DEFINE THE LOAD ON THE FLUID TO CALCULATE 
$ THE AREAS 
$ 

FORCE, 1, 20000, , 1 . 250, , 1 . 

FORCE, 1, 20001, , 5 . 000, , 1 . 

FORCE, 1, 20002, , 2 . 500, , 1 . 

FORCE, 1, 20003, ,5.000, , 1 . 

FORCE, 1,20004, ,2.500,, 1. 

FORCE, 1, 20005, , 5 . 000, , 1 . 

FORCE, 1, 20006, , 2 . 500, , 1 . 

FORCE, 1,20007,, 5. 000,, 1. 

FORCE, 1, 20008, , 2 . 500, , 1 . 

FORCE, 1, 20009, , 5 . 000, , 1 . 

FORCE, 1, 20010, , 2 . 500, , 1 . 

FORCE, 1, 20011, ,5.000, , 1 . 

FORCE, 1, 20012, , 2 . 500, , 1 . 

FORCE, 1,20013,, 5. 000,, 1. 

FORCE, 1, 20014, , 2 . 500, , 1 . 

FORCE, 1, 20015, ,5.000, , 1 . 

FORCE, 1, 20016, , 1 .250, , 1 . 

$ 

$ THE FOLLOWING TWO LOAD SETS DEFINE THE AREAS 
$ FOR ABSORPTION DAMPING 
$ 

FORCE, 2, 20000, , 1 . 250, , 1 . 

FORCE, 2, 20001, , 5 . 000, , 1 . 

FORCE, 2, 20002, , 2 . 500, , 1 . 

FORCE, 2, 20003, , 5 . 000, , 1 . 

FORCE, 2, 20004, , 2 . 500, , 1 . 

FORCE, 2, 20005, , 5 . 000, , 1 . 

FORCE, 2, 20006, , 2 . 500, , 1 . 

FORCE, 2, 20007, , 5 . 000, , 1 . 

FORCE, 2, 20008, , 2 . 500, , 1 . 

$ 

FORCE, 3, 20008, , 2 . 500, , 1 . 

FORCE, 3, 20009, , 5 . 000, , 1 . 

FORCE, 3, 20010, , 2 . 500, , 1 . 

FORCE, 3, 20011, , 5 . 000, , 1 . 

FORCE, 3, 20012, , 2 . 500, , 1 . 

FORCE, 3, 20013, , 5 . 000, , 1 . 

FORCE, 3, 20014, , 2 . 500, , 1 . 

FORCE, 3, 20015, , 5 . 000, , 1 . 

FORCE, 3, 20016, , 1 . 250, , 1 . 

$ 


LSEQ, 11, 101,1 
LSEQ, 11, 102, 2 
LSEQ, 11, 103, 3 
$ 

$ DEFINE THE EIGENVALUE EXTRACTION DATA 
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$ FOR THE FLUID 
$ 

EIGRL, 3, , , 13 
$ 

ENDDATA 

$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
$$$$$ AREA MATRIX PROGRAM - USER INPUT $$$$$ 
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
$ 

$ TOL1, TOL2 

1.. 0. 

$ IFF,FF 

1 . 1 . 

$ SEARCH REGION 
N 

$ SELECT ABSORPTION AREAS 

- 1,-1 

$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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APPENDIX C2 


Run 2 Input - Two-Dimensional Fluid/Beam Model 


NASTRAN DBSET 1= (DB01, DB15) , DBSET 2= (DB01, DB15) $ 

NASTRAN DBSET 15=(DB15) $ 

ID FLUID STRUCTURE 
SOL 70 
TIME 30 
DIAG 8,13 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 


$ $ 

$ USER INSTRUCTIONS: $ 

$ $ 

$ 1. THE USER MUST INSERT THE DMAP.DAT DATA GENERATED BY THE $ 

$ AREA MATRIX FORTRAN PROGRAM IN THE DMAP SEQUENCE. THE $ 

$ LOCATION IS SHOWN BELOW. $ 

$ $ 

$ 2 . ALSO THE USER MUST INSERT THE DMI . DAT DATA GENERATED BY $ 

$ THE AREA MATRIX FORTRAN PROGRAM IN THE BULK DATA DECK. $ 

$ THE LOCATION IS SHOW IN THE BULK DATA DECK. $ 

$ $ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
$ $ 

$ ***** BEGIN OF THE ALTER PACKAGE FOR FSI ANALYSIS ***** $ 

$ $ 


$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ 

$ INSERT DMAP.DAT HERE 

$ 

$ INSERT DMAP PROGRAM HERE 

$ 

CEND 

TITLE = FLUID-STRUCTURE INTERACTION 
SUBTITLE = COUPLED EIGENVALUE ANALYSIS 
SET 99 = 0 
SEMG = 99 
SPC = 1 

SET 66 = 10000 THRU 10016, 

20000 THRU 23216 
DISP (PLOT) = 66 
SUPER = ALL 
SUBCASE 11 

LABEL = RESIDUAL STRUCTURE 
METHOD = 31 
BEGIN BULK 
PARAM, DLOAD, -1 
PARAM, RESDUAL, -1 
PARAM, LMODES, 100 

$ 

EIGR, 31,MGIV, , , , 10, , 1 . -9, +E1 
+E1,MAX 
$ 
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$ INSERT DMI.DAT HERE 
$ 


DMI* 

DM10000 



2 

123 1 . 250000000E+00 

DMI 

DMI0000 0 

2 

1 

0 

2624 2624 

DMI* 

DMI0001 



2 

27 1 . 250000000E+00 

DMI* 

DMI0001 



8 

33 5. 000000000E+00 

DMI* 

DMI0001 



14 

39 2 . 500000000E+00 

DMI* 

DMI 00 01 



20 

45 5 . 000000000E+00 

DMI* 

DMI0001 



26 

51 2 . 500000000E+00 

DMI* 

DMI0001 



32 

57 5 . 000000 00 QE+00 

DMI* 

DMI0001 



38 

63 2 . 500000000E+00 

DMI* 

DMI0001 



44 

69 5. 000000000E+00 

DMI* 

DMI0001 



50 

75 2 . 500000000E+00 

DMI* 

DMI0001 



56 

81 5 . OOOOOOOOOE+OO 

DMI* 

DMI0001 



62 

87 2 . 500000000E+00 

DMI* 

DMI0001 



68 

93 5. OOOOOOOOOE+OO 

DMI* 

DMI0001 



74 

99 2 . 500000000E+00 

DMI* 

DMI0001 



80 

105 5. OOOOOOOOOE+OO 

DMI* 

DMI0001 



86 

111 2 . 500000000E+00 

DMI* 

DMI0001 



92 

117 5. OOOOOOOOOE+OO 

DMI 

DMI0001 0 

2 

1 

0 

2624 122 

DMI 

SUPER 0 

2 

1 

0 

2 1 

DMI* 

SUPER 



1 

1 0. OOOOOOOOOE+OO 


* l.OOOOOOOOOE+OO 

ENDDATA 
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APPENDIX C3.1 : Run 3.1 Input - Two-Dimensional Fluid/Beam Model 


$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
$$$$$ AREA MATRIX PROGRAM - USER INPUT $$$$$ 
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
$ 

$ TOL1 , TOL2 

l.fO. 

$ IFF,FF 

1 , 1 . 

$ SEARCH REGION 
N 

$ SELECT ABSORPTION AREAS 

- 1,-1 

$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 


NASTRAN DBSET 1= ( DBO 1 , DB1 5 ) , DBSET 2= (DB01, DB15) $ 

NASTRAN DBSET 15=(DB15) $ 

NASTRAN REAL=0,HICORE=300000 $ 

ID FLUID STRUCTURE 
SOL 71 
TIME 30 
DIAG 8,13 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 


$ $ 

$ USER INSTRUCTIONS: $ 

$ $ 

$ 1. THE USER MUST INSERT THE DMAP.DAT DATA GENERATED BY THE $ 

$ AREA MATRIX FORTRAN PROGRAM IN THE DMAP SEQUENCE. THE $ 

$ LOCATION IS SHOWN BELOW. $ 

$ $ 

$ 2 . ALSO THE USER MUST INSERT THE DMI . DAT DATA GENERATED BY $ 

$ THE AREA MATRIX FORTRAN PROGRAM IN THE BULK DATA DECK. $ 

$ THE LOCATION IS SHOW IN THE BULK DATA DECK. $ 

$ $ 

$ 3. IF ABSORPTION DAMPING DATA, DMIG.DAT, GENERATED BY THE $ 

$ AREA MATRIX PROGRAM EXIST THE USER MUST ALSO INCLUDE $ 

$ THIS DATA IN THE BULK DATA DECK. THE LOCATION IS SHOW $ 

$ IN THE BULK DATA DECK. $ 

$ $ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
$ $ 

$ ***** BEGIN OF THE ALTER PACKAGE FOR FSI ANALYSIS ***** $ 

$ $ 


$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ 

$ INSERT THE DMAP.DAT HERE 
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$ 

$ INSERT THE DMAP ALTER HERE 

$ 

CEND 

TITLE = COUPLED FLUID-STRUCTURE FREQUENCY RESPONSE 
SUBTITLE = 3% BEAM, 2% FLUID DAMPING 
LABEL = NO ABSORPTION - NORF=l 
$ 

SET 99 = 1,0 
SET 67 = 21216 
SET 65 = 0 
$ 

SELG = ALL 
SELR = ALL 
SEDR =65 
$ 

SPC = 1 
VELO = 67 
FREQ = 1 
SDAMPING = 1 
$ 

SUBCASE 1 
SUPER =99 
LOADSET =11 
DLOAD =21 

$ 

OUTPUT (XYPLOT) 

CSCALE 4.0 
XPAPER = 105. 

YPAPER = 80. 

XAXIS = YES 
YAXIS = YES 
XGRIDLINES = YES 
YGRIDLINES = YES 
XMIN = 0. 

XMAX = 400. 

XTITLE = FREQUENCY (HZ) 

YTITLE = PRESSURE (PSI) 

TCURVE = PRESSURE AT GRID 21216 
XYPLOT VELO / 21216 (T1RM) 

BEGIN BULK 

$ 

$ DEFINE PARAMETERS 

$ 

PARAM, NORF , 1 $ RESIDUAL FLEXIBILITY FOR FLUID 

PARAM, SMALL, 0 . 

PARAM, LMODES, 100 
PARAM, DLOAD, -1 
$ 

$ DEFINE THE DYNAMIC LOADS 

$ 

FREQ1, 1, 5. , 1 . , 395 

$ 

FORCE, 1, 10008, , 1000. , , 1 . 

LSEQ, 11,101,1 

$ 

RLOAD1, 21, 101, , ,1 

$ 


- 152 - 



TABLED1, 1, ,,,,,,, +TD1 
+TD1, 0 . , 1 . , 1000 . , 1 . , ENDT 
$ 

$ DEFINE MODAL DAMPING FOR FLUID AND BEAM 
$ 2% FOR FLUID AND 3% FOR THE BEAM 
$ NOTE: THE FLUID MODAL DAMPING IS DEFINED IN 
$ THE THIRD QUADRANT, i . e . -FREQ AND -DAMPING 
$ 

TABDMP1, 1, ,,,,,,, +TAB1 

+TAB1, -1000., -.04,0., -.04,1., + .06, 1000., + . 06, +TAB2 
+TAB2 , ENDT 
$ 

$ INSERT DMI.DAT HERE 
$ 


DMI* 

DMI0000 



2 

123 1 . 250000000E+00 

DMI 

DMI0000 0 

2 

1 

0 

2624 2624 

DMI* 

DMI0001 



2 

27 1 . 250000000E+00 

DMI* 

DMI0001 



8 

33 5.000000000E+00 

DMI* 

DMI0001 



14 

39 2 . 500000000E+00 

DMI* 

DMI0001 



20 

45 5.000000000E+00 

DMI* 

DMI0001 



26 

51 2 . 500000000E+00 

DMI* 

DMI0001 



32 

57 5 . 000000000E+00 

DMI* 

DMI0001 



38 

. 63 2 . 500000000E+00 

DMI* 

DMI0001 



44 

69 5. OOOOOOOOOE+OO 

DMI* 

DMI0001 



50 

75 2 . 500000000E+00 

DMI* 

DMI0001 



56 

81 5. 000000000E+00 

DMI* 

DMI0001 



62 

87 2 . 500000000E+00 

DMI* 

DMI0001 



68 

93 5. OOOOOOOOOE+OO 

DMI* 

DMI0001 



74 

99 2. 500000000E+00 

DMI* 

DMI0001 



80 

105 5. OOOOOOOOOE+OO 

DMI* 

DMI0001 



86 

111 2 . 500000000E+00 

DMI* 

DMI0001 



92 

117 5. OOOOOOOOOE+OO 

DMI 

DMI 0001 0 

2 

1 

0 

2624 122 

DMI 

SUPER 0 

2 

1 

0 

2 1 

DMI* 

SUPER 



1 

1 0. OOOOOOOOOE+OO 


* 1.000000000E+00 

ENDDATA 


- 153 - 



APPENDIX C3.2 : Run 3.2 Input - Two-Dimensional Fluid/Beam Model 


$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$$$$$ AREA MATRIX PROGRAM - USER INPUT $$$$$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ 

$ TOL1 , TOL2 
1 • / 0. 

$ IFF, FF 

1 / 1 . 

$ SEARCH REGION 
N 

$ SELECT ABSORPTION AREAS 

- 1,1 

$ INPUT THE NUMBER OF FREQ AT WHICH THE 
$ ABSORPTION WILL BE SPECIFIED 
1 

$ INPUT THE FREQUENCIES 

200 . 

$ INPUT THE RHO* SPEED OF SOUND FOR FLUID 
5.645 

$ INPUT THE REAL AND IMAGINARY COMPONENTS 
$ FOR THE SPECIFIC ADMITTANCE - REGION 2 
6.E-1, 1 . 2 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

PARAM, NORF, -1 $ NO RESIDUAL FLEXIBILITY FOR THE FLUID 
$ 

$ INSERT DMI.DAT HERE 

$ 


DMIG 

DAMP101 

0 

6 3 

0 



DMIG* 

DAMP101 



20008 

1 


★ 



20008 

1 

2. 657218874E-01 

5 . 314437747E-01 

DMIG* 

DAMP101 



20009 

1 


★ 



20009 

1 

5.314437747E-01 

1 . 062887549E+00 

DMIG* 

DAMP101 



20010 

1 


★ 



20010 

1 

2 . 65721 8874E-01 

5. 314437747E-01 

DMIG* 

DAMP101 



20011 

1 


★ 



20011 

1 

5 . 314437747E-01 

1 . 062887 54 9E+00 

DMIG* 

DAMP101 



20012 

1 


★ 



20012 

1 

2 . 65721887 4E-01 

5.314437747E-01 

DMIG* 

DAMP101 



20013 

1 


★ 



20013 

1 

5. 31443774 7E-01 

1 . 062887 54 9E+00 

DMIG* 

DAMP101 



20014 

1 


★ 



20014 

1 

2 . 65721887 4E-01 

5.314437747E-01 

DMIG* 

DAMP101 



20015 

1 


★ 



20015 

1 

5.314437747E-01 

1 . 062887549E+00 

DMIG* 

DAMP101 



20016 

1 


★ 



20016 

1 

1 . 328609437E-01 

2. 657218874E-01 

DTI, AERO, 0,1 






DTI, AERO, 1,0, 0,1 

* f 

0, 0, ENDREC 




DTI 

AERO 

2 

1.01256.637 



DMI* 

DMI 0000 



2 

123 

1 . 250000000E+00 

DMI 

DMI0000 

0 

2 1 

0 

2624 2624 

DMI* 

DMI0001 



2 

27 

1 . 250000000E+00 
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DMI* 

DMI0001 


8 

33 

5 . 000000000E+00 

DMI* 

DMI0001 


14 

39 

2 . 500000000E+00 

DMI* 

DMI0001 


20 

45 

5 . 000000000E+00 

DMI* 

DMI0001 


26 

51 

2 . 500000000E+00 

DMI* 

DMI0001 


32 

57 

5 . 000000000E+00 

DMI* 

DMI0001 


38 

63 

2 . 500000000E+00 

DMI* 

DMI0001 


44 

69 

5 . OOOOOOOOOE+OO 

DMI* 

DMI0001 


50 

75 

2 . 500000000E+00 

DMI* 

DMI0001 


56 

81 

5. OOOOOOOOOE+OO 

DMI* 

DMI0001 


62 

87 

2 . 500000000E+00 

DMI* 

DMI0001 


68 

93 

5. OOOOOOOOOE+OO 

DMI* 

DMI0001 


74 

99 

2 . 500000000E+00 

DMI* 

DMI0001 


80 

105 

5. OOOOOOOOOE+OO 

DMI* 

DMI0001 


86 

111 

2. 500000000E+00 

DMI* 

DMI0001 


92 

117 

5. OOOOOOOOOE+OO 

DMI 

DMI0001 0 2 

1 

0 


2624 122 

DMI 

SUPER 0 2 

1 

0 


2 1 

DMI* 

SUPER 


1 

1 

0. OOOOOOOOOE+OO 

* 

$ 

ENDDATA 

1.000000000E+00 
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APPENDIX C3.3 : Run 3.3 Input - Two-Dimensional Fluid/Beam Model 


$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
$$$$$ AREA MATRIX PROGRAM - USER INPUT $$$$$ 
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
$ 

$ TOL1, TOL2 

1.. 0. 

$ IFF, FF 

1 . 1 . 

$ SEARCH REGION 
N 

$ SELECT ABSORPTION AREAS 

- 1,1 

$ INPUT THE NUMBER OF FREQ AT WHICH THE 
$ ABSORPTION WILL BE SPECIFIED 
3 

$ INPUT THE FREQUENCIES 
2. ,3. ,600. 

$ INPUT THE RHO*SPEED OF SOUND FOR FLUID 
5.645 

$ INPUT THE REAL AND IMAGINARY COMPONENTS 
$ FOR THE SPECIFIC ADMITTANCE - REGION 1 

2. E-3, 3.E-3 

3. E-3, 6.E-3 
6.E-1, 1 . 2 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 


PARAM, NORF, -1 $ NO RESIDUAL FLEXIBILITY FOR THE FLUID 
$ 

$ INSERT • DMI . DAT HERE 
$ 


DMIG 

DAMP101 

0 

6 

3 

0 


DMIG 

DAMP102 

0 

6 

3 

0 


DMIG 

DAMP103 

0 

6 

3 

0 


DMIG* 

DAMP101 




20008 1 


★ 



20008 


1 8 . 85739631 5E-04 

1 . 771479263E-03 

DMIG* 

DAMP101 




20009 1 


* 



20009 


1 1 . 771479263E-03 

3. 542958526E-03 

DMIG* 

DAMP101 




20010 1 


★ 



20010 


1 8 . 857396315E-04 

1 . 771479263E-03 

DMIG* 

DAMP101 




20011 1 


* 



20011 


1 1.771479263E-03 

3. 542958526E-03 

DMIG* 

DAMP101 




20012 1 


★ 



20012 


1 8 . 857396315E-04 

1 . 771479263E-03 

DMIG* 

DAMP101 




20013 1 


★ 



20013 


1 1.771479263E-03 

3 . 542958526E-03 

DMIG* 

DAMP101 




20014 1 
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★ 


20014 

1 

8 . 857396315E-04 

1 . 771479263E-03 

DMIG* 

DAMP101 


20015 

1 


* 


20015 

1 

1.771479263E-03 

3. 542958526E-03 

DMIG* 

DAMP101 


20016 

1 


★ 


20016 

1 

4 . 428698157E-04 

8 . 857396315E-04 

DMIG* 

DAMP102 


20008 

1 


★ 


20008 

1 

1 . 328609418E-03 

2. 65721883 6E-0 3 

DMIG* 

DAMP102 


20009 

•1 


* 


20009 

1 

2 . 65721 8836E-03 

5.314437672E-03 

DMIG* 

DAMP102 


20010 

1 


* 


20010 

1 

1 . 328609418E-03 

2 . 6572188 36E-0 3 

DMIG* 

DAMP102 


20011 

1 


* 


20011 

1 

2. 657218836E-03 

5. 314437672E-03 

DMIG* 

DAMP102 


20012 

1 


★ 


20012 

1 

1.328609418E-03 

2. 65721883 6E-0 3 

DMIG* 

DAMP102 


20013 

1 


★ 


20013 

1 

2.657218836E-03 

5. 314437672E-03 

DMIG* 

DAMP102 


20014 

1 


* 


20014 

1 

1.328609418E-03 

2.657218836E-03 

DMIG* 

DAMP102 


20015 

1 


★ 


20015 

1 

2. 657218836E-03 

5.314437672E-03 

DMIG* 

DAMP102 


20016 

1 


★ 


20016 

1 

6. 643047091E-04 

1. 32860941 8E-0 3 

DMIG* 

DAMP103 


20008 

1 


★ 


20008 

1 

2 . 657218 87 4E-01 

5. 314437747E-01 

DMIG* 

DAMP103 


20009 

1 


★ 


20009 

1 

5. 314437747E-01 

1 . 062887549E+00 

DMIG* 

DAMP103 


20010. 

1 


★ 


20010 

1 

2. 657218874E-01 

5. 314437747E-01 

DMIG* 

DAMP103 


20011 

1 


★ 


20011 

1 

5. 314437747E-01 

1 . 0 62887 54 9E+00 

DMIG* 

DAMP103 


20012 

1 


* 


20012 

1 

2. 657218874E-01 

5. 31443774 7E-01 

DMIG* 

DAMP103 


20013 

1 


★ 


20013 

1 

5. 314437747E-01 

1 . 062887549E+00 

DMIG* 

DAMP103 


20014 

1 


★ 


20014 

1 

2. 657218874E-01 

5. 31443774 7E-01 

DMIG* 

DAMP103 


20015 

1 


★ 


20015 

1 

5 . 314437747E-01 

1 . 0 62887 54 9E+00 

DMIG* 

DAMP103 


20016 

1 


★ 


20016 

1 

1 . 328609437E-01 

2 . 657218874E-01 

DTI, AERO, 0,1 
DTI, AERO, 1,0,0,!., 

0, 0 , ENDREC 




DTI 

AERO 2 

1. 

.0 12.566 

1.0 18.850 

1.03769.911 

DMI* 

DMI0000 


2 

123 

1 . 250000000E+00 

DMI 

DMI0000 0 

2 

1 0 

2624 2624 

DMI* 

DMI0001 . 


2 

27 

1 . 250000000E+00 

DMI* 

DMI0001 


8 

33 

5. 000000000E+00 

DMI* 

DMI0001 


14 

39 

2 . 500000000E+00 

DMI* 

DMI0001 


20 

45 

5 . 000000000E+00 

DMI* 

DMI0001 


26 

51 

2 . 500000000E+00 

DMI* 

DMI0001 


32 

57 

5. 000000000E+00 

DMI* 

DMI0001 


38 

63 

2. 500000000E+00 

DMI* 

DMI0001 


44 

69 

5. 000000000E+00 

DMI* 

DMI0001 


50 

75 

2. 500000000E+00 

DMI* 

DMI0001 


56 

81 

5 . 000000000E+00 

DMI* 

DMI0001 


62 

87 

2. 500000000E+00 

DMI* 

DMI0001 


68 

93 

5. 000000000E+00 
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DMI* 

DMI0001 


DMI* 

DMI0001 


DMI* 

DMI0001 


DMI* 

DMI0001 


DMI 

DMI0001 0 2 

1 

DMI 

SUPER 0 2 

1 

DMI* 

SUPER 


* 

$ 

ENDDATA 

l.OOOOOOOOOE+OO 



74 

99 

2 . 500000000E+00 

80 

105 

5 . 000000000E+00 

86 

111 

2 . 500000000E+00 

92 

117 

5.000000000E+00 

0 


2624 122 

0 


2 1 

1 

1 

0.000000000E+00 
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APPENDIX C3.4 : Run 3.4 Input - Two-Dimensional FIuid/Beam Model 


$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
$$$$$ AREA MATRiX PROGRAM - USER INPUT $$$$$ 
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
$ 

$ TOL1 , TOL2 

1.. 0. 

$ IFF, FF 

1 . 1 . 

$ SEARCH REGION 
N 

$ SELECT ABSORPTION AREAS 

- 1,1 

$ INPUT THE NUMBER OF FREQ AT WHICH THE 
$ ABSORPTION WILL BE SPECIFIED 
3 

$ INPUT THE FREQUENCIES 
2. ,3. ,600. 

$ INPUT THE RHO* SPEED OF SOUND FOR FLUID 
5.645 

$ INPUT THE REAL AND IMAGINARY COMPONENTS 
$ FOR THE SPECIFIC ADMITTANCE - REGION 1 

2. E-3, 3.E-3 

3. E-3, 6.E-3 
6.E-1, 1 . 2 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 


NASTRAN DBSET 1= (DB01, DB15) , DBSET 2= (DB01, DB15) $ 

NASTRAN DBSET 15=(DB15) $ 

NASTRAN REAL=0, HICORE=300000 $ 

ID FLUID STRUCTURE 
SOL 71 
TIME 30 
DIAG 8,13 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 


$ $ 

$ USER INSTRUCTIONS: $ 

$ $ 

$ 1. THE USER MUST INSERT THE DMAP.DAT DATA GENERATED BY THE $ 

$ AREA MATRIX FORTRAN PROGRAM IN THE DMAP SEQUENCE. THE $ 

$ LOCATION IS SHOWN BELOW. $ 

$ $ 

$ 2. ALSO THE USER MUST INSERT THE DMI . DAT DATA GENERATED BY $ 

$ THE AREA MATRIX FORTRAN PROGRAM IN THE BULK DATA DECK. $ 

$ THE LOCATION IS SHOW IN THE BULK DATA DECK. $ 

$ $ 

$ 3. IF ABSORPTION DAMPING DATA, DMIG.DAT, GENERATED BY THE $ 
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$ AREA MATRIX PROGRAM EXIST THE USER MUST ALSO INCLUDE 

$ THIS DATA IN THE BULK DATA DECK. THE LOCATION IS SHOW 

$ IN THE BULK DATA DECK. 

$ 


$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ 

$ ***** BEGIN OF THE ALTER PACKAGE FOR FSI ANALYSIS ***** 

$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ 

$ INSERT THE DMAP . DAT HERE 

$ 

$ INSERT THE DMAP ALTER HERE 

$ 

CEND 


TITLE = COUPLED FLUID-STRUCTURE FREQUENCY RESPONSE 
SUBTITLE = 3% BEAM, 2% FLUID DAMPING 
LABEL = LINEARILY VARYING ABSORPTION - NORF=l 
$ 

SET 99 = 1,0 
SET 67 = 21216 
SET 65 = 0 
$ 

SELG = ALL 
SELR = ALL 
SEDR = 65 
$ 

SPC = 1 
VELO = 67 
FREQ = 1 
SDAMPING = 1 
$ 

SUBCASE 1 
SUPER =99 
LOADSET = 11 
DLOAD =21 

$ 

OUTPUT (XYPLOT) 

CSCALE 4.0 
XPAPER = 105. 

YPAPER =80. 

XAXIS = YES 
YAXIS = YES 
XGRIDLINES = YES 
YGRIDLINES = YES 
XMIN =0. 

XMAX = 400. 

YMIN =-30. 

YMAX =30. 


XTITLE = FREQUENCY (HZ) 

YTITLE = PRESSURE (PSI) 

TCURVE = PRESSURE AT GRID 21216 
XYPLOT VELO / 21216 (T1RM) 

BEGIN BULK 

$ 

$ DEFINE PARAMETERS 

$ 

PARAM, NORF, 1 $ RESIDUAL FLEXIBILITY FOR FLUID 
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PARAM, SMALL, 0 . 

PARAM, LMODES, 100 
PARAM, DLOAD, -1 
$ 

$ DEFINE THE DYNAMIC LOADS 

$ 

FREQ1, 1, 5. , 1 . , 395 

$ 

FORCE, 1,10008, ,1000., ,1. 

LSEQ, 11,101,1 

$ 

RLOAD1, 21, 101, , , 1 

$ 

TABLED1 , 1 , , , , , , , ,+TDl 
+TD1, 0 . , 1 . , 1000 . , 1 . , ENDT 
$ 

$ DEFINE MODAL DAMPING FOR FLUID AND BEAM 
$2% FOR FLUID AND 3% FOR THE BEAM 
$ NOTE: THE FLUID MODAL DAMPING IS DEFINED IN 
$ THE THIRD QUADRANT, i . e . -FREQ AND -DAMPING 
$ 

TABDMP1 , 1, , , , , , , ,+TABl 

+TAB1, -1000., -.04,0., -.04,1., +.06, 1000., +. 06, +TAB2 
+TAB2, ENDT 
$ 

$ INSERT DMI.DAT HERE 
$ 


DMIG 

DAMP101 0 

6 

3 

0 


DMIG 

DAMP102 0 

6 

3 

0 


DMIG 

DAMP103 0 

6 

3 

0 



DMIG* 

DAMP101 



20008 

1 


★ 


20008 


1 

8 . 85739631 5E-0 4 

1 . 771479263E-03 

DMIG* 

DAMP101 



20009 

1 


★ 


20009 


1 

1.771479263E-03 

3 . 542958526E-03 

DMIG* 

DAMP101 



20010 

1 


★ 


20010 


1 

8. 85739631 5E-0 4 

1 . 771479263E-03 

DMIG* 

DAMP101 



20011 

1 


* 


20011 


1 

1. 77147926 3E “03 

3 . 542958526E-03 

DMIG* 

DAMP101 



20012 

1 


* 


20012 


1 

8. 85739631 5E-0 4 

1 . 771479263E-03 

DMIG* 

DAMP 101 



20013 

1 


★ 


20013 


1 

1 . 771479263E-03 

3 . 542958526E-03 

DMIG* 

DAMP101 



20014 

1 


★ 


20014 


1 

8. 85739631 5E-0 4 

1 . 771479263E-03 

DMIG* 

DAMP101 



20015 

1 


★ 


20015 


1 

1. 77147926 3E-0 3 

3 . 542958526E-03 

DMIG* 

DAMP101 



20016 

1 


★ 


20016 


1 

4 . 4286981 57E-04 

8 . 857396315E-04 

DMIG* 

DAMP102 



20008 

1 


★ 


20008 


1 

1.328609418E-03 

2 . 6572188 36E-0 3 

DMIG* 

DAMP102 



20009 

1 


★ 


20009 


1 

2. 657218836E-03 

5. 314437672E-03 

DMIG* 

DAMP102 



20010 

1 


★ 


20010 


1 

1. 32860941 8E-03 

2. 657218836E-03 

DMIG* 

DAMP102 



20011 

1 


* 


20011 


1 

2 . 6572188 36E-03 

5. 314437672E-03 

DMIG* 

DAMP102 



20012 

1 


* 


20012 


1 

1.328609418E-03 

2. 657218836E-03 
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DMIG* 

DAMP102 



20013 

1 


★ 



20013 

1 

2.657218836E-03 

5 . 314437672E-03 

DMIG* 

DAMP102 



20014 

1 


★ 



20014 

1 

1 . 328609418E-03 

2 . 657218836E-03 

DMIG* 

DAMP102 



20015 

1 


★ 



20015 

1 

2 . 657218836E-03 

5. 314437672E-03 

DMIG* 

DAMP102 



20016 

1 


* 



20016 

1 

6. 643047091E-04 

1 . 328609418E-03 

DMIG* 

DAMP103 



20008 

1 


★ 



20008 

1 

2. 657218874E-01 

5 . 314437747E-01 

DMIG* 

DAMP103 



20009 

1 


★ 



20009 

1 

5 . 314437747E-01 

1 . 062887549E+00 

DMIG* 

DAMP103 



20010 

1 


★ 



20010 

1 

2. 657218874E-01 

5 . 314437747E-01 

DMIG* 

DAMP103 



20011 

1 


* 



20011 

1 

5. 314437747E-01 

1 . 062887549E+00 

DMIG* 

DAMP103 



20012 

1 


★ 



20012 

1 

2 . 657218874E-01 

5. 314437747E-01 

DMIG* 

DAMP103 



20013 

1 


★ 



20013 

1 

5. 314437747E-01 

1 . 062887 54 9E+00 

DMIG* 

DAMP103 



20014 

1 


* 



20014 

1 

2. 657218874E-01 

5. 31443774 7E-01 

DMIG* 

DAMP103 



20015 

1 


★ 



20015 

1 

5. 314437747E-01 

1 . 062887549E+00 

DMIG* 

DAMP103 



20016 

1 


* 



20016 

1 

1.328609437E-01 

2 . 657218874E-01 

DTI, AERO, 0,1 






DTI, AERO, 1, 0, 0, 1 

• f 

0, 0, ENDREC 




DTI 

AERO 

2 

1.0 

12.566 

1.0 18.850 

1.03769.911 

DMI* 

DMI0000 



2 

123 

1 . 250000000E+00 

DMI 

DMI0000 

0 

2 1 

i 

0 

2624 2624 

DMI* 

DMI0001 



2 

27 

1.250000000E+00 

DMI* 

DMI0001 



8 

33 

5. 000000000E+00 

DMI* 

DMI0001 



14 

39 

2 . 500000000E+00 

DMI* 

DMI0001 



20 

45 

5. 000000000E+00 

DMI* 

DMI0001 



26 

51 

2 . 500000000E+00 

DMI* 

DMI0001 



32 

57 

5. 000000000E+00 

DMI* 

DMI0001 



38 

63 

2 . 500000000E+00 

DMI* 

DMI0001 



44 

69 

5. 000000000E+00 

DMI* 

DMI0001 



50 

75 

2 . 500000000E+00 

DMI* 

DMI0001 



56 

81 

5. 000000000E+00 

DMI* 

DMI0001 



62 

87 

2 . 500000000E+00 

DMI* 

DMI0001 



68 

93 

5. 000000000E+00 

DMI* 

DMI 0001 



74 

99 

2 . 500000000E+00 

DMI* 

DMI 0001 



80 

105 

5. 000000000E+00 

DMI* 

DMI0001 



86 

111 

2 . 500000000E+00 

DMI* 

DMI0001 



92 

117 

5. 000000000E+00 

DMI 

DMI0001 

0 

2 1 

0 

2624 122 

DMI 

SUPER 

0 

2 1 

0 

2 1 

DMI* 

SUPER 



1 

1 

0. 000000000E+00 


* 1 . OOOOOOOOOE+OO 

$ 

ENDDATA 
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APPENDIX C4 : Run 4 Input - Two-Dimensional Fluid/Beam Model 


NASTRAN DBSET 1= (DB01, DB14, DB15) , DBSET 2= (DB01, DB14, DB15) $ 
NASTRAN DBSET 15= (DB14, DB15) $ 

NASTRAN REAL=0,HICORE=300000 $ 

ID FLUID STRUCTURE 
TIME 30 
DIAG 8,13 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 


$ $ 

$ USER INSTRUCTIONS: $ 

$ $ 

$ CASE CONTROL DECK - $ 

$ $ 

$ SET i = FLUID GRID ID’s THAT AREA ON THE SURFACE OR $ 

$ = ALL $ 

$ DISP = i $ 

$ $ 

$ SET j = SPECIFY FREQUENCIES AT WHICH THE GRID POINT $ 

$ PARTICIPATION FACTORS ARE TO BE CALCULATED $ 

$ OFREQ = j $ 

$ $ 

$ BULK DATA DECK - $ 

$ $ 

$ DMIG,PDOF,0, 6,1,0 $ 

$ DMIG, PDOF, GIDi, CID, , GIDi, CID, 1 . 0 $ 

$ $ 

$ GID IS THE GRID POINT FOR WHICH THE PARTICIPATION FACTORS $ 

$ ARE TO BE CALCULATED $ 

$ CID IS THE COMPONENT NUMBER $ 

$ $ 

$ IF ONE WANTS ONLY THE REAL DISPLACEMENT OUTPUT SET THE $ 

$ PARAM, REAL, 0 $ 

$ $ 


$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ 

$ INSERT THE DMAP PROGRAM HERE 

$ 

CEND 

TITLE = GRID POINT PARTICIPATION FACTOR CALCULATION 
SUBTITLE = FOR FLUID GRID 21216 
$ 

$ INSERT SET. DAT HERE 

$ 

SET 911 = 20000, 20001, 20002, 20003, 20004, 20005, 20006, 

20007, 20008, 20009, 20010, 20011, 20012, 20013, 

20014 20015 20016 

DISP = 911- 

SET 66 = 40., 270. 

OFREQ =66 
SUBCASE 11 

LABEL = LOAD ON THE STRUCTURE (SUPERELEMENT 1) 
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BEGIN BULK 

P ARAM, EQUAL, 0 

DMIG, PDOF, 0, 6, 1, 0 

DMIG, PDOF, 21216, 1, , 2121 6,1,1. 

ENDDATA 



APPENDIX D : INPUT DECKS FOR SECTION 8 
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APPENDIX D1 


Run 1 Input - Vehicle Acoustic Analysis 


** R U N 1 ** 


ID FLUID STRUCTURE 
SOL 63 
TIME 500 
DIAG 8,13 

READ 9 $ FSI.63 ! ! ! NOTE ! ! ! THIS MUST BE IMMEDIATELY BEFORE CEND 
CEND 

TITLE = FLUID-STRUCTURE INTERACTION 
SUBTITLE = NORMAL MODES ANALYSIS 
SEALL = ALL 
SPC = 111 
$ 

$ DEFINE ALL THE SUPERELEMENT ID'S EXCEPT 0 

$ 

SET 65 = 

DISP (PLOT) = ALL 
SUBCASE 1 
SUPER = 1 
METHOD = 1 

$ 

$ DEFINE THE PARTITION VECTOR FOR THE FLUID 

$ 

SET 66 = 100000 THRU 101000 

$ 

SUBCASE 11 

LABEL = RESIDUAL STRUCTURE MODEL 
LOADSET =11 
PARAM, COUPMASS, 1 
PARTN = 66 
METHOD =11 
SUBCASE 12 

LABEL = FLUID MODEL 
METHOD =12 
BEGIN BULK 
$ 

$ DEFINE THE STRUCTURE MODEL 

$ 

$ 

$ DEFINE THE FLUID MODEL 

$ 

LSEQ, 11, 101,1 
LSEQ, 11, 102, 2 
LSEQ, 11, 103, 3 
LSEQ, 11, 104, 4 
$ 

$ DEFINE THE EIGENVALUE EXTRACTION DATA FOR THE 
$ FLUID 


^6* Ji^JNENIIONAIJS 8DWK 
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$ 

EIGRL, 12,-1.,, 10,1 

SPC1 111 23456 100000 THRU 101000 

$ MN12 EP2 ACOUSTIC MODEL UPDATE 2 
PSOLID 100001 100001 2 

PSOLID 100002 100002 2 

PSOLID 100003 100002 2 


MAT 9 

100001 

8317. E8 



+MMT1A 

+MMT1A 

+MMT1B 

+MMT1B 

8317. E8 


8317. E8 

7.058+3 

MAT 9 

100002 

7994. E7 



+MMT2A 

+MMT2A 

+MMT2B 

+MMT2B 

7994. E7 


7994. E7 

7.058+3 

GRID 

100001 


2023.15 443.730 519.240 



GRID 

100947 


2912.00 

400.000 

460.590 



CPENTA 

100001 

100001 

100207 

100212 

100508 

100203 

100216 

+EACU674 

100109 

100120 






CHEXA 

100675 

100001 

100107 

100102 

100947 

100946 

100109 

100104+EACU675 







+EACU675 

100122 

100120 






CPENTA 

100676 

100001 

100947 

100102 

100103 

100122 

100104 


$ 

$ 

$ DEFINE THE PRESSURE LOAD ON THE FLUID TO 
$ CALCULATE THE AREAS 
$ 

P LOAD 4 1 100001 1. 100212 


PLOAD4 1 100676 1. ' 100122 

$ 

$ THIS PRESSURE LOAD DEFINES THE ABSORPTION AREA FOR THE ROOF 
$ 

PLOAD4 2 100002 1. 100207 


PLOAD4 2 100569 1. 100236 

$ 

$ THIS PRESSURE LOAD DEFINES THE ABSORPTION AREA FOR THE FLOOR 
$ 

PLOAD4 3 100084 1. 100302 


P LOAD 4 3 100676 1. 100103 

$ 
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100585 


100105 


100203 


100190 


100413 


100299 


100104 



$ THIS PRESSURE LOAD DEFINES THE ABSORPTION AREA FOR THE TRUNK LID 


$ 

PLOAD4 

4 

100041 

1 . 

100353 

100555 

PLOAD4 

4 

100605 

1 . 

100465 

100390 

ENDDATA 
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APPENDIX D3.1 : Run 3.1 Input - Vehicle Acoustic Analysis 


★★★★★★★★★★★★★★★A 

** R U N 3.1 ** 


NASTRAN DBSET 1= (DB01 , DB02 ) , DBSET 2= (DB01, DB02) $ 

NASTRAN DBSET 15=(DB02) $ 

BUFFSIZE=9217 $ 

ID ACUSTIC, TWO 
TIME 1000 
DIAG 8 
$ 

$ *** INSERT SOL 71 DMAP ALTER HERE 

$ 

READ 9 $ FSI.71 ! ! ! NOTE ! ! ! THIS MUST BE IMMEDIATELY BEFORE CEND 
CEND 

TITLE = FREQUENCY RESPONSE CALCULATION 
SUBTITLE = 3% MODAL DAMPING 
$ 

SET 99 = 0 
SELG =99 
SEDR =99 
$ 

SET 66 = 100597 
VELO (PHASE) = 66 
$ 

FREQ = 1 
SDAMP = 1 
$ 

SUBCASE 11 

LABEL = DRIVE SHAFT UNBALANCE 
LOADSET =33 
DLOAD = 201 
OUTPUT (XYPLOT) 

SEPLOT 0 
CSCALE 4.0 
XPAPER=105 . 

YPAPER=8 0 . 

XAXIS=YES 

YAXIS=YES 

XGRID LINES=YES 

YGRID LINES=YES 

YLOG = YES 

XDIVISIONS = 9 

XTITLE = FREQUENCY (HZ) 

YTITLE = PRESSURE 

TCURVE = PRESSURE AT DRIVERS EAR (GRID 100597) 

XYPLOT VELO / 100597 (T1RM) 

BEGIN BULK 

$ 

$ DEFINE PARAMETERS 

$ 

PARAM, SMALL, 0 . 


- 173 - 

*ank pREC£D1NG page blank mot filmed 


PARAM, LMODES, 2000 
PARAM, DLOAD, -1 
PARAM, RESDUAL, -1 
$ 

$ DEFINE LOADS 

$ 

FREQ1, 1, 10 . , 1 . , 90 

$ 

DAREA,201, 32818, 3,1. 

DAREA, 202, 9261, 3,1. 

DAREA,203, 35174, 3,1. 

DAREA, 204, 11927, 3,1. 

$ 

RLOAD1, 21, 201, , , 1 
RLOAD1, 22, 202, , , 1 
RLOAD1, 23, 203, , , 1 
RLOAD1, 24, 204, , , 1 
$ 

DLOAD, 201, 1 . , 1 . , 21 . 1 . , 22 . 1 . , 23 
+,1.,24 

DLOAD, 202, 1 . , -1 . , 21 . 1 . , 22, 1 . , 23 
24 

TABLED 1, 1 

+,0., 3. 536+10, 1000., 3. 536+10, ENDT 

$ 

$ DEFINE DAMPING 

$ 

TABDMP1, 1 

+,-1000., -.06,0., -.06, 1.-6, +0.06, 1000., +.06 
+,ENDT 
$ 

$ INSERT DMI.DAT HERE 

$ 

include (DMI) 

ENDDATA 
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APPENDIX D4 


Run 4 Input - Vehicle Acoustic Analysis 


★★★★★★★★★★★★★★ 

** R U N 4 ** 
★★★★★★★★★★★★★★ 


NASTRAN DBSET 1= (DB01, DB02, DB15) , DBSET 2= (DB01, DB02, DB15) $ 

NASTRAN DBSET 15= (DB02, DB15) $ 

NASTRAN BUFFSIZE=9217 $ 

ID FLUID STRUCTURE 
TIME 30 
DIAG 8,13 

READ 9 $ FSI . 71R ! ! ! NOTE ! ! ! THIS MUST BE IMMEDIATELY BEFORE CEND 
CEND 

TITLE = AUTOMOBILE ACOUSTIC ANALYSIS - MODAL DAMPING ONLY 
SUBTITLE = PARTICIPATION FACTOR CALCULATION FOR DRIVERS EAR 
$ 

$ INSERT SET. DAT HERE 

$ 

include (SET) 

DISP = 911 

SET 912 = 64., 86. 

OFREQ = 912 
SUBCASE 11 

LABEL = DRIVE SHAFT UNBALANCE 
BEGIN BULK 
DMIG, PDOF, 0, 6,1, 0 
DMIG, PDOF, 100597, 1, , 100597, 1, 1 . 

ENDDATA 
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